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ABSTRACT 


HIGH ENERGY SLIT APERTURE SPECT AND SIMPLIFIED INVITRO 
METHODS FOR THE DOSIMETRY OF POSITRON EMITTING RADIOTRACERS 

by 

Mark Christopher Wrobel 
Captain, USAF 
1996 

198 pages 

PhD, Environmental Health Sciences 
University of Michagan 

Co-Chairs: James E. Martin, W. Leslie Rogers 

The dosimetry of new positron emitting radiopharmaceuticals is initially estimated 
using animal tissue and organ biodistributions assessed invitro. Such methods are time 
and labor intensive and can have limited accuracy. This research investigated two 
alternative methods by which biodistribution can be obtained from the laboratory rat: 
invitro organ measurements using a reduced sacrifice technique, and invivo 
measurements using single photon emission computed tomography (SPECT). 

As an alternative to a four time-point sacrifice method, a two time-point method 
was evaluated as a means to determine the organ cumulated activity of * ’C labeled 
radiopharmaceuticals. Residence times calculated using two time-points acquired during 
the first half-life of **C were either equivalent or larger than those resulting from using 
four sacrifice times. Correction factors were required for the urinary bladder and 
gallbladder when using this simplified technique due to delayed uptakes. 

Invivo assessments were performed with SPRINT, a full ring detector SPECT 
system using a slit aperture to obtain a 3-to-l object to image magnification ratio. 



Acceptable resolution for 511 keV photons was achieved using a high energy parallel 
slice collimator and a novel technique to correct for penetration of the slit aperture by 
high energy photons. The resulting system resolution was approximately 4.5 mm axially 
and transaxially. System sensitivity was ~55 cpm/pCi, a consequence of high resolution 
collimation and poor intrinsic detector efficiency. 

Qualitative and quantitative high energy imaging was performed using two- 
dimensional and three-dimensional phantoms and live rats injected with 
[**F]fluorodeoxyglucose. Phantoms displayed well resolved hot and cold objects, and 
animal images showed good contrast of the brain and heart regions. Hot phantom regions 
and animal organs were quantified with an accuracy of ±8% of the known activity, with a 
standard deviation of ±15%. Accurate determination of region activity required 
substantial activity concentrations (~100 pCi/ml) and imaging times of one hour. 

Analysis of animal biodistributions using a reduced sacrifice protocol provides 
results in good agreement with and generally conservative to results using all sacrifice 
intervals. Animal biodistributions assessed invivo with slit aperture SPECT can show 
good accuracy for organs with substantial uptakes, but the method is hindered by very 
poor system sensitivity. 
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CHAPTER I 


INTRODUCTION 


Problem Statement and Research Objectives 

Absorbed dose calculations are required to evaluate the risks posed by the use of 
new positron emission tomography (PET) radiotracers in medical studies. Human 
absorbed doses are initially estimated using animal tissue and organ biodistributions 
assessed invitro and invivo’. The focus of this research is to investigate two methods by 
which biodistribution information can be obtained from the laboratory rat: invitro organ 
measurements using a reduced sacrifice technique and invivo organ measurement using 
single photon emission computed tomography (SPECT). 

Determining the biodistribution of PET radiotracers using conventional tissue 
counting techniques is labor intensive and requires the sacrifice of numerous animals. 
Conventional methods require four or more temporal measurements out to 3 or 4 half- 
lives of *'C, with four to six animals required per sampling interval. Animal sacrifice 
data for several '*C labeled radiopharmaceuticals such as epinephrine, raclopride, and 2a- 
tropanyl benzilate ([1 lC]TRB) suggest that only one or two sacrifice intervals soon after 
injection could provide a preliminary assessment of organ cumulated activity that is in 
good agreement with conventional multiple time sacrifice methods. 

Invivo assessment of radionuclide organ activities via single photon emission 
computed tomography (SPECT) is a desirable alternative to conventional sacrifice 
methods if it can be shown to accurately determine organ specific uptake and clearance. 
The advantage of such a method is the ability to trace the biodistribution pattern over time 
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within the same animal. Further, invivo assessment should prove more accurate than 
conventional methods for organs which have complex biodistributions such as the urinary 
bladder and gallbladder. This work proposed to determine the efficacy of assessing 
animal biodistributions of 511 keV emitting radiotracers using a newly developed SPECT 
system^ as an alternative to conventional animal sacrifice methods. 

Objectives 

1) To evaluate for **C labeled radiopharmaceuticals whether a limited sacrifice 
method can be used as a simple but efficient estimator of human organ doses for short¬ 
lived PET agents. 

2) To evaluate the potential and limitations of using a full-ring single photon 
emission computed tomography system to non-invasively determine the biokinetics of 
positron emitting labeled radiopharmaceuticals in the rat. 

Background 

In the development of new radiotracers, a single site in the body, such as the brain, 
and a given metabolic process, such as glucose metabolism, are of scientific or medical 
interest to the investigator. However, radiation doses are received from all tissues of the 
body following the systemic introduction of a radiotracer. Estimation of the absorbed 
dose delivered to particular organs and tissues for a radiotracer is needed to determine the 
effective dose to the patient, establish the optimal diagnostic regimen to minimize dose, 
establish limits on administered activity so that regulatory limits are not exceeded and, 
as far as practical, to compare risks attributable to the radiation with the possible benefits 
of the investigation. In order to determine the dose received by a given tissue or organ, 
the complex relationship between the physical characteristics of the radionuclide used and 
the biological behavior of the ligand and its metabolites must be known. 
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In general, the dosimetry of a radiopharmaceutical requires the following types of 
information: 

(a) the activity administered 

(b) the physical half-life of the radionuclide, 

(c) the type, energy and frequency of particulate and electromagnetic emissions; 

(d) the identity and hence location, mass and shape of target organs, 

(e) the temporal and spatial distribution of activity in source regions, 

(f) the fraction of radiation energy emitted in the source region which is absorbed 
in the target region. 

The first three involve only physical data, and are generally well-known. The 
latter three however require biological information, are often difficult to establish and can 
show significant variation between species and from subject to subject. In particular, 
measuring the temporal biodistribution of a new radiotracer is the greatest challenge in 
assessing dosimetry since data for d) and f) are available for a variety of “standard 
human” geometries'*. Since radiation dosimetry requires biodistribution information, 
investigators must use animals, particular the laboratory rat, to establish approximate yet 
conservative measures of the activity present in various tissues as a function of time post 
injection. Historically, animal data has included total body radioactivity measurements, 
the assay of specimens of blood and excreta, and the invitro analysis of whole organs 
following sacrifice. Within the last several years, pinhole collimated SPECT and high 
resolution PET have also been proposed as a means for invivo quantification of organ 
activity in small laboratory animals 5 . 6 , 7 , 8 , 9 , 10 . 11,12 ^ focus of this work will be to evaluate 
the limitations and benefits of pinhole imaging for high energy gamma-rays. From invivo 
or invitro serial measurements a “model” for the variation with time of activity in various 
organs can be constructed. This information then allows estimation of human dose using 
the Medical Internal Radiation Dosimetry formalism. 
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Physical Concepts and the MIRD Formalism 


The Medical Internal Radiation Dosimetiy (MIRD) schema has been broadly 
adopted for the assessment of internal dose resulting from the internal application of 
radioactivity in medicine*. The MIRD method has been described fully by Loevinger and 
Berman*^. The fundamental concept of this schema is that absorbed dose calculations are 
performed for pairs of regions, one of which is called the source region, and the other the 
target region. The activity in the source region emits radiation which imparts energy to 
the target region. Source and target regions can be equal, overlap, separate, or one region 
can entirely encompass another region. Certain quantities depend on interactions between 
a source and target region, and to distinguish between the two, the dependence is shown 
as (rk<—rh), where rk is the target and rh is the source. 

The absorbed fraction in a volume from a region is defined as the energy 
imparted to a target volume v from source region r divided by the energy emitted by the 

source region: 

^(v <— r) 

The absorbed fraction is defined only for target regions that are volumes. The 
specific absorbedfraction, O, is defined for a target volume v as the quotient of the 

absorbed fraction and the mass of v: 

^(v <— r) 


0 (v <-r) = - 


m„ 


In general, the specific absorbed fraction in any target region rk from any source 
region rn is defined as: 

In the MIRD schema, the model is used to designate the assumptions used in the 
dose calculation. This typically will include two specific sub-models: an anatomic 
model which provides in formation on the spatial and structural assumption of the human. 
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and a kinetic model which includes assumptions relative to the biodistribution including 
uptake and clearance of an administered substance. The absorbed dose of interest is the 
energy imparted per unit mass to specified regions of the anatomic model. The regions of 
the anatomic model are assumed to be uniform and homogenous, unless macroscopic 


variations in density or structure are specifically included. 

Characterization of a source region requires knowledge of the decay mechanisms 
of the radionuclide. The number of decays which a radionuclide undergoes in a given 
source region is designated cumulated activity and is normally expressed in microcurie- 
hour (pCi-hr), which equals 1.332x10* nuclear transformations. The unit of energy is the 
rad-gram (rad-g), which equals 6.242x10’ MeV. In general, the mean energy emitted per 
unit cumulated activity for type i radiation is: 


V //Ci - hry 


2.13i7,£, 


where nj is the number of particles of type i emitted per transformation and Ej is the 
energy of the particle in MeV. 


MIRD Equations for the General Case 

Consider a pair of source and target regions. Energy is emitted from source region 
rh and some fraction of the energy is absorbed by a target region rk. Let the total activity 
in rh equal Ah, the mean energy emitted per nuclear transformation be A, and the specific 
absorbed fraction be 0(rk'<-rh). Then the rate of energy absorption per unit mass of the 
target region is the mean absorbed dose rate: 

There is generally more than one type of radiation associated with a radionuclide, 
as is the case for positron emitters, and the above equation can be expanded to a more 
general form to account for each radiation emitted: 
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i 

If this expression is integrated over some time interval of interest, the total 
absorbed dose to a given target is obtained. It is generally assumed that Oi is independent 
of time, so the equation for mean dose becomes: 

D{rt <-fh)= A (h/2)E ^ h) 

i 

where A h(ti,t 2 ) is the integral of Ah(t) from ti to t 2 . For a given radionuclide and a 
given source target geometry, the quantities Ai and Oi can be uniquely determined, and a 
quantity S can be defined as: 

% ^ 'i.)=Z (^k <- h) 

i 

The quantity S is designated the S-value and equals the mean absorbed dose per 
unit cumulated activity, with typical units of rad/pCi-hr. In actual application, there can 
be many source regions irradiating a given target. The total mean absorbed dose to a 
target region rk is: 

h 

h 

For application to models in medicine and biology, it is frequently assumed that 
the relative spatial distribution of activity within a given source region (i.e. a tissue or 
organ) is constant over time, while the total activity in a source region is expected to vary 
with time. Under these conditions, the quantities which relate the energy emitted by a 
given source to the energy absorbed by a given target are independent of time, and the 
S-value is a constant. 

When the assumption is made that the source and target organs are in a large 
homogeneous absorbing material where edge effects can be ignored and activity is 
uniformly distributed in source regions, the criteria for the uniform isotropic model are 
obtained. Under such conditions, a reciprocal relationship exists between any pair of 
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source-target organs: the specific absorbed fraction O and the mean absorbed dose per 
unit cumulated activity S are independent of which region is designated the target or 
source. 

The uniform isotropic model is used in this research. S-values are thus constant 
and have been extracted for reference man'* from the MIRD-III program*'*. Source 
regions will typically be either whole organs or whole organ contents for which activity 
content must be determined as a funetion of time. 

Outline of Dissertation 

This dissertation describes research on methods of assessing radiotracer 
biodistributions in the laboratory rat. Chapter II discusses the methods and findings 
associated with a reduced sacrifice technique for determining rat biodistributions and 
compares the results with conventional multiple time-point methods. In Chapter III, the 
concepts of SPECT are introduced and the salient characteristics and expected 
performance of SPRINT are described. Chapter IV examines conventional SPECT image 
reconstruction using direct and iterative methods. The conventional filtered back- 
projection for parallel and fan-beam geometries is presented, and statistical iterative 
reconstruction techniques are examined with focus on the space-alternating generalized 
expectation-maximization (SAGE) algorithm. 

Chapter V focuses on the slice collimator requirements of SPRINT, with 
emphasis on high energy imaging. This chapter discusses those design, construction 
methods, and performance characteristics of parallel and fan-beam slice collimation that 
are necessary for achieving a desired axial resolution with SPRINT. In Chapter VI, the 
issues of quantitative SPECT imaging using SPRINT are introduced, with focus on new 
research conducted to correct for aperture penetration of high energy gamma rays. 

Chapter VO discusses other factors considered for conducting high energy imaging using 
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SPRINT and then presents qualitative and quantitative results of 511 keV imaging of 2-D 
and 3-D phantoms and live rats injected with [**F]fluorodeoxyglucose. Finally, Chapter 
VIII presents summary observations and identifies areas of future work. 
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CHAPTER II 


A PROPOSED REDUCED SACRIFICE TECHNIQUE FOR THE DOSIMETRY 
OF ^^C LABELED RADIOPHARMACEUTICALS 

Human absorbed dose from carbon-11 ('*C) labeled radiopharmaceuticals is 
initially estimated using animal tissue and organ biodistributions’, typically in the rat. 

For a new ”C labeled radiotracer, this determination will usually include measuring 
organ activities at four or five time points during the first three or four physical half-lives 
of ”C, using four to six animals per time-point. Organ time-radioactivity profiles are 
obtained from the measured organ radioactivity, and cumulated activity, residence time, 
and organ absorbed dose per unit administered radioactivity are calculated. 

Determination of animal radioactivity biodistributions for ”C labeled radiotracers 
is not trivial. The short 20 minute half-life of ”C requires a rapid sacrifice protocol and 
the use of a significant number of animals to provide high precision and accuracy in the 
determination. These issues have been recognized by Gatley^, who has recommended 
eliminating animal measurements and basing dose estimates on human cardiac output 
models, and Schaumann^ who has recommended limiting sacrifice data for long lived 
’'’C-labeled agents. One of the objectives of this study was to evaluate the use of a 
simplified technique which employs only two sacrifice intervals to assess residence time. 
Animal biodistribution data were obtained for ten ”C labeled radiopharmaceuticals, and 
dosimetry estimates were compared between using all sacrifice intervals and a simplified 
technique. Such an alteration yielded acceptable estimates of human dosimetry. 

The accurate assessment of urinary bladder wall dose and gallbladder wall dose 
are often overlooked when evaluating new short-lived radiotracers. This is partly due to 
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the difficulty of measuring the residence time for these organs as well as the perception 
that because of delayed filling and the short physical half-life of they will be of lesser 
dosimetric importance than organs which receive first-pass depositions fi'om the blood 
pool. The urinary bladder contents and estimated gallbladder contents residence times in 
the rat for ten ^’C radiotracers were estimated using measured organ activities and a set of 
conservative assumptions based on ICRP-53 recommendations'^. Results indicate that the 
bladder organs can be critical in determining the limiting dose of new labeled 
radiotracers for human patients. 


Theory 

Given a single bolus injection of radioactivity, Ao, into the blood pool, the 
radioactivity in organ or tissue S is modeled using the following expression: 

A n+m ^ \ ^ . 

E ^ 7 Z 3 [exp- + ^p)t - exp- {Aj + Ap)t\ \ 

^0 j=n+\ (=1 [ ~ J 

where Fg is the fractional distribution to organ or tissue S, aj is the fraction of Fg 
eliminated with biological removal constant aj is the fraction of Fg taken up with 
biological uptake constant Aj; n is the number of elimination components; m is the 
number of uptake components and Ap is the physical decay constant^. This behavior is 
diagramatically presented in figure 2.1. A bolus of activity Ao is injected into the transfer 
compartment, equivalent to the blood pool or systemic circulation in which instantaneous 
uniform mixing assumed. Pathways are then one-way determined by time-invariant 
transfer constants. 

For absorbed dose calculations in nuclear medicine, and for very short half-life 
tracers in particular (T m < 20 minutes), several simplifying assumptions to the above 
model may be applied. First, uptake is either based on a single component or 
instantaneous uptake is assumed'^. Further, the effective removal of radioactivity from an 
organ is described by a single component, and for very short half-life radiotracers can be 
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taken to be equal to the physical decay constant"^. With single uptake and excretion 
components, the above model becomes: 

= .. JJ T [exp- (/I, + Ap )^ - exp- (A j + )?] 

Aq {Aj A^) 

The rat organ time-radioactivity profiles for the ’’C-labeled compounds listed in 
table 2.1 demonstrate that the maximum organ radioactivity is found at the time of first 
sacrifice, which is either two or five minutes post injection. Radioactivity then decreases 
according to either a mono-exponential or bi-exponential curve with inflection 
consistently occurring after the second sacrifice interval, 10 to 15 minutes post injection. 
The trend in organ activity during the first 15 minutes post injection indicates either 
initial drug washout with an effective removal constant (X, + Xp) being greater than Xp, or 
delayed organ uptake with an effective removal constant less than Xp but greater than 0. 
Sampling points after this period provide the “tail” of the biodistribution to four or five 
physical half-lives of the radiopharmaceutical. Exceptions to this finding are the urinary 
bladder and gastrointestinal track which have complex profiles because of delayed filling 
and emptying. Removal of radioactivity in the “tail” region was observed to be mono¬ 
exponential and either predominantly or solely due to physical decay. 

These findings suggest that the above model may be revised to use two discrete 
time intervals. The first interval is characterized by single uptake and clearance 
components with physical decay and is identical to the original model. The second 
interval is characterized by a single clearance component weighted by physical decay,with 
Xk < Xp (figure 2.2). A suitable model for first compartment source organs with a bolus 
injection of a radiotracer then becomes: 

4(0 = 4^. X ^ [exp-(A, +A^)t- exp -(Aj + )t] t < 15 min 

4(0 = 4(l5)exp-(A^+Ap)t t>15min 
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Table 2.1. PET Radiotracers Evaluated 


Name or Abbreviation 

Chemicai Name and Behavior 

nrqEpj 

[^''Cjepinephrine (hydrophilic) 

["CJmHED 

["CJmetahydroxyephedrine (hydrophilic) 

["C]PHE 

["CJphenylephrine (hydrophilic) 

[”C]MTBZ 

[’^CJmethoxytetrabenazine (lipophilic, amine (3°)) 

[”C]NMPB 

["C]N-methyl piperidyl benzilate, (lipophilic, amine (3°)) 

["C]TBZ 

f ^C]tetrabenazine (lipophilic, amine (3°)) 

["CJTRB 

[”C](+)-2a-tropanyl benzilate (lipophilic, amine (3°)) 

["CJFNZPAM 

["CJflunitrazepam (lipophilic, amide) 

["C]PK11195 

f^C]1-(2-chlorophenyl-N-methyl-N-(1-methylpropyl)-3- 
isoquinoline carboximide (lipophilic, amide) 

I”C]RAC 

["CJraclopride (variable behavior) 


If this model is valid, and assuming A,k is small, two sacrifice intervals should 
then be sufficient to describe the biodistribution. This expected behavior allows us to 
present the hypothesis that differences will be small between biodistributions assessed 
from multiple sacrifice data and that from a simplified method using two sacrifice times, 
one very early post injection and the second after one physical half-life of 


Organ S 



Figure 2.1. Two-Compartment Kinetics Model for Organ Radioactivity 

Determination. 






Injected 

Activity 

Ao 


15 


Organ S 
t<15 min 



Figure 2.2. Simplified Model Proposed for Labeled Radiotracers. 


Organ S Bladder 

Kidney Urinary Bladder 

Liver Gallbladder 



Figure 2.3. Three-Compartment Kinetics Model for Urinary Bladder and 
Gallbladder Radioactivity Determination. 


Urinary Bladder and Gallbladder 

Simplified models have been developed to predict radioactivity in the urinary 
bladder"^ ^ ^ and gallbladder ^ ^ ^ using two, three or more catenary compartments 
(figure 2.3). The model demonstrates a delayed uptake path to the bladder organ, which 
then empties in discrete time intervals. The ICRP model assumes a 3.5 hour voiding 
interval for the urinary bladder and 3 hours for the gallbladder. Since these times are long 
compared to the physical half-life of ‘*C, biological removal is often ignored. However, 
such models generally over-simplify the physiological processes involved and the 
complex relationship between urine or bile flow rate, emptying periods, and the volume 
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present in the bladder which actually occurs . A three compartment model implies 
that sampling be performed at a minimum of six time-points, assuming that measured 
time-radioactivity profiles can be simplified to only three transfer coefficients. In 
practice, four to five time-points have been used to estimate bladder residence times, with 
data reflecting significant variability between both animals at a given time point, as well 
as between compounds. Under such conditions, the simplifying assumptions made for 
solid first-pass organs are not valid for the bladder organs. 

Methods 

To evaluate whether a reduced sacrifice protocol is sufficient to estimate 
cumulated activity, the ten ''C labeled radiotracers in table 2.1 were examined. 
Biodistributions were evaluated at four sacrifice times with four to five Sprague-Dawley 
rats (mixed sex) sacrificed per interval. Animals were injected IV (femoral vein) with up 
to 600 pCi of the subject radiotracer, prepared by the University of Michigan PET 
Facility. Subgroups were then sacrificed at of 2 to 5, 10 to 15, 30 to 45 and 60 to 90 
minutes post injection and whole organs dissected out, weighed, and assessed for in a 
sodium iodide well counter. These intervals will be referred to as Tl, T2, T3 and T4, 
respectively. Organ cumulated activities and residence times were calculated from this 
data, and a statistical analysis conducted to determine standard deviations and confidence 
intervals. The results provided a reference for comparison with the proposed simplified 
technique. The reduced sacrifice method used the above data, but only one or two 
sacrifice intervals were employed to characterize each radiotracer and organ. Cumulated 
activity and residence times assessed from the simplified method were compared using a 
Student’s t-test and the percent difference from the conventional multiple point sacrifice 
method. 



17 


Five permutations of available data were tested. Three permutations used two 
sacrifice intervals: T1 and T2, T1 and T4, and T2 and T4, while the remaining two used a 
single sacrifice interval: either T1 or T2. Cumulated activities assessed from these 
simplified methods were compared using a t-test and the percent difference from the 
extensive multiple point sacrifice method. 

The assessment of cumulated activity was determined using the method of Smith, 
Brownell and Ellett*^ modified for direct integration to determine the cumulated activity 
between each available data point. This point-to-point integration method was deemed 
appropriate because of the limited tissue distribution data available. A limitation of this 
method is that it does not address the radioactivity profile for time before the first 
sacrifice interval, thus excluding a possibly important fi'action of the total cumulated 
activity. The most conservative estimate of this profile is obtained by assuming an 
exponential curve extrapolated back from the first two available data points. This curve 
will be characterized by a y-intercept which gives a conservative value for the organ 
radioactivity at time equals zero, and an effective removal constant equal to the slope of 
the curve. In addition, the method excludes consideration of any natural trend in the 
time-radioactivity profile, such as being predominantly exponential. For the specific case 
of carbon-11 compounds where physical decay dominates the time-radioactivity profile, 
an exponential fit between sacrifice intervals is the most reasonable. 

For the animal activities at each sacrifice interval, the sample average and sample 
standard deviation were calculated. These were then used in the following expression to 
determine the cumulated activity in the organ from the time of injection to infinity: 
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A(O^oo) = 


‘'eff(7i->7-2) 


"[l ^eff(ri-»r2)^2)] 


(A -a{-^^ 


+(A,-A,) 


( 


-ln(4X 

(r.-r,) ^ 


Un(^,)-to(^,); 

+1.44 xTx A, 


The first term on the right of this expression represents the cumulated activity 
from the time of injection to the second sacrifice interval. This assumes an instantaneous 
uptake in the organ of Aq equal to the y-intercept of an exponential fit of the first two 
sacrifice intervals, and a slope of Xefr. The next two terms assess the cumulated activity 
between the second and third and then the third and fourth sacrifice intervals based on an 
integration between these points, and again assuming an exponential fit. The last term 
represents the cumulated activity from the final sacrifice interval to infinity based on 
physical decay as the removal mechanism from the organ. 

The above calculation requires determination of Ao and Xeff derived from the first 

two sacrifice activities and times, and are equal to: 

/ 

=exp 


^7; xln(^ 2 )-r 2 xln(^,)^ 


'^eff ^ 2 ) ~ 


T,-T, 

\n{A,)-\n{A,) 


T -T 
■'2 ■'1 


For determination of the cumulated activity using a two or single sacrifice interval 


method, the expressions for cumulated activity reduce to: 

A 


A(O-^oo) = 




-Jl - exp(->^, 




)] +1.44xr^ x4 


A (0^ 00 ) = 1.44 xTpx'^A^ exp(^^2j)j 


(2 points) 
(1 point) 


The source organs considered in this evaluation are listed in table 2.2. They were 
collected and measured as part of established protocols for determination of the 
biodistribution of a new radiotracer. Once the residence times in all source organs were 
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established, they were corrected using factors recommended by Roedler*'^ to accoimt for 
the difference in organ to total body weight proportions between standard rat and standard 
man. Appropriate S-values (rad/uCi-hr) for all important source-target pairs were 
obtained from the MIRDOSE III program*^, with the S-values for the rest of body 
assessed using the method of Coffey and Watson*^. The dose to a given target organ was 
then determined from the conventional MIRD formalism . Absorbed dose for all 
target organs listed in table 2.2 were determined in this manner. 

Table 2.2. Important Source and Target Organs Considered for PET 

Radiotracers 

Source Organs (organs harvested as part of sacrifice protocols) 

Adrenals Kidneys Ovaries Spleen Rest of Body 

Brain Liver Pancreas Testes 

Heart Lung Small Intestine Urinary Bladder Contents 

Target Organs (organs for which absorbed dose was calculated) 

Adrenals* Kidneys Pancreas Testes Red Marrow^ 

Brain Liver Spleen Ovaries Bone SurfacesJ 

Heart Wall Lungs Uterusf Total Body 

Urinary Bladder Wall Gallbladder Wall 

*Target organs highlighted were critical organs for one or more of the 
radiopharmaceuticals considered. 

t Uterine dose is assessed to approximate fetal dose, and assumes no transplacental 
transport of the radiotracer or its metabolites. 

X Dose to bone surfaces and bone marrow are estimated using remainder of body 
residence time. 

Urinary Bladder and Gallbladder 

Dosimetry for the urinary bladder and gallbladder required special consideration. 
Urinary bladder residence time was assessed in ligated male rats from the measured 
bladder contents. Female data was obtained from the radioactivity assessed as “missing” 
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in rat at the time of sacrifice. Missing radioactivity was that not accounted for in the 
measured organs and carcass, and was assumed to have been excreted by the animal 
during the waiting period between injection and sacrifice. The two measures were 
combined to calculate an effective cumulated activity, resulting in a conservative estimate 
of bladder residence time and subsequent absorbed dose per unit injected dose. For the 
purposes of this paper, a conservative estimate of a parameter indicates a value larger 
than the true value. 

The determination of a dose to the small intestine wall or gallbladder wall per unit 
injected radioactivity is impossible to assess in the rat since the animal has no gallbladder 
and excretion from the liver is shunted directly to the small intestine. The ICRP indicates 
that for most very short-lived radiopharmaceuticals, gastrointestinal residence times can 
be ignored for dosimetric purposes'*, since clearance from the human gallbladder is 
considered to occur in bolus amounts, with the first clearance time at three hours post 
injection which is long compared to the half-life of carbon-11. It is assumed the 
radioactivity directly shunted to the small intestine results in a small overall dose to the 
intestine because of delayed uptake and its large mass, and subsequently small S-value. 

A conservative determination of gallbladder wall dose was made using the ICRP-53 
biliary excretion model which assumes that thirty percent of the biliary excretion would 
be deposited in the gallbladder'*. Thus, thirty percent of the measured small intestine 
residence time was used as an estimate for the gallbladder residence time. 

Results 

The study included four to five animals per time point, with 13 animal organs 
evaluated at each time point for each of the ten radiotracers. The most prevalent 
biodistribution observed (65% of cases) was characterized by a rapid maximum uptake at 
Tl, followed by rapid clearance to the second sacrifice interval. After this time, removal 
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was equal to or somewhat greater than that expected by physical decay. This 
biodistribution is demonstrated in figure 2.4, and was most characteristic of the 
hydrophilic compounds such as [“C]epinephrine and [“C]mHED. This behavior was 
also seen for [**C]NMPB (lipophilic, amine), [**C]FNZPAM (lipophilic, amide) and 
[**C]raclopride, with the biodistribution exhibiting significant biological removal after 
the second sacrifice interval. Delayed organ uptakes were demonstrated for the liver and 
testes of all hydrophilic compounds and for most organ uptakes of [**C]TBZ, [**C]MTBZ 
and [**C]TRB, accounting for 20% of the cases examined. The radioactivity measured in 
the carcass and remaining tissues typically followed a mono-exponential profile of pure 
physical decay, accounting for 15% of the cases. 

Because biodistributions varied dramatically during the first 15 minutes post 
injection, the only acceptable reduced sacrifice permutation required use of the two 
sacrifice intervals measured within the first physical half-life of * 'C. This method had the 
advantage of preserving information during the period when biological uptake and 
clearance played a significant role determining the time-radioactivity profile. After the 
second time of sacrifice, physical decay was assumed to be the only removal mechanism. 
All other attempted permutations failed to adequately model the measured biodistribution 
resulting in considerable error in the estimated cumulated activity. 
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Organ Activity versus Time Post Injection 
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Figure 2.4. Biodistribution of [“CJEpinephrme in Rat Kidneys. Top Line 
Indicates Removal By Physical Decay Only After the Second Sacrifice Interval. 

The residence times determined using all four sacrifice intervals and only the first 
two sacrifice intervals are compared in table 2.3 for the three solid organs receiving the 
highest absorbed dose The residence time assessed from the simplified method varied 
from being 50% larger (conservative) to 15% smaller than that assessed using the all 
available time points. In the vast majority of cases, a conservative measure of residence 
time resulted as compared to using all data. For those cases where a positive percent 
difference was measured (i.e. the cumulated activity assessed from all data was greater 
than that assessed using the reduced method), the difference was small and typically there 
was no statistical difference between the values as indicated by a Student’s t-test. In less 
than 10% of the cases was the measured residence time using the new method lower than 
the assed using the conventional method, up to a maximum difference of 15%. 
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Table 2.3. Rat Residence Times for the Three Solid Organs Receiving the Largest 

Absorbed Dose 


Drug 

Source 

Organ 

All Data Assessed 

Th(S) ±<TTh(s)* 

Limited DataTI &T2 

Th(S) ±C7Xh(S)* 

%Difft 

EPl 

Heart 

32 ±2.7 

37 ±5.0 

-16% 


Kidney 

88 ±40 

95 ±36 

-7.8% 


Liver 

420 ±39 

440 ±76 

-4.8% 

mHED 

Heart 

39 ±3.3 

41 ±4.7 

-4.8% 


Liver 

370 ±37 

360 ±49 

0.35% 


Kidney 

63 ±48 

69 ±41 

-11% 

PHE 

Ovaries 

1.1 ±0.19 

1.1 ±0.18 

-3.2% 


Kidney 

41 ±5.0 

44 ±5.8 

-7.8% 


Heart 

10 ±0.74 

12 ±1.7 

-17% 

MTBZ 

Ovaries 

3.1 ±0.53 

2.9 ±0.90 

6.2% 


Adrenals 

3.5 ±0.63 

3.4 ±1.1 

1.0% 


Pancreas 

32 ±3.2 

37 ±6.6 

-18% 

NMPB 

Ovaries 

2.5 ±0.39 

3.3 ±0.47 

-29% 


Kidney 

68 ±13 

72 ±15 

-6.6% 


Brain 

74 ±6.2 

77 ±7.4 

-4.5% 

TBZ 

Adrenals 

2.9 ±0.48 

3.3 ±0.68 

-16% 


Ovaries 

2.1 ±0.26 

2.1 ±0.15 

-1.2% 


Pancreas 

65 ±5,8 

62 ±7.5 

-16% 

TRB 

Lung 

63 ±6.2 

75 ±11 

-17% 


Brain 

80 ±6.1 

93 ±12 

-16% 


Heart 

18 ±1.0 

23 ±1.3 

-24% 

FNZPAM 

Adrenals 

1.5 ±0.54 

2.4 ±0.91 

-55%$ 


Liver 

160 ±27 

210 ±45 

-28% 


Ovaries 

0.34 ±0.09 

0.39 ±0.08 

-15% 

PK11195 

Adrenals 

10 ±2.3 

9.5 ±4.1 

5.0% 


Heart 

72 ±3.4 

86 ±3.7 

-18% 


Ovaries 

5.8 ±1.8 

6.7 ±1.7 

-15% 

RAC 

Ovaries 

2.1 ±0.47 

2.3 ±0.65 

-8.9% 


Kidney 

55 ±8.2 

67 ±9.6 

-20% 


Liver 

230 ±15 

270 ±21 

-20% 


* The standard deviation of the measured residence time in seconds. 

fThe percent difference between the measured residence resulting from the two methods, 
where a negative percent difference indicates a conservative (higher) result using the 
simplified technique 

X Anomalous difference based on small organ activities and difficulty in precise 
determination of adrenal organ activity. 
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Table 2.4. Estimated Human Absorbed Dose per Unit Administered Activity for the 
Three Solid Target Organs Receiving the Largest Dose 



Target 

All Data Assessed 


Limited Data T1 & T2 



Drug 

Organ 

mGy/MBq ±ct* 

rad/mCi ±cy* 

mGy/MBq ±a* rad/mCi ±a * 

%Difft 

EPI 

Heart Wall 

0.016 ±0.0012 

(0.060 ±0.0043) 

0.019 ±0.0024 

(0.069 ±0.0087) 

-15% 


Kidneys 

0.012 ±0.0043 

(0.045 ±0.016) 

0.013 ±0.0045 

(0.049 ±0.017) 

-7.9% 


Liver 

0.010 ±0.0009 

(0.039 ±0.0034) 

0.011 ±0.0017 

(0.041 ±0.0065) 

-5.2% 

mHED 

Heart Wall 

0.019 ±0.0018 

(0.073 ±0.0066) 

0.021 ±0.0022 

(0.076 ±0.0082) 

-4.9% 


Liver 

0.0092 ±0.0007 

(0.034 ±0.0025) 

0.0092 ±0.0011 

(0.034 ±0.0042) 

-0.19% 


Kidneys 

0.0091 ±0.0047 

(0.034 ±0.017) 

0.010 ±0.0051 

(0.037 ±0.019) 

-9.8% 

PHE 

Ovaries 

0.0073 ±0.0007 

(0.027 ±0.0027) 

0.0071 ±0.0008 

(0.026 ±0.0028) 

2.8% 


Kidneys 

0.0063 ±0.0006 

(0.023 ±0.0022) 

0.0068 ±0.0007 

(0.025 ±0.0027) 

-7.5% 


Heart Wall 

0.0060 ±0.0005 

(0.022 ±0.0018) 

0.0070 ±0.0008 

(0.026 ±0.0031) 

-16% 

MTBZ 

Ovaries 

0.015 ±0.0023 

(0.055 ±0.0084) 

0.014 ±0.0038 

(0.052 ±0.0141) 

4.8% 


Adrenals 

0.0145 ±0.0025 

(0.054 ±0.0092) 

0.015 ±0.0039 

(0.054 ±0.0144) 

-0.25% 


Pancreas 

0.0089 ±0.0007 

(0.033 ±0.0027) 

0.010 ±0.0015 

(0.038 ±0.0054) 

-16% 

NMPB 

Ovaries 

0.013 ±0.0017 

(0.048 ±0.0063) 

0.016 ±0.0020 

(0.059 ±0.0074) 

-22% 


Kidneys 

0.0096 ±0.0017 

(0.036 ±0.0063) 

0.010 ±0.0019 

(0.038 ±0.0071) 

-6.6% 


Brain 

0.0095 ±0.0008 

(0.035 ±0.0029) 

0.010 ±0.0009 

(0.037 ±0.0034) 

-4.9% 

TBZ 

Adrenals 

0.012 ±0.0018 

(0.045 ±0.0068) 

0.014 ±0.0025 

(0.051 ±0.0092) 

-14% 


Ovaries 

0.011 ±0.0010 

(0.039 ±0.0036) 

0.011 ±0.0007 

(0.039 ±0.0024) 

-0.06% 


Pancreas 

0.0084 ±0.0009 

(0.031 ±0.0033) 

0.0095 ±0.0015 

(0.035 ±0.0057) 

-13% 

TRB 

Lung 

0.012 ±0.0011 

(0.043 ±0.0042) 

0.014 ±0.0018 

(0.050 ±0.0065) 

-16% 


Brain 

0.010 ±0.0007 

(0.039 ±0.0027) 

0.012 ±0.0016 

(0.045 ±0.0058) 

-15% 


Heart Wall 

0.010 ±0.0006 

(0.038 ±0.0022) 

0.013 ±0.0006 

(0.046 ±0.0024) 

-21% 

FNZPAM 

Adrenals 

0.0070 ±0.0019 

(0.026 ±0.0071) 

0.010 ±0.0033 

(0.037 ±0.012) 

-44% 


Liver 

0.0045 ±0.0005 

(0.017 ±0.0020) 

0.0056 ±0.0010 

(0.021 ±0.0038) 

-23% 


Ovaries 

0.0041 ±0.0004 

(0.015 ±0.0016) 

0.0044 ±0.0004 

(0.016 ±0.0015) 

-7.2% 

PK11195 

Adrenals 

0.039 ±0.0079 

(0.14 ±0.029) 

0.037 ±0.015 

(0.14 ±0.054) 

4.2% 


Heart Wall 

0.036 ±0.0015 

(0.14 ±0.0055) 

0.043 ±0.0018 

(0.16 ±0.0065) 

-17% 


Ovaries 

0.027 ±0.0067 

(0.10 ±0.025) 

0.032 ±0.0073 

(0.12 ±0.027) 

-13% 

RAC 

Ovaries 

0.011 ±0.0021 

(0.040 ±0.0076) 

0.012 ±0.0028 

(0.043 ±0.010) 

-7.5% 


Kidneys 

0.0083 ±0.0010 

(0.031 ±0.0036) 

0.0097 ±0.0012 

(0.036 ±0.0045) 

-17% 


Liver 

0.0063 ±0.0004 

(0.023 ±0.0013) 

0.0073 ±0.0005 

(0.027 ±0.0018) 

-15% 


*The standard deviation of the measured organ dose per unit injected radioactivity in 
mGy/MBq (mrad/mCi). 

fThe percent difference between the calculated absorbed dose per unit injected dose 
resulting from the two sacrifice methods. 






25 


The comparison of absorbed dose per unit injected radioactivity for the three solid 
organs receiving the highest dose from each radiotracer is shown in table 2.4. Of the 13 
organs for each of the ten radiotrcaers evaluted, the average percent difference between 
the reduced sacrifice technique and the extensive technique was 8%. The range of 
differences between the two methods was from 55% larger (conservative) to 19% smaller. 
In over half the cases, there was less than a 10% difference between the two methods, and 
in over 85% of the cases there was less than a 20% difference between the two methods. 
Overall the reduced sacrifice method lead to a conservative measure of absorbed dose per 
unit injected activity in 85% of the cases, excluding the gallbladder wall dose and the 
urinary bladder wall dose. In only 15% of the cases was a positive percent difference 
measured, and these differences were within the standard deviation of the calculated 
organ doses. 

Urinary Bladder and Gallbladder 

Residence times for urinary bladder contents and estimated gallbladder contents, 
and dosimetry for the gallbladder wall and urinary bladder wall are presented in tables 2.5 
and 2.6, respectively. The standard deviation in urinary bladder residence times were on 
the order of 5% to 50%, and were reasonable considering the expected variability 
between male bladder contents and the radioactivity assessed as missing from the female 
rats. The standard deviation for the “gallbladder contents” ranged from 10% to 20%. 

The bladder data indicated a consistent negative bias between the multiple time-point and 
reduced sacrifice methods for many of the radiotracers, on the order of thirty to forty 
percent (i.e. the reduced sacrifice method provided a smaller measure of residence time 
than the current method). This is largely attributable to the delayed filling of these organs 
which is not assessed when using a reduced number of sampling points. Similar results 
are observed for the estimated gallbladder residence time and gallbladder wall dose. A 
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maximum negative bias of 20% was measured between the four time-point and two time- 
point methods, again indicative of the delayed uptake in the “gallbladder” from the liver- 
biliary pathway. 


Table 2.5. Rat Residence Times for the Urinary Bladder and Gallbladder 


Drug 

Source 

Organ 

All Data Assessed 

Th (s) ± a Th (s) * 

Limited Data T1 & T2 

Th (s) ± a Th (s) * 

Bladder 

Content 

%DifftMethod* 

EPl 

U. Bladder 

240 ±79 

250 ±20 

-5.7% 

1 


Sm. Intest. 

122 ±24 

110 ±23 

9.7% 


mHED 

U. Bladder 

260 ±23 

160 ±51 

39% 

2 


Sm. Intest. 

n.a. 

n.a. 



PHE 

U. Bladder 

450 ±27 

300 ±26 

33% 

3 


Sm. Intest. 

82 ±18 

73 ±7.5 

11% 


MTBZ 

U. Bladder 

18 ±0.96 

17 ±0.37 

3.7% 

2 


Sm. Intest. 

95 ±12 

96 ±17 

-1.3% 


NMPB 

U. Bladder 

390 ±37 

270 ±73 

32% 

4 


Sm. Intest. 

320 ±29 

267 ±46 

15% 


TBZ 

U. Bladder 

110 ±25 

71 ±32 

36% 

4 


Sm. Intest. 

66 ±5.8 

62 ±7.2 

5.1% 


TRB 

U. Bladder 

138 ±30 

87 ±60 

37% 

4 


Sm. Intest. 

n.a. 

n.a. 



FNZPAM 

U. Bladder 

360 ±58 

400 ±71 

-9.7% 

3 


Sm. Intest. 

n.a. 

n.a. 



PK11195 

U. Bladder 

330 ±100 

350 ±178 

-6.9% 

1 


Sm. Intest. 

n.a. 

n.a. 



RAC 

U. Bladder 

n.a. 

n.a. 




Sm. Intest. 

160 ±29 

130 ±50 

20% 



* The standard deviation of the measured residence time in seconds. 

t The percent difference between the measured residence resulting from the two methods. 

X Method by which animal urinary bladder radioactivity was determined: 1) Males used 
measured bladder contents, females used unrecovered radioactivity, 2) Bladder 
radioactivity based on male bladder contents, essentially 100% recovery in females, 3) 
Urine not collected, bladder radioactivity estimated from missing radioactivity at time of 
sacrifice, females only, 4) Urine not collected, bladder radioactivity estimated from 
missing radioactivity at time of sacrifice, males and females 

n.a. = not available 
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Table 2.6. Absorbed Dose per Unit Administered Activity for the Gallbladder Wall 

and Urinary Bladder Wall 


Drug 

Target 

Organ 

All Data Assessed 

mGy/MBq ±a* 

rad/mCi ±c* 

Limited Data T1 & T2 

mGy/MBq ±ct* rad/mCi ±a * 

%Difft 

EPI 

UB Wall 

0.046 ±0.015 

(0.17 ±0.055) 

0.049 ±0.0037 

(0.18 ±0.014) 

-6.1% 


GB Wall 

0.026 ±0.0046 

(0.098 ±0.017) 

0.024 ±0.0045 

(0.090 ±0.017) 

7.6% 

mHED 

UBWall 

0.050 ±0.0044 

(0.19 ±0.016) 

0.031 ±0.0096 

(0.12 ±0.035) 

38% 


GB Wall 

n.a. 

n.a. 

n.a. 

n.a. 


PHE 

UB Wall 

0.085 ±0.0050 

(0.32 ±0.018) 

0.058 ±0.0050 

(0.22 ±0.018) 

32% 


GB Wall 

0.018 ±0.0035 

(0.068 ±0.013) 

0.017 ±0.0015 

(0.062 ±0.0054 

8.6% 

MTBZ 

UB Wall 

0.006 ±0.0002 

(0.021 ±0.0008) 

0.006 ±0.0001 

(0.021 ±0.0005 

\ 

-0.82% 


GB Wall 

0.021 ±0.0023 

(0.079 ±0.0085) 

0.022 ±0.0033 

/ 

(0.081 ±0.012) 

-2.3% 

NMPB 

UBWall 

0.074 ±0.0069 

(0.28 ±0.026) 

0.051 ±0.014 

(0.19 ±0.051) 

31% 


GB Wall 

0.063 ±0.0056 

(0.23 ±0.021) 

0.054 ±0.0089 

(0.20 ±0.033) 

14% 

TBZ 

UB Wall 

0.023 ±0.0047 

(0.085 ±0.018) 

0.016 ±0.0059 

(0.058 ±0.022) 

32% 


GB Wall 

0.016 ±0.0011 

(0.058 ±0.0042) 

0.015 ±0.0015 

(0.056 ±0.0054 

\ 

3.7% 

TRB 

UB Wall 

0.028 ±0.0056 

(0.10 ±0.021) 

0.019 ±0.011 

/ 

(0.069 ±0.042) 

34% 


GB Wall 

n.a. 

n.a. 

n.a. 

n.a. 


FNZPAM 

UB Wall 

0.070 ±0.011 

(0.26 ±0.040) 

0.076 ±0.013 

(0.28 ±0.049) 

-9.3% 


GB Wall 

n.a. 

n.a. 

n.a. 

n.a. 


PK11195 

UB Wall 

0.064 ±0.019 

(0.24 ±0.069) 

0.068 ±0.033 

(0.25 ±0.12) 

-6.1% 


GB Wall 

n.a. 

n.a. 

n.a. 

n.a. 


RAC 

UB Wall 

n.a. 

n.a. 

n.a. 

n.a. 



GB Wall 

0.035 ±0.0055 

(0.13 ±0.021) 

0.028 ±0.0096 

(0.10 ±0.036) 

19% 


*The standard deviation of the measured organ dose per unit injected radioactivity in 
mGy/MBq (mrad/mCi). 

fThe percent difference between the calculated absorbed dose per unit injected dose 
resulting from the two sacrifice methods. 

n.a. = not available 


Discussion 

The biodistributions observed for ten carbon-11 labeled radiopharmaceuticals 
suggested that a time-dependent model could be used for solid organs that adequately 
described the variable initial uptake and clearance of the radiotracer followed by mono- 
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exponential clearance. Application of a reduced sacrifice technique based on this model 
made use of two sacrifice times at 2 or 5 minutes and 10 or 15 minutes post injection. 
Assessment of residence time and absorbed dose per unit injected dose from this 
technique provided results in good agreement with using extensive sacrifice methods. 
This proposed simplification in biodistribution assessment provides several advantages 
over that proposed by Gatley*^, which determined a theoretical upper limits of dose based 
on cardiac output models. Use of only modeled biodistribution can be overly 
conservative and does not necessarily provide good agreement with the radiotracer’s 
measured biodistribution. With the use of animal data, conservatism is maintained, and 
actual organ activities are relied upon for dose estimation. 

The magnitude of a given organ dose was most dependent on the measured 
residence time for that organ and its self-irradiation S-value. This is a consequence of the 
dose from ”C being principally due to energy deposition by the positron, with relatively 
little energy deposited by annihilation radiation^®. S-values for organ self-irradiation are 
subsequently three to five orders of magnitude higher then when a target organ is 
irradiated by a different source organ. This observation allows some flexibility in the 
protocols used to establish organ residence times for positron emitters since external 
source organs will contribute negligible dose to a given target. Organs which are 
dosimetrically unimportant can be assessed quickly using a reduced technique, and 
protocols tailored to focus on critical organs such as the gonads or bladder organs. 

In general the liver, “remainder of body” and urinary bladder contents had the 
largest residence times regardless of compound and were typically one or two orders of 
magnitude higher than other organs assessed. However, because of the inverse 
relationship between organ mass and S-value, smaller body organs were often found to 
have the highest absorbed dose per unit injected radioactivity. Solid organs which 
consistently demonstrated the highest dose per unit injected radioactivity included the 
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liver, kidneys and heart for hydrophilic compounds, and the ovaries, adrenals, kidneys 
and lung for lipophilic compounds, dependent on the designed receptor site. 

Two organs which demonstrated consistently large cumulated activities were the 
urinary bladder contents and “estimated” gallbladder contents, with high absorbed doses 
resulting for the urinary bladder wall and gallbladder wall, respectively. Comparing the 
absorbed dose data between the bladder organs and solid organs, it is apparent that the 
urinary bladder is frequently the critical organ limiting administered dose. The 
exceptions include ["CJraclopride for which the gallbladder is limiting and ['‘CJMTBZ 
where the ovaries are limiting. However, considerable difficulty was encountered in 
determining a time-radioactivity profile for these organs. For the case of the urinary 
bladder, only two or three male rats are currently used per time-point. Because the 
measured urine radioactivity can vary significantly between these animals due to 
physiological differences, a valid residence time was difficult to establish. The current 
method does not address any diluting volume of urine at the time of ligation, and the 
bladder may be so expanded from urine accumulation at the time of sacrifice that it may 
affect measured kidney radioactivity. Further, the ligation may fail or the animal can 
remove the ligation if not properly restrained causing a toss of data. Instead, urinary 
bladder residence time was conservatively assessed from a combination of the available 
male data and the radioactivity assessed as “missing” in female rats at the time of 
sacrifice. The conservatism introduced by using only missing data from the female is 
observed for [^^C]phenylephrine and [**C]FNZPAM where no male data were acquired. 
The assessed bladder residence times for these two cases were considerably larger than 
any other radiotracer considered. Ideally, the use of customized metabolic cages for 
collection of excreted material is desired, but these can greatly complicate the 
measurement process. Excluding the use of metabolic cages, the current method provides 
results which show surprisingly small variance and lead to conservative estimates of 
human radiation dose. Similar observations are made in assessing the residence time of 
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the gallbladder contents where radioactivity was observed to vary markedly between 
radiopharmaceutical category (lipophilic verses hydrophilic) as well as between animals. 
The dose calculated for the gallbladder wall was found to be significant for most 
radiopharmaceuticals examined, second only to the dose to the urinary bladder wall. The 
additional complication of rat bladder dosimetry is its applicability to the human. The 
well known metabolic differences between rat and human contribute additional 
conservatism to the measured bladder wall doses . 

Because of the variability in bladder organ time-radioactivity profiles, the 
described reduced sacrifice technique requires particular consideration when applied to 
bladder organs where uptake is clearly delayed and organ concentration is not readily 
determined by a simple model. Observation of tables 2.5 and 2.6 indicates that correction 
of the measured urinary bladder residence time by 60% will result in good agreement 
between the two methods, with a range in bias from 0 to 50% (conservative). A similar 
correction of 25% can be applied to gallbladder residence times assessed with the two- 
point method to correct its negative bias, providing a range of differences from 0 to 22% 
(conservative). If necessary, refinement of the estimated urinary bladder dose can be 
performed using metabolic cages and application of the dynamic bladder model'^. 
Unfortunately, the short half-life of prevents a similar method for invitro verification 
of the gallbladder wall dose. Recent works have discussed the application of high 
resolution PET imaging for small animal work^^ Invivo quantification of urinary 

bladder and gallbladder contents may someday prove to be a better method of measuring 
organ radioactivity non-invasively with improved time resolution and accuracy compared 
to sacrifice methods. Alternatively, bladder data are perhaps best determined “after the 
fact” from actual human subjects based on invivo quantification and invitro measurement 
of samples. Sinee the urinary bladder was commonly the limiting organ in terms of 
allowable injectable dose, collection of actual human data is important to validate the 
conservative and error prone animal data. 




31 


Recent recommendations of the ICRP^'^ have supported the use of effective dose 
in the determination of risk verses benefit in biomedical research. The current tissue 
weighting factor for the gonads is four times that of the urinary bladder^^, placing greater 
significance in the gonadal dose in determining limiting patient dosage. Should present 
regulatory bodies adopt ICRP recommendations, increased emphasis will be placed in the 
accurate determination of both gonadal dose and bladder doses, since these organs were 
found to receive the maximum effective dose overall. Since both organs are assessed 
with limited precision in mixed sex animal studies, validation of these organ 
biodistributions in higher animals or during initial human trials will have increased 
importance. 


Conclusions 

Measurement of residence time can be performed conservatively and quickly 
using a reduced sacrifice method focusing on two intervals within one half-life post 
injection. The method essentially reduces the time and animal sacrifice requirements by 
half, while still maintaining necessary accuracy. The applicability of a reduced sacrifice 
method can be extended to other short or very short half-life radiopharmaceuticals, such 
as those labeled with 10 minute half-life N. Considering the inaccuracies in using rat 
data to characterize human biodistributions, the use of this reduced sacrifice method 
provides results which are quite acceptable. 

The measured biodistributions of both hydrophilic and lipophilic compounds 
indicates that the urinary bladder wall and gallbladder wall can be dosimetrically 
important. A correction factor of 1.6 is required for urinary bladder residence times and 
1.25 for estimated gallbladder residence time assessed using the reduced sacrifice 
technique due to the delayed filling of these organs. Using the regulatory limits of 3 rads 
to the gonads and lens of eye and 5 rads to other organs , the smallest injectable 
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radiopharmaceutical dosage was found to be 410 MBq (11 mCi) for ['*C]FNZPAM, 
where the bladder wall dose was limiting. All other radiotracers have considerably higher 
permissible activities. These doses should be sufficient for initial evaluation of a new 
radiotracer on current two-dimensional imaging systems. The advent of septaless three- 
dimensional imaging PET systems will provide significant improvements in sensitivity, 
by as much as a factor of six to ten . This will subsequently allow researchers to 

investigate the performance of new radiotracers with doses well below the conservative 
dose limits assessed from the reduced sacrifice method. Investigators should verify 
bladder time-radioactivity profiles in human subjects during initial imaging trials to refine 
the recognized conservative bladder measurements made in the rat. 
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CHAPTER III 


NON-INVASIVE ASSESSMENT OF ANIMAL BIODISTRIBUTIONS 

Invivo assessment of the organ activity of a radiotracer is desirable because it can 
trace biodistribution within the same subject, relieving problems associated with inter¬ 
animal variability. The method could specifically facilitate the assessment of organs with 
complex biodistribution patterns such as the gallbladder and urinary bladder, as well as be 
useful for the dosimetry of tumors where a sequential history of tumor uptake and 
removal in the same animal is required. Historically, several methods have been used for 
invivo assessment of organ burdens including use of a single external detector (e.g. 
thyroid burden and bladder activity assessment), use of rectilinear scans, and Anger 
camera scintigraphy. Conjugate-view counting remains a common tool for invivo 
quantification. These conventional methods however suffer from interference by 
backgroimd radiation sources extraneous to the organ of interest, i.e. from other tissues 
overlying, underlying and adjacent to the organ being counted’. 

The development of emission computed tomography (ECT) 30 years ago^ has 
provided a methodology for obtaining estimates of three dimensional activity 
distributions in a given tissue. Two methods of ECT are currently in use, single photon 
emission computed tomography (SPECT) and positron emission tomography (PET). In 
the prior case, nuclides emitting single gamma rays are measured via planar scintigraphy 
at multiple angles around the patient. A general SPECT system is shown in figure 3.1. 

In SPECT imaging, a patient is administered a gamma ray emitting radiopharmaceutical. 
The drug is designed to preferentially localize in a given tissue or to follow a certain 
physiologic process. Photons emitted by the radionuclide in the patient are detected by a 
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sodivim iodide (Nal) gamma camera, where lead collimation and detector position 
electronics are used to localize the ray along which the gamma ray was emitted. The 
camera thus forms a scintigraph, or projection image of the pharmaceutical distribution in 
the patient. Collimation can be interchanged on the gamma camera to achieve optimal 
sensitivity and resolution for the energy of the gamma emitter and the target tissue under 
consideration. The camera is rotated about the patient in a series of steps, and a new 
planar image is collected at each location. The collected images are reconstructed 
through a process of Fourier analysis, convolution techniques, or iterative reconstruction 
techniques to arrive at a transaxial image (or slice view) of the patient. A three- 
dimensional image of the radioactive soxirce distribution can be obtained by stacking and 
further manipulating the reconstructed slices. In PET, a positron emitting 
radiopharmaceutical is injected into the patient, and annihilation radiation is detected by a 
ring of inorganic scintillators (typically bismuth germinate, or BGO) which are connected 
to a complex coincidence circuit. Only annihilation photons arriving at transposed 
detectors within a specific time interval (< 10 nsec) are recorded and used for producing 
the reconstructed transaxial image. A critical difference between SPECT and PET is the 
use of conventional lead collimation to restrict the field of view in SPECT systems. In 
PET, collimation is primarily accomplished electronically through the coincidence 
counting circuit^. A significant advantage of PET over SPECT is that electronic 
collimation is at least an order of magnitude more sensitive than mechanical collimation, 
£ind sensitivity is not a function of spatial resolution. 

If scatter and attenuation in the patient are temporarily ignored, the gamma ray 
emissions counted via SPECT or PET are proportional to the line integral of activity 
viewed by the camera. If resolution distance is small (i.e. high resolution) compared to 
the size of uniform objects being imaged, the reconstructed image will have a voxel count 
rate which is proportional to the concentration of activity in the organ or tissue. This 
invivo measurement method thus removes the effect of confounding overlying and 
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underlying activity common in planer imaging, and permits dynamic visualization of a 
given physiological process in the same patient. 
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Figure 3.1. A Generalized SPECT System Featuring a Parallel Hole Collimator 


Accurate quantification of radionuclide distribution by emission computed 
tomography is dependent on numerous factors, most notably system resolution and 
sensitivity. The system extrinsic resolution determines the minimum size object which 
can be quantified accurately and good system resolution is a requisite to maximize image 
contrast and organ definition. Poor system resolution causes bias in uptake estimates and 
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degrades image quality. In general, system resolution depends on 1) geometric and 
electronic resolution of the detector collimator-system, 2) the adequacy of the sampling 
geometry, and 3) the range of the positrons (approximately 360 mg/cm^ for maximum 
energy positrons). This latter effect cannot be controlled and introduces a loss of at 
least one to two millimeters in resolution. The intrinsic detector resolution and the 
collimator resolution combine (approximately in quadrature) to define the limiting system 
resolution. However, there is an inverse relationship between sensitivity and resolution. 
Projection image data are inherently deficient in counts due to limits on injected activity, 
the specificity of the compound for the target tissue, the counting time and the limited 
system sensitivity. This can result in large statistical fluctuations in the imaged photons, 
which in turn results in noise in the reconstructed image. High sensitivity is desired to 
minimize image noise, but this is typically only accomplished by increasing the detector 
thickness for high energy photons or increasing collimator sensitivity each of which 
results in a commensurate reduction in resolution. Accurate quantification also requires 
attention to factors which impact measurement of true line integrals. These include 
effects of attenuation within the patient, errors due to accidental coincidence in PET, 
effects of scattered radiation in the object, environmental background and septal or 
aperture penetration of the measured gamma rays through the lead collimation'*. 

Current systems used clinically for PET and SPECT have neither sufficient 
sensitivity nor resolution to accurately quantify time dependent radionuclide distributions 
in small laboratory animals, such as the rat. Work by Rogers et al. at the University of 
Michigan has focused on development of a small animal SPECT system in order to 
accelerate and simplify the design and evaluation of new radiopharmaceuticals. This 
research will focus on characterizing the performance of this instrument for high energy 


imaging. 





Figure 3.2. Picture and Exploded View of SPRINT-II Gantry, Slice Collimator and 
Slit Aperture. The cylinder at the center of SPRINT represents the approximate 

field-of-view. 
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Characteristics of SPRINT-II: A Second Generation Single Photon Ring 

Tomograph 

The development of new single photon emitting radiotracers for brain imaging 
such as 123j_i]y[p^ and 99xc-HMPAO, and for new tumor-specific therapy 

agents such as and various iodine labeled monoclonal antibodies have 

generated the requirement for a dedicated single photon imaging system with better 
sensitivity and resolution than is available with general purpose imaging systems. Work 
by Rogers et. al. at the University of Michigan has focused on development of a 
stationary detector single photon ring tomograph (SPRINT-II) designed specifically for 
human brain imaging^. Revision of the collimator design and other characteristics of the 
original system has expanded its use for the imaging of small laboratory animals. The 
full ring design of SPRINT provides two specific advantages in SPECT. First, multiple 
projection data can be simultaneously sampled improving sensitivity and temporal 
resolution, while reducing temporal blurring in dynamic studies. Second, the fixed ring 
geometry maintains a rigid center of rotation and reduces electronic instability of the 
detectors, providing higher quality projection data. 

Design and Construction 

An exploded figure of the SPRINT-II gantry as configured for high-energy animal 
imaging is shown in figure 3.2. For animals, the illustrated system has an 8 cm (3.51") 
diameter field of view. This diameter is defined by the dimension of a rotating multi-slit 
aperture ring, currently configured with 6 slits at 60° intervals as depicted in figure 3.3. 
Projection data is acquired by stepping the aperture ring over a 60° arc which provides a 
full 360° projection set. The slit edges are made of depleted uranium (3.175 mm (0.125") 
on a side), embedded in lead for improved resolution at high energies. Slit spacing is 
cam-controlled and can be adjusted radially to change resolution. In order to enhance 
resolution for small animals and simultaneously utilize the available detector area, the slit 
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aperture placement is designed to increase magnification of object projections onto the 
detectors in inverse proportion to the size of the animal. The small animal aperture has a 
6.5 cm slit radius and a 25.1 cm detector radius, providing an image-object magnification 
factor of three to one. The maximum field-of-view is defined by the projection of the 
object through adjacent slits, such that no overlap occurs of the rays in the detector plane. 
This is shown in figure 3.3. The salient characteristics of SPRINT are presented in table 
3 . 1 ^ 

Encircling the slit collimator is a slice collimator, a parallel or converging plate z- 
axis collimator which is fabricated from alternating layers of air equivalent closed cell 
foam and lead foil. Design, construction and performance of this collimator is covered in 
Chapter Five. The aperture and slice collimator work in tandem to determine the system 
extrinsic resolution. In-plane image resolution is determined by the slit aperture which 
provides collimation along the circumference (x-axis of the detector array) of the system. 
Axial image resolution will be limited by the slice collimator which provides collimation 
along the z-axis (central axis) of the system. 

Surrounding the slice collimator are eleven two dimensional sodium iodide 
(Nal(Tl)) camera modules which use maximum likelihood position estimation logic’. 
They are arranged in a 50 cm diameter ring with a scintillator packing fraction of 96 
percent. Each detector module is nominally composed of 44 Nal(Tl) bars 3 mm wide, 13 
mm deep and 15 em long. The bars lie parallel to the central axis of the gantry and the 
slit aperture. The 13 mm depth Nal(Tl) bars are comparable to the thickness of older 
conventional Anger cameras, designed optimally for the detection and resolution of 140 
keV (Tc-99m) photons. Photoelectric detection efficiency for 511 keV photons in 13 mm 
of Nal(Tl) is estimated to be between only 7 percent, compared to 94% for 140 keV 
photons. 
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Table 3.1. SPRINT-II Salient Characteristics 


No. Detector Modules 

11 

PMTs per Modules 

20 

Module Dimensions 

14 cm X 15.2 cm 

Detector Radius 

25 cm 

Slit Aperture Radius 

6.5 cm 

Magnification 

3-to-l 

Number Slits/Angular 

Separation 

6/60° 

Slit Width 

0-4 mm 

Slice Collimation 

Low-Energy Parallel (140 keV) 

Medium Energy Converging (364 keV) 

High Energy Parallel (511 keV) 

Axial FOV 

9 cm 

In-Plane FOV 

9 cm 


System resolution is dependent on the combined effects of inherent detector 

resolution, and the imaging geometry. The maximum resolution obtainable for SPRINT 

can be estimated with the help of Figure 3.4. 

Given a point source emitting low energy photons at the center of the field of 

view, it will project an image onto the detector given by: 

( 5 ^ + Sj^) 

K =-^- 

For 511 keV imaging, a slit width of 3.175 mm (0.125”) has been used'. This 
results in Ra = 12.3 mm. This geometric resolution must then be convolved with the 
intrinsic resolution provided by the detector module, which is approximately equal to the 


' The choice of this slit width and discussion of penetration effects will be presented in 
Chapter VI. 




45 


detector bar width of 3 mm. Assuming spatial invariance and a Gaussian response 
function, the resolution at the detector is given by: 



Figure 3.4. Calculation of Theoretical In-Plane System Resolution 


R = + Rj = \2.6mm 

When this is projected back through the slit aperture onto the object space, the 

system resolution is determined as r, equal to: 

SgR . . 

r = —— = AAmm 

So 

Thus, assuming the object has a sufficient signal to noise ratio in the reconstructed 
image, the minimum size object which can be resolved is approximately 4.4 mm. A 
similar analysis can be used to determine the axial resolution of SPRINT. Since axial 
resolution is dependent on the slice collimator, the achievable axial resolution will be 
dependent on the collimator used. From figure 3.5, the collimator resolution is calculated 
as: 

^ SJS.+Sj,) 


which when convolved with the detector resolution of 4 mm, results in a predicted axial 
resolution of 2.3 mm for the medium energy converging collimator, 5.4 mm for the high 
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energy collimator and 9.3 mm for the low-energy collimator. Slice collimator 
performance will be reviewed in more detail in Chapter 5. 



Based on a predicted system resolution of 5 mm for high energy imaging, it is 
informative to compare resolution with average organ dimensions in the typical 
laboratory rat to predict some of the limitations with SPRINT II assessment. Table 3.2 
provides the typical dimensions and weights of organs in an adult rat with an adult weight 
of 200 to 250 grams. From the table, it is observed that quantification of activity in the 
adrenals, ovaries and pancreas would be very difficult due to either their small size or 
complex shape. Larger organs should have good potential for quantitative assessment, 
but measurement of heterogeneity of uptake within larger organs is still limited by system 
resolution. 

Since rat organs are not isolated from activity in adjacent tissues, activity 
concentration in the central region of the organ is the best achievable measure available 
from animal SPRINT images. Total organ activity must then be assessed by measuring or 
estimating rat organ volume, and assuming that organ concentration is homogeneous. 
Several possible alternatives exist to make this determination including using average 
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organ weights or volumes from "standard" rat corrected to the weight of the animal used 
in the study, estimating organ volumes from high resolution MR images, and sacrificing 
the animal at the end of imaging and measuring the organ volume directly. 

For organs which are well isolated from other activity in the body, e.g. the heart 
and brain for radiotracers which show preferential uptake in these organs and little uptake 
in the adjacent tissues, it is possible to use large regions of interest around the entire 
organ to estimate total organ activity. However, such estimates are sensitive to system 

resolution limits and partial volume effects. 

Table 3.2. Anatomical Data for the Laboratory Rat* 


Organ 

Size (mm) 

Weight (gm) 

Comments 

Adrenals 

4-5.5 mm long. 

0.02 -0.022 gm 

Right covered by 


3-4.5 mm wide. 


caudate lobe of liver. 


2.5-3 mm thick 


Both near kidneys. 

Brain 

8-10 mm high. 

1.5 - 2.2 gm 


(cerebrum) 

12-14 mm wide, 

15-20 nun wide 



Eye 

4 mm diameter 

0.12 gm (ea) 


Heart 

20 mm long, 

10-15 mm wide, 

10 mm deep 

0.9 -1.0 gm 


Kidney 

15-25 mm long, 

10-15 mm wide, 

10 mm deep 

0.73 -1.2 gm 


Liver 

20-40 mm wide. 

8.3 gm for 200 

4 lobes, uneven 


2-10 cm deep 

gm rat 

structure 

Lung 

20 mm high 

0.93-2.13 gm 

uneven structure. 


10 cm diameter 


surrounds heart 

Ovaries 

3-5 mm diameter 

0.06 gm 


Pancreas 


0.5 -1 gm 

very uneven structure, 
follows liver contour 

Spleen 

30-50 mm long, 

10 mm diameter 

0.75-1.3 gm 


Testes 

Up to 20 mm long, 

Up to 14 mm diam. 

2.0 - 3.5 gm 


Urinary 

Empty; Pear Shaped, 



Bladder 

10 mm long, 5 mm wide 
Full: Round, 

15-25 mm diameter 









48 


References to Chapter III 


* Budinger TF, Derenzo SE, Huesman RH. Role of tomography in providing radionuclide 
distribution and kinetic data. Third International Radiopharmaceutical Dosimetry 
Symposium, Oak Ridge National Laboratory, HHS-Publ (FDA) 81-8166; 60-74,1980. 

^ Kuhl DE, Edwards ED. Image separation isotope scanning. Radiology 1963;80:653-62. 

^ Sorenson JA and Phelps ME, Physics in Nuclear Medicine, 2nd Ed. Philadelphia: W.B. 
Saunders Company; 1987;391-423. 

Jaszczak RJ, Greer KL, Coleman RE. SPECT quantification of regional radionuclide 
distributions. Fourth International Radiopharmaceutical Dosimetry Symposium, Oak 
Ridge Associated Universities, CONF-851113-(DE86010102), 1985;82-96. 

^ Rogers WL, Clinthome NH, Shao L, Chiao P, Ding Y, Stamos JA, Koral KF. SPRINT 
II: A second generation single photon ring tomograph. IEEE Trans Med Imag. 
1988;7:4:291-297. 

^ Rogers WL. ASP: A small animal tomograph for radiotracer evaluation. University of 
Michigan, Proposal research summary. 1992. 

’ Clinthome NH, Rogers WL, Shao L, Koral KF. A hybrid maximum liklihood computer 
for scintillation cameras. IEEE Tran. Nuc Sci 1987;NS-34:97-101. 

* Hebei R, Stromberg MW, Anatomy of the Laboratory Rat, New York; Williams and 
Wilkins: 1976. 



CHAPTER IV 


TWO DIMENSIONAL IMAGE RECONSTRUCTION FOR SPRINT 


Data acquisition in SPRBSfT is performed by incrementally rotating the slit 
aperture n times through a total arc of sixty degrees. The object is projected 
simultaneously through the six slits onto the detector ring for each position of the 
aperture, such that 6n projections are collected at 60/n degree increments. The resulting 
series of two dimensional scintigrams must be processed to reconstruct the transverse 
section of the original source distribution. The solution to calculating a density function 
from its projections was first presented by Radon in 1917 in the form of an integral 
equation.’ Other approaches have since been presented including Fourier techniques and 
iterative solutions. This chapter will discuss the two general algorithms which have been 
implemented for image reconstruction with SPRINT: conventional filtered backprojection 
and statistical iterative reconstruction using the space alternating generalized expectation- 
maximization (SAGE) algorithm. 

Filtered Backprojection Algorithms 

Consider the simple case of a point source suspended in the center of an imaging 
system. Acquired projection images of the point map the location of the point in the 
direction parallel to the detector, however, the source could lie at any point along the line 
perpendicular to the profile. Since the depth of the activity is unknown, a first 
approximation of the source distribution is obtained by projecting the data from each 
projection back across an image field (figure 4.1). This simple operation is called 
backprojection. The resultant crude image of the point source is called a laminogram and 
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it demonstrates a build-up of intensity around the source creating a “star” effect. The star 
is caused by the complete projection of the point as a ray (figure 4.1 (B)). If an infinite 
number of views were acquired of this point, the backprojected image would demonstrate 
a blurring proportional to l/r, where r is the distance from the point source. If the system 
is linear and spatially invariant, it can be shown that the backprojected image is equal to 
the true image convolved with a l/r blur and to reconstruct the true image, this l/r 
blurring must be removed. The most common method for removing this blur is by 
filtering the Fourier transform of the individual projections, and then backprojecting the 
filtered projections. To present why and how this is accomplished, the central section 
theorem is introduced^. 



Figure 4.1. (A) Acquisition of Point Source Projection at First Detector Location. 
(B) Backprojected Distribution Of Source Into All Picture Elements Along 
Acquired Projection’s Representative Angle. Cross section of source distribution is 
represented by the criss-cross of backprojected ray-sums. 
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The Central Section Theorem and Parallel-Beam Filtered Backprojection 

Given a projection in the x direction of the source distribution f(x,y), the resulting 
projection will be g(y), such that: 

giy)= (4.1) 

This projection represents an array of line integrals in the x direction. According 
to the central section theorem, the Fourier transform of a parallel projection of a two- 
dimensional function f(x,y) with respect to a radial distance R taken at angle 9 gives a 
slice of a two-dimensional Fourier transform of the same function F(u,v). subtending an 
angle 0 with the u axis in the Fourier domain. This is depicted in figure 4.2. Since the 
transform of each projection forms a radial line in F(u,v), we can fill F(u,v) by acquiring 
many projections of the object and taking their transforms. Once filled, F(u,v) can be 
inverse transformed to reconstruct the original distribution f(x,y). The process can be 
illuminated further using figure 4.3. A source distribution f(x,y) is observed by a camera, 
forming an array of line integrals, each being a distance R from the origin, with the 
camera rotated at angle 0. This forms a projection given by: 



Figure 4.2. Illustration of the Central Section Theorem. 
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Pg{R) = I j/(x,>')J(xcos 6 ’ + >’sin 6 >- R)dxdy 


(4.2) 


where P 0 (R) is the projection information in the 0 direction. The integration takes place 
along line x cos0 + y sin0 = R or, in polar coordinates, r, r cos(0-(t))=R. The delta line 
in the integrand sifts out the desired line from f(x,y), forming an effective line integral. 
We recognize Pe(R) to be a series of one dimensional projection measurements taken as a 
function of R at different angle 0. 



Figure 4.3. Projection of a Two-Dimensional Function 

The derivation of the central section theorem is as follows. The two-dimensional 


Fourier transform of a source distribution is given as: 
F(u,v) = J|/(. ^yy<2.(ux^'0')dxdy 


(4.3) 


This can be expressed in polar coordinates F(u,v)=F(p,0), where u=p cosj3 and v=p sinyff: 
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F{p,f3) = 


(4.4) 


Alternatively, if the Fourier transform of the original projection is taken with respect to R 
at a given angle 0: 

FM = 


= \P,(R)e-‘'‘’^dR 

“00 

= J[ 11/(^! y)^ix cos J3 + y sinp- R)dxdy^e~'^’^^’^dR 


(4.5) 


It is clear that equations 4.4 and 4.5 are equivalent, and therefore the central limit 
theorem can be stated as: 

F(p,/3)^__,^:5,{P,iR)} (4.6) 


Thus the Fourier transform of a projection at angle 0 forms a line in the two- 
dimensional Fourier plane at the same angle, such that P=0. 

The central section theorem indicates that after filling the entire F(p,0) space with 
transforms of the projections at all angles, the reconstructed source distribution is 
provided by the two-dimensional inverse Fourier transform. In practice, a two- 
dimensional Fourier transform is computationally demanding £uid requires various 
interpolations and coordinate transformations. Further, since the Fourier domain is in 
polar coordinates, we observe that the information density is more sparse in the periphery 
of the Fourier image than at the center. This introduces larger interpolation errors in the 
high frequency aspects of the image, resulting in image degradation. To overcome these 
difficulties, and to resolve the original problem associated with the 1/r blur, filtered 
backprojection is used. The derivation of this algorithm begins from the previously 
derived result that the inverse Fourier transform of the Fourier image produces the desired 
source distribution: 

= J lF(u,v)exp^i27t(iix + vy)^dudv 


(4.7) 
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Conversion from Cartesian coordinates into polar coordinates of the frequency domain, 
where u=p cos 0, v=p sin 0, and dudv = p dpd0, allows the following derivation: 


2;roo 

/(jc,y) = J I F{p,6) exp[/2;r/?(xcos 9 + ysinO)\pdpdO 
0 0 
TTCO 

= IJ F{p, 6) exp[/2;7p(x cos ^ + y sin 9)\pdpd9 + 
00 


JI F{p, 9+n) exp[/2;T/?(x cos(^ + n')+y sin(^ + n'))]pdpci9 
00 


m 00 


= J 
0 


j F(p,6)\p\ exp[/2;ry9(x cos9 + y sin 6)\dp u6 


-00 


n\ 00 


f{x,y) = J j J 

oi-ooL 


\F(.p.e^e‘^’*dp 


ilrfd 


d{x cos ^ sin 0 - R)dR^0 


(4.8) 


We observe that this solution composes three operations. The term in the inner 
bracket represents a filtering operation in the frequency domain which removes the 1/r 
blur. It was earlier presented that the backprojected image equaled the true image 
convolved with a 1/r blur. Based on the convolution theorem, the Fourier transform of 
the backprojected image equals the Fourier transform of the true image times the Fourier 
transform of 1/r. It can be shown that the Fourier transform of 1/r is equal to 1/p in the 
frequency domain. Thus the Fourier transform of the true image can be obtained by 
multiplying the Fourier transform of the backprojected image by p: 


Backprojection = True * (1 / r) 

^{Baclqjrojection) = 3{True) x 3(1 /r) (4.9) 

3{True) = ^{Backprojection) x |p| 

From the central section theorem, this |p| filter can also be applied to each Fourier 
transformed projection to achieve the same filtering effect. 




55 


The next operation in (4.8) represents the backprojection operation of each filtered 

projection onto the x,y space along the line x cos0 + y sin6 = R. 

00 

|g(i?)t^(xcos^ + ysin^-f?)i//? (4.10) 

—00 

The final integrand represents the backprojection operation for all angles 0 over the range 
from 0 to 71. Thus the three consecutive operations are filtering, backprojection, and 
moving to the next projection to repeat the process. These operations are carried out in 
the polar frequency domain, permitting a simple filtering operation and preservation of 
high frequency aspects of the image. These three operations can be better seen by re¬ 
writing equation 4.8 as: 


fix,y)= J- 
0 
n 

V 

0 


CX) 


J 

—OOL“°o 




S(x cos 0 + ysin6- R)dR 


\de 


j 3/j^3j|P^(i?)||yC>|J<^(xcos^ + y sin^- 


(4.11) 


—00 


Here it is clearly observed that the function backprojected at all angles is a filtered 
version of the projection Pe(R), where the filter provides |p| weighting. 


Convolution Backprojection 

The back-projection function observed in (4.11) can be written as: 

again using the convolution theorem of Fourier transforms. This operation indicates that 
the multiplying by |p| in the frequency domain is equivalent to convolving the inverse 
transform of |p| in the space domain. This introduces the convolution back-projection 
algorithm, where each projection Pe(R) is convolved with a function c(R) and then 
immediately backprojected. The function c(R) is chosen to correct the 1/r blur and 
provide any additional filtering which may be required. The convolution-backprojection 
algorithm thus becomes: 
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1 


00 

J[P^(i?)*c(/?)]«5(x cosd + y sind- R)dR 


0 


[—00 


\d0 


Converging and Slit-Aperture Filtered Backprojection 


(4.13) 


The derivation just presented is suitable for cases where the line integrals are 
perpendicular to the surface of the detector, as when a parallel hole collimator is used 
with a gamma camera. In contrast, the SPRINT system makes use of a slit aperture which 
creates a two dimensional pinhole geometry with equally spaced rays and is similar to 
fan-beam type geometries. Figure 4.4 depicts the difference in detected rays experienced 
with parallel, fan-beam (converging) and pin-hole (slit) type collimators. In the fan-beam 
and pin-hole geometries we observe that while the field-of-view is diminished, the 
projection of the object onto the detector is magnified, assuming the object is sufficiently 
close to the collimator. This results in both increased sensitivity and resolution for the 
smaller objects being imaged. There are two specific categories of converging/pinhole 
systems for which simple filtered back-projection algorithms exist: one where the rays 
measured at the detector are equally spaced in angle, the other where they are equally 
spaced in distance on the detector. As an example of each, positron emission tomography 
systems normally have detectors equally spaced in angle. For conventional gamma 
cameras equipped with fan-beam collimators and for the case of SPRINT where a slit 
aperture is combined with a detector that is segmented into equally spaced bars, the rays 
are equidist£int. This is observed in figure 4.5. The following derivation for fan-beam 
reconstruction follows that of Zeng and Gullberg^ modified for application to SPRINT. 
Additional information can be found in the work by Berthold and Horn'*. 
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Pinhole (Slit) 

Parallel Collimator Fan-Beam Collimator Collimator 



Figure 4.4. Comparison of Field-of-Views and Ray Paths for Three Collimators 


Pinhole (Slit) 

Fan-Beam Collimator Collimator with 

Segmented Nal Detector 



I I ' \ 

/ I '\ \ \ 


Figure 4.5. Comparison of Fan-Beam and Slit Aperture Showing Equidistant Rays 

(not to scale). 

Consider the fan beam geometry presented in figure 4.6. Let the center of rotation 
be at distance D from the detector, and the focal length of the collimator is f The fan 
beam is measured at equal intervals along the camera face, s, with the camera at angle p 
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resulting in projection Rp(S). The rays of the fan beam are specified at angle a. In figure 
4.6 B, we consider the projection being mapped to a parallel coordinate system along 
Dr-D2’ which passes through the center of rotation. We note a projection Py(t) in this 
new frame is equivalent to Rp(S) in the original frame. From figure 4.6 B, the following 
relationships are noted and used in the derivation; 

f S 

X = a + p = p + arctan - 


\D. 


t = Dsma = Scosa = 


SD 


dtdX = 




dSdp 


(4.14) 


From the parallel case of convolution-backprojection, let the source distribution be 


represented by the following: 


n 


fix,y)=\ 

0 


I [r^ (S)* cosp + y sin p - 


(4.15) 


Substituting the relationship Py(t) for Rp(s) and expending the convolution: 


n 


f(x,y) =f 

0 


00 

\^P^{t)c{t- T)^5{xco^y + ysiny - t)dt 


—00 


\dy 


n 


0 


00 


(4.16) 


I P {t)c(x cosy +ysmy - r)dT\dy 


—00 


The conversion to polar coordinates can be made by substituting x = r cos 6 and y =r sin 0 


71 


/M)=f 

0 


00 


j P (r)c(r cos(x -6)- t)dr my 


(4.17) 


-00 


Consider reconstructing over an angular range from 0 to In for this geometry. Similarly, 
the radial range is limited from -tmto tm. This converts 4.17 to: 


1 


Itz 


^ 0 


j P (r)c(r cos(^ -y)- '^)dT 


-C 


\dy 


(4.18) 
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- s„ 

2;T-arctan— r 

1 D f ns ^ S 

fir,e) = - I j I P,,^^-^ 2 -^Jc(rcos(^ + arctan—-^) 

-arctan— 

D 


dS)dp 
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r coa p-9 + arctan 


D 


(R a\ ^ 
= r cos(p - 6) 






(D + rsin(/?-^))^ 
(S'-S)UD 


+ S^ 


Substituting these variables and the new filter function result back into 4.20: 


1 


2^ 


^ 0 




-S^ 


(S-S')- 


UD 




dS 




Let: 


/>' = 


pUD 




This reduces the filter function to: 




+S 

D^+S^ 


2 « 

- i\p'\^ 


i2)r(S'-S)p'j I 


dp' 


c{S'-S) 


(4.22) 


dp (4.23) 


(4.24) 


(4.25) 


The resulting fan-beam convolution-backprojection reconstruction algorithm resolves to: 


In 


/('•,«)= j 


00 


D 


0 ^ +5'^) 




ic(5-5') dS 


\dp 


(4.26) 


Equivalently, the filtered backprojection form of the algorithm becomes: 


2n 


f(x,y)= J 

0 




00 

\ 


|, ^ „Rf(S) 


S(x cos 9 + y sinO- S)dS 

—CX) 


J 

2 



\dp 


(4.27) 

Comparing the parallel beam filtered backprojection algorithm (4.11) with (4.27), 
we make the following observations. For fan-beam reconstruction, the original projection 
data is weighted by the factor cos a, equivalent to D/Vd^+S^ . Also, the ramp filter is 1/2 
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that used for parallel reconstruction. Finally, the backprojection is performed after 
weighting each filtered fan beam projection by 1/U^ where U is a function of r, p, and 0. 

Slit (Pinhole) Reconstruction Factors for SPRINT 

The slit aperture of SPRINT introduces a requirement to correct measured 
projection data for the effects of angular incidence onto the slit as well as the detector^. 
Consider figure 4.8 A, which demonstrates the fluence of photons emitted from a point 
source which then impinge on a slit. As the angle of incidence increases, the effective 
width of the slit decreases by the factor cos 0. This causes a reduction in the measured 
flux which must be corrected by dividing the measured projection data by cos 0. 

A similar effect is observed in figure 4.8B, which shows the obliquity factor 
introduced as a given flux impinges on the detector. As rays traverse from their source to 
the detector, there is an increase in the length of the ray path as the angle of incidence 
increases. In conventional pin-hole collimated systems this introduces an inverse square 
falloff with distance, and causes a cos^0 dependence with incidence angle. In contrast, 
SPRINT maintains a 2-D geometry defined by the joint response of slit aperture and slice 
collimator. The subsequent decrease in projection intensity is proportional to 1/r 
introducing a cos 0 dependence. The rays upon arriving at the detector are again incident 
at a given angle to its surface. As this angle increases, there is an increased loss of 
intensity as the fluence spreads out to strike the detector surface. This effect introduces 
an additional cos 0 dependence. The net effect is a cos \\i cos^0 correction required for 
the measured projection data. 
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Figure 4.8. Obliquity Factors Effecting Measured Fluence Rates 
The Filter in Filtered Backprojection 

The p filter has a simple ramp shape in the frequency domain as shown in figure 
4.9. Because of its shape it is frequently described as a ramp filter. It is not necessary nor 
desirable to have the filter continue on to very high frequencies since useful information 
from the image is only as high as the system geometric resolution and angular sampling 
permits . The cut-off frequency v, applied to the ramp filter is determined from sampling 
theorem^ which states that to recover spatial frequencies in a signal up to a maximum 
frequency Vmax, the signal must be sampled with a mininum linear distance d: 

d<-^ (4.28) 

This means that the highest spatial frequency which can be recovered from the data must 
be sampled at a m inimum of two points per cycle. Coarser sampling does not allow 
higher frequencies to be recovered accurately and the result is aliasing artifacts in the 
reconstructed image. Linear sampling in SPRINT is once per 3 mm in the image plane, 
which equates to 1 m m in the object plane due to system magnification. This limits the 
resolvable frequency to 0.5/mm, termed the Nyquist frequency of the system. The cut-off 
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frequency of the ramp filter v, is chosen to be greater than the Nyquist frequency to insure 
the image maintains high frequency components. A cutoff frequency of 1/mm is 



Frequency, v 


Figure 4.9. Ramp Filter to Correct for 1/r Blur in Frequency Domain 

Since the p filter is a high-pass filter, high frequency noise in the image will be 
amplified by its application up to the cutoff frequency. To reduce this effect, low-pass 
filters are typically applied in conjunction with the ramp filter to reduce this high end 
noise. Typical examples of low pass filters include Hamming, Hanning, Butterworth and 
Shepp-Logan. Figure 4.10 demonstrates the effect of applying the Hamming and 
Hanning filters with the ramp filter. Here we observe the effect of the ramp filter is to 
exactly suppress the 1/r blur and very low frequency background components of the 
image. In contrast, the Hanning/Hamming filters suppress high frequency noise resulting 
from the statistical variation in counts experience between neighboring pixels. The filters 
are implemented through simple multiplication of their Fourier transform with the 
transformed projection data in the frequency domain. 





Amplitude 
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Application of Low-pass Hanning and Hamming Filters with High-Pass 

Ramp Filter 
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Figure 4.10. Combined High-Pass Ramp Filtering with Low-Pass Hamming or 

Hanning Filters 
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STATISTICAL RECONSTRUCTION ALGORITHMS 


Filtered backprojection techniques are widely used for image reconstruction 
because of their relative simplicity to implement as well as the speed in which they can 
reconstruct a fairly good quality image. Despite these benefits, there are potentially 
significant advantages to using discrete estimation techniques for image reconstruction 
that are based on a statistically defined objective function. Iterative algorithms can 
readily include an accurate model of the imaging system which may include a spatially 
variant geometry, and the physical effects of scatter and attenuation. Through a 
geometrically and electronically accurate system model, it becomes possible to 
“deconvolve” geometric resolution to some degree. By incorporating biasing factors in 
the model (i.e. background radiation, etc.) the accuracy of the reconstruction can be 
improved. Further, these methods do not suffer from the uniform sampling requirements 
imposed by Fourier methods making them more capable of handling missing projection 
data. Because the reconstruction algorithm addresses the statistics of photon counting, 
reconstructed images can have more favorable noise characteristics than achieved with 
FBP. This is particular important in SPECT, where projection data typically 
demonstrates high statistical variance from low counts due to the limited radioactivity 
used and the wide variance in statistical accuracy of different projection elements. The 
chief shortcoming of these algebraic techniques which impedes their widespread adoption 
is that they are iterative, causing them to be much slower and computationally more 
demanding then Fourier methods. 

In matrix notation, the issue is to solve the inverse problem: 

Y=AX (4.29) 

where Y is the set of measured projection data, and X is the source distribution to be 
determined. A is a matrix of weights, computed for a given imaging system geometry, 
the elements of which describe the mean contribution of each pixel to each ray-sum 
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measured by a particular detector element. These weights reflect the probability of 
detecting a photon emitted from a given source voxel by a given detector element. The 
weights matrix is readily adapted to include such factors as the solid angle of the detector, 
detector spatial resolution, attenuation and scatter within the object, cosmic background 
and the collimator response function. 

Maximum Liklihood Estimation and the Expectation Maximization (EM) 
Algorithm 

A specific approach to solving this inverse problem is based on the Poisson nature 
of the projection data, in which liklihood is the objective function to be maximized. The 

likelihood relationship is written: 

i(l) = log likelihood(7|l) (4.30) 

The most common class of algorithm used to resolve this problem are termed 
expectation-maximization (EM) algorithms; using iterative techniques to find the 
maximum likelihood estimate. 

Lange and Carson^ present the original derivation of this algorithm. The observed 

data in an experiment can be represented as a random vector Y. Let Y have the density 

function (or likelihood) Ly(Y, A,), where X is some vector of parameters to be estimated. 

In general, it is difficult to maximize Ly(Y, X) with respect to X. A possible solution to 

this problem is to conceive of a richer or larger sample space for Y where optimization 

problems are easier to resolve. Thus, the EM algorithm postulates a “complete data” 

random vector “X” such that Y is a function of h(X). Furthermore, X is assumed to have 

a density function Lx(X, A.) with respect to a measure |a(X). Using these assumptions, 

Lx(Y,A-) can be recovered by integration: 

L^{Y,A)= (4.31) 

All density functions under consideration are discrete. This reduces the above expression 
to: 




68 


L,(Y,A) = '£L,<,X,A) (4.32) 

Each iteration of the EM algorithm consists of two steps. In the expectation (E) 

step, one performs the conditional expectation: 

E(lnLJX,A)\Y,r) (4.33) 

where A," denotes the current vector of parameter estimates. In the maximization (M) 
step, this conditional expectation is maximized with respect to X, X" being held constant, 
to give the new estimate of The essence of the EM algorithm can be explained by 
considering the function: 

A ^ ) = .E(ln 4 (X, A)\y, 1”) - In I/F, 1) (4.34) 

Regardless of the value of X", H(X, X") attains its maximum when X = X". The EM 
algori thm exploits this fact by choosing X = X"'^' to maximize the relationship in 4.33. It 
follows that: 

ln4(T,l”^')>In4(r,/l) (4.35) 

In other words, the EM algorithm is designed to increase the log liklihood at each 
iteration. 

Properties of the Poisson Distribution 

A non-negative, integer-valued random variable Z follows a Poisson distribution 
if: 

P{Z = k) = e-^^ (4.36) 

k\ 

for some p > 0. The Poisson distribution has the following characteristics: 

1. Z has the mean E(Z) = p. 

2. If Zi, Z 2 ,... Zm are independent Poisson variables, then Z= ZZj is also a 
Poisson variable. 

3. Suppose that a random number of Z particles is created by a Poisson process 
with mean p. Let each particle be independently distributed to one of m categories, the 
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kth of which occurs with probability pk- If we let Zk represent the number of particles 
which enter category k, Zk follows a Poisson distribution with mean ppk. Furthermore, 
the random variables, Zi,... Zm are independent. 

The relevance of these characteristics is directly applicable to SPECT imaging, 
where the number of photons recorded in a given element of the detector module are 
independent random variables and the number of photons recorded for the various 
projection measurements are independent and random. Because of this independence 
property, the likelihood over all projections reduces to the product of the separate 
projection likelihoods. 

The EM Algorithm in Emission Computed Tomography 

Consider the SPRINT II system depicted in figure 3.3. The field of view is 

divided into small pixels in a Cartesian grid. Let Zj represent the source intensity in pixel 

j, and cy be the probability that a photon leaving pixel j reaches the ith projection 

detector. Assume the emission source is viewed by N detectors, and let Ny denote the 

number of emissions from the jth pixel that are detected by the ith detector. From the 

Poisson nature of the source distribution, the mean number of detected photons 

originating from pixel j during the ith projection is: 

= A,jAj (4.37) 

and 

Ny ~ Poisson^AyAj^ 

Here Atj is the time interval for collecting a projection, and Ay are considered known 
constants. The intensities Xj are to be estimated, which we will abbreviate as X. 

Define I, as the set of pixels which contribute to projection I, and Jj as the set of 
projections to which pixel j contributes. Now let Xy, j eh, be the random number of 
photons that are emitted by pixel j and contribute to projection i. As mentioned above. 
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the mean of Xjj is AijX,j. Also let Yi be the total number of photons recorded for projection 
i. From these assumptions: 

Poissonj^ AyZj | (4.38) 

In view of the characteristics of the Poisson distribution, Xy and Yj have Poisson 
distributions. Yj being the observed data for projection i, Xy constitutes the unobserved 
but complete data set. This equation represents our system model. 

To specify the E-step of the EM algorithm, we need the complete data log 
likelihood: 

lnL,iXa) = ln( 41 ,)-lnX,!) (4.39) 

/ ye/,. 

Taking the conditional expectation of Xy with respect to Yj and the current vector of 
parameter estimates A," yields: 

E[\viL,{X,X)\Y,X^) = SZ {-4^y + 4 ln(4^y)} + ^ (4-40) 

;• ye/, 


where: 

n„ = e[x,\x.x-) 




(4.41) 


and R does not depend on the new A. 

The M-step is now trivial. Partial differentiation of equation 4.40 with respect to 
Aj yields: 


dX 


■E{\nL^{X,X)\Y,X) 


/&/y /e/y 


which is set equal to zero, and solved: 


T.N, 


'in+\ _ 

‘ “XT 


X” 


-z 


AjY, 




(4.42) 


/e/y /eJy keli 

for the jth component of the new parameter vector A"'*''. This expression has a fairly 
simple intuitive explanation. Ny is the best current estimate of the number of detected 
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photons that originate from pixel j during the ith projection. As above, the sum ^ N-j is 

iejj 

the estimated number that originate from pixel j summed over all projection that pass 
through pixel j. This estimated number is set equal to the expected number ^ AyAj and 

ieJj 

the resulting equation solved for A.j. 

Equation 4.42 has the following additional benefits. First, the non-negativity 
constraint, Xj" > 0 is automatically satisfied, assuming Xj® is positive. Second, each 
iteration verifies that: 

/ j el i i 

This indicates that the current expected number of photons stays fixed at the total 
observed number, thus preserving counts. 

The Space-Alternating Generalized EM Algorithm for Penalized Maximum- 
Likelihood Image Reconstruction 

EM algorithms have two important drawbacks which prevent their widespread use 
for routine image reconstruction. The first problem is that EM algorithms converge very 
slowly, particularly when random background events such as scatter, dark current, and 
background radiation are included in the algorithm. Second, because of the blurring 
effect of photon detection and position estimation the reconstruction process is ill- 
conditioned. Consequently, images that are reconstructed purely by using the maximum 
likelihood criterion are unacceptably noisy. To overcome this, a smoothness penalty or 
regularization parameter must be introduced into the objective fimction. Unfortunately, 
the M-step in the EM algorithm becomes intractable when one includes smoothness 
penalties, making its application difficult. 

n 

To overcome these issues, Fessler and Hero developed the space-alternating 
generalized EM (SAGE) algorithm which attempts to exploit the structure of the log- 
likelihood to hasten convergence and allow regularization. By comparison, the EM 
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algorithm simultaneously updates all pixels in the image at each iteration, whereas the 
SAGE algorithm implements a sequential pixel update. This method eliminates coupling 
between pixels at this step and enables use of local smoothness penalties. Further, the 
slow convergence of EM algorithms is tied to the massive size of their complete-data 
spaces^. The sequential update methodology of SAGE allows the use of small “hidden- 
data” spaces for image reconstruction, which leads to fast convergence. The actual 
details of the SAGE algorithm are complex, and beyond the scope of this work. The 
following briefly summarizes the salient features of SAGE as compared to the EM 
algorithm. 

We begin with a more descriptive system model than presented in equation 4.38 
for the EM case. Specifically: 

i; = X ^ ~ Poissonj^ A,jXj -h r, | (4.43) 

where p represents the cumulative random background rates stemming from dark current, 
background radiation, scatter, and aperture/collimator penetration. The conditional 
expectation of the log-likelihood is of the form: 

E{\nL,(X\X'\Y,l’)-= +z)+N„ ln{Aj{Aj +Zj)'j-fiP(A) 

(4.44) 

where a penalty term is subtracted from the expectation, and zj is introduced to 
incorporate the random background. When comparing to equation 4.40, the expectation 
step is conducted over a “hidden” but complete data space X'’ which includes all ray-sums 
that pass through the pixel to intersect detector element i. This isolates (or decouples) the 
pixel being evaluated from the remainder of the estimated image permitting faster 
convergence and inclusion of a smoothness term. The expectation result is reduced by a 
quadratic smoothness penalty, PP(A,), that weights the pixel currently being evaluated by 
the difference in values between it and the eight pixels which neighbor it. The larger this 


z, = tnm^, 

A. 
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difference is, the larger the penalty. The factor p can he varied to achieve different levels 
of smoothing in the reconstructed image. A factor Zj is added to the pixel estimate, A,j, 
which provides inclusion of the random background. In contrast to conventional EM, 
there is no summation over j and i, as this expression is maximized and evaluated 
sequentially for each pixel. The M-step is performed by taking the partial derivative of 
the expectation with respect to Xj, resulting in a quadratic expression from which is 
calculated. The projection estimate, Ny, is then updated before proceeding to the next 
pixel, at which point the iteration loop is repeated. A more detailed presentation of the 
SAGE algorithm derivation can be found in the reference by Fessler.*® 

Generation of the SPRINT System Weights, A 

The accuracy of the reconstructed image from any iterative algorithm, SAGE 
included, is directly dependent on the accuracy of the model used to determine the system 
weights matrix. Specifically, the geometry of the system must be well defined, and the 
detection efficiency, or relative detection efficiency of each detector element must be 
known. The system geometry of SPRINT is complex due to the combined geometry 
defined by the slit aperture and slice collimator, the obliquity factors of the aperture and 
detector, and the effect of rays which intersect partial voxels or partial detector elements. 
We will briefly examine each of these issues, and discuss their handling before 
concluding. 

To simplify the determination of system weights, SPRINT is considered a two- 
dimensional geometry, with weights generated for a single image slice and then assumed 
applicable for all slices. This is a good approximation since the slice collimation defines 
a consistent geometry axially. The size of the weights matrix is determined by the 
number of projections, the number of detector elements, and the number of image pixels. 
The SPRINT gantry contains 488 three mm wide Nal(Tl) bars in its 11 modules, of which 
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459 are useable for high energy imaging'. A single projection measurement consists of 
the recorded counts from approximately one-sixth of these bars, with six projections 
being acquired at each position of the slit aperture. The total number of projections 
required is derived from the sampling requirements necessary to obtain the desired image 
resolution. Bracewell*' has shown that the number of angular samples per 360° required 


to give uniform resolution and freedom from angular aliasing is approximately: 
JV = ^ 


(4.45) 


where D is the diameter of the aperture ring equal to 13 cm, and rs is the desired system 
resolution. If rs is taken as 4 mm, which is somewhat less than the system resolution 
calculated in Chapter III, then approximately 200 projections would be required, 
necessitating 34 steps of the aperture (200/6). For simplicity, a forty-step acquisition 
protocol is thus used, which generates 240 projections. Since there are 459 measured 
elements in every set of six projections, there are 18,360 rows in the weights matrix 
which describe the probability that a photon emitted from a given voxel is recorded by a 
given element. 

The nine centimeter wide field of view of SPRINT is segmented into a 100x100 
pixel field, each pixel then being 0.88 mm on a side. The voxel depth (slice thickness) 
was chosen to be 3 mm, equal to the pixel size in the projection data. Each column of the 
weights matrix is associated with a given pixel, so the final weights matrix is large, 

18,360 rows by 10,000 columns. The geometric detection probability expressed by each 
weight is to the solid angle subtended by a given detector pixel as seen by each image 
pixel. The probability is approximately equal to: 



‘ This will be elaborated upon in Chapter VII. 
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where 1 is the effective width and h the effective height of the detector pixel observed, 
taken as 3 mm (refer to figure 3.3). The distance d is that from the image pixel to the 
detector pixel. Weights are determined sequentially for a given detector element and all 
image pixels which can be “seen” by that element. Detector pixel centroids were 
determined from the detector response measured by rotating a highly collimated slit 
source around the detector ring. Primary image pixels include those intersected by a ray 
passing through the center of the detector pixel and the center of the slit. Secondary 
image pixels lie off-axis of the central ray and contribute (at least partially) to a given 
detector element response. Dependent on the image pixel location, the effective detector 
width, 1, seen by the pixel will be greater than, equal to or less than the actual element 
width. The detector response is weighted accordingly to reflect the effective pixel width 
seen by an image pixel (figure 4.11). 



Figure 4.11. Weight Generation for Image Voxels Observed by a Detector Element 
(Solid Lines). Effective detector width, 1, determined from part of detector element 

observed by image pixel (dotted line). 

Weights are stored in the ASPIRE format to take advantage of the sparse data 
structure, and ASPIRE routines were used to perform vector operations^ 
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CHAPTER V 


DESIGN, CONSTRUCTION AND PERFORMANCE OF ONE-DIMENSIONAL 
LEAD FOIL PARALLEL PLATE AND FAN-BEAM COLLIMATORS FOR 

SPRINT 


SPRINT imaging requires construction of precision slice collimators which are 
specific to the gamma-ray energy being imaged and the desired object (organ) size to be 
examined. The first requirement then for SPRINT imaging of 511 keV photons is 
establishing a system axial resolution which is significantly less than rat organ sizes of 
interest. This chapter describes the design criteria and methods developed for 
constructing parallel and converging ring collimation for SPRINT and ASP* The 
measured performance of the constructed collimators is also presented. 

Design criteria and Early Construction Methods 

As described previously, collimation of SPRINT detector modules is provided by 
a rotating multi-slit ring aperture which collimates radiation azimuthally combined with a 
stationary slice collimator which provides axial collimation. Slice collimators are 
composed of alternating annuli lead foil and air equivalent foam, and can have either 
parallel or converging designs. SPRINT was designed for human brain, primate brain 
and small canine brain imaging of ^^""Tc labeled compounds at 8 mm resolution and 
tumor imaging of *^*I labeled compounds in small animals at 3.5 mm resolution. This 
work considers biodistribution assessment of positron labeled radiotracers in the rats, 
requiring a 5 mm resolution. Each application requires a unique collimator to achieve the 
specified resolution. To image 140 keV ^^"’Tc gamma rays and 511 keV annihilation 
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radiation, two parallel plate collimators were constructed, respectively. A converging 
collimator with 2:1 magnification was constructed to achieve the very small resolution 
desired for imaging ^^*1. Designs were developed for ^^"’Tc and collimators using 
ray-tracing techniques. Monte Carlo simulations were used to design the 511 keV 
collimator and confirm the design for the converging collimator^. 

Figure 5.1 shows the comparative performace of several collimator designs 
simulated for 511 keV imaging. The designs are specified hy the thickness of lead and air 
equivalent spacer in thousandths of an inch. The three sets of bars show the relative 
sensitivity between designs obtained by summing over the response fimctioon FWHM 
and FWTM, as well as the peak amplitude of the response. A tradeoff is observed up to 
the 60x80 or 30x60 designs when either increasing lead foil thickness or decreasing the 
spacer thickness leads to an increase in resolution but decrease in sensitivity. Other 
considerations in design selection included spatial variance of the response function and 
the availability of spacer material in desired widths. 
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Figure 5.1. Design comparison of collimator performance from Monte Carlo 
simulations. Each design simulated in the x-axis shows the lead foil thickness (mils) 

by the air gap thickness (mils). 
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The final designs chosen were a compromise between resolution, sensitivity and 
spatial variance. The low energy collimator used 0.23 mm lead foil and 2.0 mm spacers 
with a septal length of 6.8 cm. The medium energy converging collimator required 
tapered slices, with 0.50 mm lead foil and 0.9 mm spacers at the collimator face with an 
11.6 cm length. Finally, the high energy collimator used 0.76 mm lead foil and 1.5 mm 
spacers, again with a length of 11.6 cm. Schematic cross-sections of the low-energy 
collimator and I converging collimators are shown in figure 5.2. The high energy 
collimator is similar in design to the low energy collimator, but with exterior dimensions 
and septal length equal to the converging collimator. The figure emphasizes the major 
difference in construction patterns and lead foil spacing for the two collimator types. The 
47.4 cm focal length of the converging collimator restricts the field of view to half that of 
the parallel plate collimators, but provides a 2-to-l axial magnification factor. 



Figure 5.2. Cross-Sectional Views of Low Energy Parallel and Converging 
Collimators. Parallel designs provide a 15 cm axial field of view (FOV). 
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Construction of these collimators proved to be very challenging. Tolerances could 
not exceed 50 microns between individual slices to obtain adequate azimuthal and axial 
uniformity. Early efforts at construction used Styrofoam sheet bonded to pure lead foil 
using epoxy compounds. Styrofoam is not readily machinable nor can Styrofoam sheet 
be accurately molded, so that is was virtually impossible to produce accurate spacers of 
the desired thickness. Although expanded Styrofoam beads can be molded with high 
precision, the process is difficult and was not attempted. Attempts were made to mold 
expanding urethane foam to the desired width and taper. Severe difficulties were 
encountered in releasing the material from the mold and maintaining dimensional 
tolerance and uniform density. 

Other complications were encountered using pure lead and liquid adhesives. Pure 
lead foil is highly malleable and too easily wrinkled or dented during the fabrication 
process. The bond strength of epoxy resins for both Styrofoam and lead foil was ver 
poor, and it was impossible to maintain high tolerances on epoxy thickness. 

These problems were overcome through the careful selection of materials and 
design of a precisely controlled construction process. The method could repeatedly 
produce either flat and tapered collimator slices with a 50 micron tolerance over a wide 
range of foam and lead thicknesses. 

Collimator Materials 

The foam product chosen was a rigid polymethacrylimide foam, trade name 
Rohacell, which had a number of highly desirable mechanical characteristics including 
high dimensional stability and strength, and excellent bonding and radiation transmission 
characteristics'*. The foam can be ordered in sheets over a wide range of thickness’ 
appropriate for collimator construction, and is readily cut and machined to tolerances of 
50 microns. Thicknesses of foam sheets obtained from the supplier had a tolerance of 0.1 
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mm . The foam density was 0.052 g/cm^ with linear attenuation coefficients of 0.0076 
cm'*, 0.0055 cm'*, and 0.0048 cm'* for 140 keV, 364 keV and 511 keV gamma rays, 
respectively. This results in an approximate 5% scatter and attenuation probability for all 
three collimators. 

The high malleability of pure foil was reduced by purchasing a lead alloy with 2% 
antimony and 0.5% tin by weight. This foil was stiff enough that sheets could be handled 
and cut to size without any significant deformation of the surface or edges. The lead 
remained sufficiently malleable that any surface contours were readily flattened during 
the bonding process. Foil could be ordered at the desired width and precise thickness 
(<25 micron tolerance) required for a given collimator design. 

DOW DAF 899 thin film adhesive was used to bond the lead foil and foam^. This 
dry, tack-free adhesive is a hot-melt film which formed a strong bond to both lead and 
foam at a melting temperature of 100° C. The product has a nominal thickness of 0.1 
mm with a tolerance of ±0.025 mm. After bonding, joint thickness was approximately 
0.05 mm. 


Collimator Contraction Methods 


Parallel Collimators 

Foam sheets were vacuumed to provide a clean surface for bonding. The sheets 
were precut using a square template slightly larger than the diameter of the finished 
collimator. Thin film adhesive squares were precut 2 mm smaller than the foam to 
prevent the adhesive from bonding to a heated press. 

Lead foil was dispensed from a roller assembly onto a smooth PVC cutting 
surface with extreme care taken to avoid deforming the lead. The lead foil was cut to the 
size of the foam sheets using a hand-held rotary fabric cutter with tungsten blades (figure 
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5.3). The precut thin film adhesive and polymethacrylimide foam were carefully stacked 
and aligned on top of the lead and the stack was sandwiched between a pair of thin 
aluminum sheets for transportation and uniform compression. The assembly was moved 
directly to a heated press which was preheated to 100° C and compressed for thirty to 
forty five seconds depending on the thickness of lead used. The press with an assembly 
ready for compression is shown in figure 5.4. 

The bonded assembly was moved to a milling machine mounted vacuum chuck 
shown in figure 5.5 for final cutting. Accurate cutting of the inner and outer diameters of 
each slice was a critical step in this construction process. The method required circular 
patterns to be cut in each slice with 100 micron tolerance, without damaging or feathering 
the foil edges for a range of foil thickness from 0.23 to 0.76 mm. The vacuum chuck 
featured a surface that was flat to within 50 microns, and had a series of surface 
perforations through which a vacuum could be drawn. PVC rings were embedded in the 
surface at the specified inner and outer diameters of the collimator to provide a cutting 
surface. 

The lead-foam slice was placed on the vacuum chuck, lead-side down, and a 
vacuum applied to hold the piece firmly in place. Slice radii were cut using tungsten 
blades from a rotary fabric cutter. The cutter assembly was mounted in the milling 
machine spindle as shown in figure 5.6. The cutter was lowered onto the PVC cutting 
surface using the spindle height adjustment, locked into place, and rotated by hand to 
make the cut. 

Finished slices were stacked on a surface table with two right angle fixtures used 
to assure precise slice alignment. Between each pair of slices a piece of thin film 
adhesive was added which was precut to the slice dimensions. The stack in progress for 
the low energy collimator is shown in figure 5.7. When the desired stack height was 
reached, aluminum endplates and compression plates were added and the entire assembly 
compressed using large C-clamps (Fgure 5.8). The required width of the finished 
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collimator was 15 cm. Each adhesive layer thickness decreased 50 microns during 
bonding, and this needed to be accounted for in order to obtain the proper finished 
dimensions. The resulting Tc-99m collimator had 59 slices, while the 511 keV collimator 
had 63 slices. Final bonding of the assembly was accomplished by placing the collimator 
in a large glass drying oven at a temperature of 100° C for several hours to allow uniform 
heating. To protect the exposed lead edges during handling, several thin coats of epoxy 
resin were applied to the exterior of the collimator. It was important that the epoxy used 
did not contain solvents that would attack the foam. 

Fan-Beam Collimators 

The design of a converging geometry slice collimator introduced additional 
complexity into the construction process because of the requirement to mill a radial taper 
in the foam surface of each slice. Each annulus required a 0.2331° taper to obtain 2 to 1 
axial object magnification from the finished collimator. Each slice was also a conical 
form in the stack. Consequently, each symmetric pair of slices had a unique inner and 
outer diameter, and required a precise azimuthal wedge to be removed from the annulus 
so that the slices would stack tightly together without stretching the lead-foam ring. 

Bonded lead-foam slices were placed lead side down on a radially tapered vacuum 
chuck which was mounted on a rotary table on the milling machine. The angle of taper, 
0.2331 °, was gradual enough so that bonded slices were held in place firmly by an 
applied vacuum. On one side of the rotary mill table, a tooling fixture was mounted to 
the milling table which would cut the inner and outer diameter of the slice, again using 
rotary cutters. Adjustments were available on this fixture for the angle of cut, the depth 
of cut, as well as the radius of the cuts. On the opposing side of the rotary table a fixture 
was mounted which allowed a slotted plate to be pivoted onto the slice surface. The slot 
provides a template to produce radial cuts in the slice so that a precise angle wedge could 
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be cut from the annulus. A fly-cutter was used to taper the foam and a hood assembly 
was used to the collect the foam dust generated. 

The taper was defined by the radial taper of the vacuum chuck so that the fly 
cutter need not be tilted. The slice was rotated azimuthally under the fly-cutter a full 
360° degrees using the rotary mill table, with a fast spindle speed and slow rotation of the 
table critical to obtaining a smooth uniform surface in the tapered foam. The fly cutting 
operation is shown in figure 5.9. The side-mounted rotary cutters were next used to cut 
the slice radii (figure 5.10). Finally, the wedge cut was made using a hand-held rotary 
cutter and the slit fixture discussed above (figure 5.11). 

The final dimensions of our fan-beam collimator required 78 total slices. Slices 
were constructed in symmetric pairs, each with unique radii and removed wedge angle, 
and stacked separately on two endplates machined with a 9.125° angle to match the top 
and bottom slice angles. The cut seam was sealed with a narrow strip of thin aluminum 
tape and then tack soldered to prevent the seam from splitting during compression and 
heating. Between each prepared slice a layer of thin film adhesive was placed which was 
also pattern cut to the same conical shape as the slice. For the center twenty slices of the 
collimator, the required radius change was less than a millimeter and the required wedge 
angle was less than 0.3 ° so that no wedge cut was required and the radii were held 
constant. When all slices were constructed, one stack of slices was inverted onto the 
other stack with right angle fixtures used to assure flush alignment of the slice stack. 
Finally, the top endplate and compression plates were added, and the entire stack 
compressed and bonded as described for the parallel collimators. The three completed 
collimators are shown in figure 5.12, with the fan-beam collimator shown inserted into 
the SPRINT-II gantry. 
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Figure 5.6. Mounted Rotary Cutters Shown with Completed Collimator Slice 
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Figure 5.8. Completed Collimator Stack Under Compression, Prepared For Heating 
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Figure 5.10. Side Mounted Rotary Cutters Making Inner and Outer Radii Cuts 
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Figure 5.11. Slice Prepared for Radial Cut Using Side Mounted Template 



Figure 5.12. Completed Collimators: High Energy, Converging and Low Energy 
from Left to Right, with Converging Collimator Mounted in SPRINT. 
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Performance Results 

Intrinsic Collimator Uniformity 

Evaluation of collimator uniformity is performed by collecting a detector 
uniformity correction map for each of the 11 detector modules with no collimator. Next, 
images of a uniform line source are acquired with each slice collimator. These are 
corrected using the detector correction map to obtain intrinsic collimator non-uniformity. 
Figure 5.8 shows the resulting contrast-enhanced flood image for the I converging 
collimator and 511 keV collimator acquired with 140 keV gamma rays to intensify 
collimator artifacts. The images are composed of 11 sub-images, one for each module 
around the circumference of the SPRINT gantry. Summary measures of uniformity for all 
three collimators are provided in Table 5.1. Standard deviation in uniformity was 
calculated from the variance in pixel counts within a region of each module which 
excluded detector artifacts. 

The ^^"’Tc collimator demonstrated the best uniformity, with no gross artifacts 
visible in the image. The I collimator had the greatest variance, with visual artifacts 
including several striations and a large region of slightly increased transmission. Close 
observance of this image also indicates point areas of increased penetration associated 
with the seam connection of each slice. The 511 keV collimator had good measured 
uniformity, but shows obvious striations across its width when contrast is enhanced. The 
fractional non-uniformity corresponds directly to variation in lead foil spacing so that the 
maximum non-uniformity of 14% for the collimator indicates a 0.13 mm (inner 
surface) to 0.20 mm (outer surface) variation in spacing. The non-uniformity in the 
parallel collimators is in good agreement with the variation in bond thickness and the 
manufacturer’s specified 0.1 mm tolerance in foam thickness. 



91 



B. Uniformity of 1-131 converging collimator 



C. Uniformity of F-18 parallel collimator 

Figure 5.13.Uniformity of Collimators: Images are Shown for Each of the 11 
Detectors Comprising the SPRINT Detector Ring. Contrast is enhanced to 40% of 
full intensity range. Structure at detector ends lies outside the useful axial field-of- 

view 
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Table 5.1. Uniformity of Collimators at Different Energies 



Uniformity of collimators, standard deviation offlood images 

Collimator 

yyniTc 

BTi 

18 p 

’^■"Tc Parallel 

5% 



Converging 

14% 

10% 

7% 

511 keV Parallel 

10% 

8% 

6% 


Intrinsic Resolution 

Intrinsic resolution is determined from an iterative unfolding of the detector point 
response function from the joint system point response (collimator and detector)^. 
Detector response is measured using a highly collimated slit source translated along the 
detector surface. The system response is measured using a point source translated along 
the central axis of the gantry. Sample intrinsic collimator response functions acquired for 
the three collimators at adjacent 0.5 mm increments are shown in figure 5.14-5.16. The 
five responses shown cover an interval corresponding to the separation of two adjacent 
lead foil rings. Specifically, slice separation was 2.5 mm for the ^^"'Tc collimator, 1.9 
mm for the converging collimator (outer edge), and 2.4 mm for the 511 keV 
collimator. The measured intrinsic collimator resolution (FWHM) compared to the 
geometric resolution of each collimator is shown in Table 5.2. Range values provided in 
the table represent the mean resolution achieved for measurements taken between septa 
compared to those taken under septa. Full-width at tenth maximum (FWTM) data are 
also provided in the table. System axial resolution is one to two millimeters wider than 
the intrinsic collimator resolution. 
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-Series 1 

-Series2 

.Series3 

-Series4 

-SeriesS 


Figure 5.14. Low Energy Parallel Collimator Intrinsic Response Functions 
Measured at 0.5 mm Intervals Along Center of SPRINT Field-of-View 



-Series 1 

-Series2 

.SeriesS 

-Series4 

-SeriesS 


Figure 5.15. Medium Energy Converging Collimator Intrinsic Response Function 
Measured at 0.5 mm Intervals Along the Center of SPRINT Field-of-View 



Figure 5.16. High Energy Parallel Collimator Intrinsic Response Functions 
Measured at 0.5 mm Intervals Along the Center of SPRINT Field-of-View 
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Geometric collimator resolution generally agreed well with that measured, with 
the ‘^’l collimator and 511 keV collimator having slightly better performance than 
expected. The plotted response functions demonstrate the spatially variant response of 
these collimators, with a minimum response occurring at locations directly behind lead 
septa. This effect is most pronounced for the I collimator which had an almost 
complete loss of response when a septa was aligned directly over the source. In 
comparison, the 511 keV collimator had a reduced spatial dependence because of poorer 
intrinsic resolution as well as septal penetration. 


Table 5.2. Intrinsic Collimator Resolution. Range Values Indicate Average Values 
Between Septa and Average Values Under Septa. 




Axial Resolution 

Collimator 

Geometric 

99niTc 

,3.1 

18p 


Resolution 




99mTc 

8.3 mm 

8.5 ± 2.9 mm FWHM 



Parallel 


(range 8.0 -14.8 mm)* 





13.6±2.0mmFWTM 



Bli 

2.0 mm 

1.1 ±.19 mm FWHM 

1.9 ±.45 mm FWHM 

2.4 ± .38 mm FWHM 

Converging 


1.8±2.7mmFWTM 

3.1±.56mmFWTM 

4.2± 1.2mmFWTM 

511 keV 

3.5 mm 

2.5 ±0.91 mm FWHM 

2.9 ± 1.3 mm FWHM 

3.3 ±1.7 mm FWHM 

Parallel 


(range 2.2 to 5.0 mm) 

(range 2.4 to 7.0 mm) 

(range 2.5 to 7.3 mm) 



4.5 ± 1.6mmFWTM 

4.7± l.SmmFWTM 

5.6 + 2.1 mmFWTM 


System Sensitivity 

Sensitivity is measured for all three collimators using a point source at locations 
showing a moderate, high and low response based on uniformity maps. A 20% energy 
window is used centered on the photopeak. Results are shown in Table 5.3. The ^^""Tc 
collimator when used with a slit aperture set for 2.9 mm FWHM in-plane resolution 
provided a system sensitivity that is comparable to current two-head SPECT systems. 
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using low-energy all-purpose collimators^. ^^'"Tc when used with either the 
converging or 511 keV parallel collimator demonstrated a factor of four reduction in 
sensitivity due to increased collimator length, thicker septa and smaller slice spacing. An 

ini 

additional factor of three reduction is experienced when I is used with the high energy 
collimators due to a reduced photoelectric cross-section (16% detection probability), and 
slightly smaller abundance of the detected gamma ray. The 12 mm thick sodium-iodide 
detectors demonstrate only a 7% photoelectric interaction probability for 511 keV 
photons. However, some sensitivity is gained for positron emitters due to the 200 percent 
abundance of annihilation radiation. Although the two gamma-rays are correlated, the 
poor detection sensitivity makes it unlikely that both gamma rays will be detected. 

Overall there is an approximate factor of 10 reduction in system sensitivity between the 
^^"’Tc collimator and the two high energy collimators. 

Table 5.3. Sensitivity of Collimators as Measured With a Point Source Centered in 
the SPRINT-II Gantry. Values in parentheses represent system sensitivities 


obtained with a 2 mm wide slit aperture. 


Sensitivity in cpm/pCi/deg and cpm/pCi for 2 mm slit 
(2 mm slit results in 2.9 mm in-plane resolution for ^^'"Tc) 

Collimator 

99mTc 



^^■"Tc Parallel 

‘^*1 Converging 

511 keV Parallel 

61+5.2% cpm/|xCi/deg 

(670+5.2% cpm/jiCi) 

17+13% cpm/|iCi/deg 

(180+13% cpm/pCi) 

17+12% cpm/pCi/deg 

(180+12% cpm/pCi) 

6.2+10% cpm/pCi/deg 

(67+10% cpm/pCi) 

5.5+5.7% cpm/pCi/deg 

(59+5.7% cpm/pCi) 

7J±1.S% cpm/pCi/deg 

(83+7.8% cpm/jxCi) 

7.3+4.8% cpm/pCi/deg 

(79+4.8% cpm/pCi) 



Summary 


Novel construction techniques and materials were used to build precision slice 
collimators for SPRINT-II. These methods proved successful over a range of material 
dimensions, and for the production of both parallel and converging geometries. Observed 
collimator performance generally met design objectives. The ^^"’Tc collimator showed 
excellent uniformity and sensitivity, with moderate spatial variance and adequate 
resolution. The and 511 keV collimators demonstrated excellent resolution, but 
suffered from larger variance in uniformity, increased spatial variance and greatly reduced 
sensitivity when compared to the ^^‘"Tc collimator. Measured resolution and sensitivity 
are as expected for the design specifications of each collimator. Variance in uniformity is 
directly dependent on measured variance in foam and bond thickness between slices 
which can be reduced for the parallel collimators by machining the foam rather than 
accepting the manufacture’s tolerances. Although the methods developed are for a 
specific imaging system, they should find general application in the construction of foil 
collimators for other purposes . 
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CHAPTER VI 


CORRECTING FOR APERTURE PENETRATION BY HIGH ENERGY 

GAMMA-RAYS 

SPRINT’S in-plane resolution is determined by a rotating slit aperture, as 
discussed in Chapter III. It is desirably to establish an image resolution which is at least 
equal to or less than that achieved axially. By using a slit opening of 3.125 mm, the 
predicted geometric resolution is approximately 4.4 mm FWHM, comparable to the 
system resolution measured axially (Chapter V). Unfortunately a problem specific to 511 
keV photon imaging is penetration of these gamma rays through the aperture edges which 
results in a subsequent loss of resolution. This chapter will focus on an innovative 
method developed to compensate for aperture penetration in high energy slit aperture or 
pin-hole SPECT. 

Background 

In reconstructing of images acquired with pin-hole and slit aperture geometries 
one usually assumes the measured projection data is exclusively from gamma rays that 
pass through the slit. As gamma ray energy increases, the aperture edges are less 
effective attenuators, resulting in projection data being contaminated by photons that pass 
through the aperture edges as well as through the collimator body. As a result, the system 
point response function is significantly broadened with a subsequent loss in image 
resolution and contrast. Ideally, this problem would be overcome by simply using more 
effective attenuators in the construction of the aperture. Unfortunately, at the photon 
energies of iodine-131 (364 keV, 637 keV) and positron emitters (511 keV), edge 
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penetration is significant even for depleted uranium (DU) or gold. These materials are 
the best attenuating materials that are obtainable. For SPRINT imaging 511 keV photons, 
approximately half the measured projection data stems from photons that penetrate the 
hybrid lead-DU small animal aperture. The effect is graphically depicted in figure 6.1 
which shows a comparison of the aperture penetration predicted by ray-tracing to that 
measured for an '*F line source centered in the SPRINT gantry. For comparison, the 
normalized profile of a measured ^^"’Tc line source is also shown. 

Predicted vs Actual Slit Penetration at 511 keV (2 mm slit) 



0 5 10 15 20 25 30 

Pixel 

Figure 6.1. Demonstrated Aperture Penetration Measured for 511 keV Annihilation 

Radiation 

Due to this effect, high energy planar and SPECT imaging are conventionally 
performed using parallel, converging or fan-beam collimators that are designed with thick 
septa and long lengths to provide effective collimation. For example, the Siemens 511 
keV parallel hole collimator is designed with 0.25 cm septa, 0.34 cm holes, and has a 
length of 7.5 cm*. Such collimators are heavy, suffer from poor intrinsic sensitivity and 
resolution, have severe spatial variance in their response, and parallel collimators do not 
magnify the object as do pin-hole apertures. Thus pin-hole or slit aperture SPECT is an 
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attractive option for high resolution animal imaging if aperture penetration can be 
corrected. 


Methods 


Initial Attempts to Correct for Penetration 


Projection data acquired from *^F line sources indicates that the penetration 
component of the projection appears as a smooth background added to true signal. This 
was observed to be similar in effect to scatter radiation originating from within the 
measured object in low-energy SPECT and PET. Hence, methods successfully used for 
scatter correction were explored as to their efficacy in correcting aperture penetration. A 
direct compensation technique that could be applied to existing projection data is the 
convolution-subtraction technique^ originally applied to SPECT and later to PET image 
reconstruction^. Application of this method begins with the premise that projection data 
is composed of the true signal, Psig(x), and penetration Ppen(x): 


P = P +P 

measured Signal Penetration 


( 6 . 1 ) 


The projected penetration distribution is then assumed to be of the following 
analytical form: 

D 

h.ne = I Pn,.asured * " t ,t)dt (6.2) 

-D 

Where h(x,t) is the convolution kernel describing the penetration as a function of 
X (image plane) and t (detector plane). If an accurate kernel can be determined, the true 
projection signal is estimated by: 

D 

Psignal ~ Pmeasured ~ \Pmeasured*^(^ ~ ^d)dt (6-3) 

-D 

The projection data from line sources translated across the field-of view were 
measured to empirically determine the blurring (convolution) function (figure 6.2). It is 
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observed that the tails of these distributions are highly spatially variant, presenting an 
immediate obstacle to implementing this simple correction method. The tails are foimd 
to be a function of both source location within the field of view, as well as projection 
angle. The resulting three-dimensional kernel is complicated further by being 
discontinuous in x where rays pass directly through the slit opening and the kernel 
becomes zero. Because of these factors, implementation of this technique was not 
pursued. 



Figure 6.2. Spatially Variant Penetration Tails Observed for Line Sources Moved 

Throughout the Field-of View 

These same issues prohibit an iterative solution to this problem, where aperture 
penetration is modeled and incorporated into the system response matrix. The first 
complication is that for each source location and for each aperture position the probability 
of a photon penetrating the aperture and interacting in a given detector element must be 
determined. This would require an extensive Monte Carlo assessment of the SPRINT 
imaging geometry that is both difficult and time consuming. Second, a non-sparse system 
matrix results from the penetration tails extending across the projection space (figure 6.2). 
This infers a finite probability of aperture penetration for all detector pixels from all 
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source pixels thus causing the system matrix to be “full” or complete. For SPRINT, this 
would require a 390 MByte weight matrix, which is too large for RAM storage on most 
computer platforms and substantially lengthens image reconstruction. 

A Proposed Solution: Correction Through Blocking 

The developed method attempts to measure the aperture penetration by “blocking” 
the slit aperture. We observe that the measured projection data is composed of both true 
signal and aperture penetration (Equ. 6.1). If the “penetration” aspect of the projection 
data can somehow be measured directly or estimated from measurements, then the true 

projection signal is estimated as: 

P = P - P 

Signal measured Penetration 

/ ^ ^ (6.4) 

p —i p i p 1— P 

^Signal Signal ' '^Penetration J -^Penetration 

In this case, blocking the aperture provides projection data which serves as an 
estimate of the penetration background. This can then be subtracted from the open slit 
projection to obtain a true net signal. Experiments were conducted on the effect of 
blocking the aperture slit with a depleted uranium bar (3.125 mm x 6.25 mm) and 
measuring the resulting projection for an on-axis line source. The results are shown in 
figure 6.3. The blocked slit data has the expected form of the penetration signal, except 
that it does not drop to zero at the center of the projection. This is because the block is 
not completely opaque to 511 keV gamma rays. The figure also shows the net projection 
data resulting when the blocked slit projection data is subtracted from the open slit data. 
The normalized projection profile of a ^^""Tc source is presented for reference. The net 
projection data is observed to have a resolution that approaches that for the ®^'”Tc data, 
with the penetration tails significantly reduced. The method’s efficacy was assessed for 
extended objects by making a series of off-axis source measurements (figures 6.4 and 
6.5). Again, the method demonstrates good improvement in resolution and tail-reduction 



103 


for the net responses obtained by subtracting the matched open and blocked slit 


measurements. 



Figure 6.3. Projection Profiles for Open Slit and Blocked Slit Apertures 



-Open. Far-Left -Open, Md-Left -Open, Center Open, Md-Rght-Open, Far-Rght 

-Block, Far-Left -Block, Md-Left -Block, Center .Block, Md-Rght-Block, Far-Rght 


Figure 6.4. Comparison of Open and Blocked Slit Profiles for Various Location 

Across the Field of View 
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Figure 6.5. Net Line Source Projections Obtained from Difference Between Open 

and Blocked Slit Data 

Correcting for the Non-Ideal Block Response 

The devised algorithm was applied to imaging a cylinder filled with a uniform 
distribution of activity. Reconstructed images however demonstrated a significant 
depression in counts in the central region of the object. Two factors contribute to this 
effect: the finite attenuating properties of the depleted uranium block and the finite 
thickness of the block. First, 6.2 mm of DU permits approximately 10% of incident 511 
keV photons to be transmitted (based on total interaction cross-sections'^). These photons 
contribute to the measured blocked projection causing an overestimate of the aperture 
penetration where, by definition, penetration should be zero. The second factor causes a 
“penumbral” effect between the aperture edge and block, where photons incident on the 
aperture edge can also be incident on the DU bar. The result in an underestimation of the 
penetration data since the block introduces additional attenuation of the aperture edges. 

















105 


The spatial dependence of these two effects are better observed in Figures 6.6 and 
6.7, where the ray traced profiles through the DU block are presented for four point 
sources in the field of view. For points along the central axis, the block acts as a rect (fl) 
function, attenuating the incident radiation with little penumbral effect. In contrast, far 
off-axis locations approach a A function. An asymmetric A function is reached at a 26.5° 
angle from normal, associated with an off-axis point of 33 mm. The figure also 
demonstrates the severity of the penumbral effect, depicted as the overlap region between 
solid (DU block) and dashed (aperture edge) lines. The effect changes as the source 
location moves towards the edge of the field-of view. 



Figure 6.6. Ray-Tracing from Point Source Locations in Field of View Through DU 

Block and Aperture Edges 
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On-Axis (1, Center) Mid-Left (2,20 mm off-axis) 



0.476 

12.1 





0.758 




Mid-Far Left (3, 30 mm off-axis) Far-Left (4, 40 mm off-axis) 



Figure 6.7, Demonstration of Spatially Variant Penetration and Penumbral Effects 
Observed with Actual Block (Dimensions are in/mm) 


To correct for these effects, an inverted Blackman window was applied as a 
spatial weight to the measured penetration projection data, such that: 


Blackman Weight(i) = \-k* 0.42-0.5• cos 


lK(j - 1) 

(JV-i) 


+ 0.08 • COS] 


47r{i - 1 ) 

(JV-l). 


where i is a projection element, N is the number of elements in a projection (91 or 92), 
and k is an empirically derived constant equal to 0.18. The constant was derived to 
achieve a uniform profile across a uniform object. The resulting function required a 
minimum weighting of 0.82 in the center of the projection, up to 1.0 at the projection 
edges. Results are shown in figure 6.8. Reconstruction of the uniform object using the 
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uncorrected penetration data shows a central depression (figure 6.8A). Filtered back- 
projection reconstruction with the weighted penetration data (figure 6.8B) demonstrates 
recovery of the central area to produce a uniform field. For comparison, the image 
without penetration correction is shown in figure 6.8C, which has approximately twice 
the counts of the corrected image. This is in good agreement with measured projection 
data indicating penetration accounts for approximately 50% of the measured count rate. 




Figure 6.8. Effect of Imperfect Aperture Blocking. Figure A and Bottom Profile 
Indicate Effect of Directly Subtracting Measured Penetration from Open Slit Data. 


Data Acquisition and Image Reconstruction 

To measure aperture penetration as part of a normal tomographic acquisition, the 
slit aperture was modified so that every other slit was blocked by a depleted uranium bar. 
Available material dictated the bar size of 3.125 mm by 6.25 mm, and 150 mm long. 
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This bar size then determined the minimum open slit width of 3.125 millimeters. A 
tomographic acquisition sequence was then implemented which rotated the aperture 80 
steps, 1.5 degree per step, sweeping the aperture in a 120” arc. This doubled the number 
of projections acquired compared to conventional low-energy imaging, providing an 
open slit and blocked slit projection for each aperture position. Post-processing of the 
acquired raw projection data was conducted to form two complete data sets of open and 
blocked slit data. The blocked slit data were then weighted to correct for the effect of 
penetration and penumbra. The modification of the aperture design and the implemented 
sampling protocol are shown in figure 6.9. 



Figure 6.9. Slit Aperture Arrangement and Sampling Protocol for High Energy 
Imaging. Edges of Slit are Edged with DU. Aperture Ring Rotates While Detector 

Ring is Stationary. 


Reconstruction of images was conducted using both filtered back-projection and 
the SAGE iterative algorithm^ In applying FBP, the open-slit projection data was 
preprocessed by subtracting the weighted blocked-slit projection data and then applying 
suitable uniformity correction Pg. 
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Where "w" represents the Blackman window. High frequency noise in the 
reconstructed image was smoothed slightly using a Butterworth filter with a 0.9 cm'* 
frequency cutoff. 

For iterative reconstruction, the SAGE algorithm was applied, with the weighted 
blocked slit projection data used as the random background term: 

Yopen ~ Poisson[P^(Aj.J + )} (6.7) 

The system model employed weights generated using the geometric slit opening, 
noted as Ate, that excluded any aperture penetration effects. Regularization was 
performed to optimize image quality, typically in the range of P=2'^^ for low count 
images to 2'^* for high count images (reference Equ. 4.44). 

Results 

1 S 

Penetration correction was assessed using F line sources arrayed across the field- 
of-view. Reconstruction performed using FBP without penetration correction is shown in 
figure 6.10A, with the associated horizontal profile shown in figure 6.1 OB. Qualitatively, 
the image suffers from poor resolution with long tails introduced on each line as well as a 
tangential and radial dependence in resolution. The same lines reconstructed with 
penetration correction are seen in figure 6.1OC, with a matched profile in 6.1OD. This 
image shows a significant improvement in resolution across the field of view, and 
resolution is now more isotropic. There is also a notable reduction in the tails of each 
source. Quantitatively, the non-penetration corrected image had a central resolution of 
8.8 mm FWHM and 18 mm FWTM due to large penetration tails. Images reconstructed 
with penetration subtraction resulted in a 5.7 mm FWHM and a 10.8 FWTM resolution. 

If the blocked slit projection data are not spatially weighted and used directly for 
correction, a slightly better resolution of 5.0 mm FWHM and 8.6 mm FWTM is achieved 
due to over-correction. Images of lines reconstructed using SAGE (P=2' ) demonstrate 
significant additional improvements in image quality. Qualitatively, backprojection 
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artifacts are removed and resolution approaches the calculated theoretical limit of 4.4 
mm, with a FWTM of 8.0 mm. 

Sensitivity assessed after correction for these line sources was only 66 ± 5.5 
cpm/pCi compared to 127 ± 6.0 cpm/pCi assessed for the high energy collimator 
with a slit width of 3.175 mm (Chapter 5). The approximate halving of sensitivity is 
again attributed to the dual measurement correction method. 



B D 

Figure 6.10. Six Line Source Array Reconstructed Using FBP, Without and With 

Penetration Correction. 

An **F filled micro-Jaszczak phantom consisting of holes ranging from 1.5 to 4.0 
mm in diameter was also imaged. Figure 6.11A shows a ramp filtered FBP 
reconstruction without penetration correction. The image was created with 87.6 million 
counts, but resolution and contrast are poor. Figure 6.1 IB illustrates the same image 
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reconstructed with penetration correction from 40.1 million net counts. Figure 6.11C 
demonstrates the effect of reconstructing using penetration correction and the SAGE 
algorithm (P=2'^^). In uncorrected images, the 4 mm line set was barely discernible and a 
uniform background level was present. Images produced with correction resulted in clear 
observation of the 3 mm holes and negligible background, with the SAGE reconstruction 
demonstrating reduced noise and ray artifacts. 



Figure 6.11. Images of pJaszczak Phantom, Hot Lines from 4.0 mm to 1.5 mm. 


1 s 

Finally, two F line sources (2 mm wide) formed into a “U” and “M” were 
imaged without and with penetration correction, respectively. The resulting FBP images 
are shown in figure 6.12, with the “U” demonstrating poorer resolution and a diffuse 
background as compared with the reconstructed “M”. 



Figure 6.12. U of M, Without and With Penetration Correction. 
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Discussion 

To evaluate in greater detail why a Blackman window was effective in correcting 
the observed depression in central object counts, we conducted a series of Monte Carlo 
simulations for 511 keV gamma-rays imaged with SPRINT®. Simulations were 
performed by Dr. Wilderman of the University of Michigan’s Department of Nuclear 
Engineering. An accurate representation of the SPRINT aperture was simulated using 
three designs: with a 3.125 mm open slit, with the slit blocked by the depleted luanium 
bar, and with the slit blocked by an infinitely thin perfect absorber. Point sources were 
then simulated for each of these geometries both on-axis and at four off-axis locations. 

To simplify implementation, the simulations did not include the detector response. 

Figure 6.13 shows the difference in detected photons falling within a ±10% energy 
window and incident at the detector surface between the DU blocked aperture and 
perfectly blocked aperture for the on-axis and two off-axis points. In these profiles we 
note the two separate effects attributable to the finite size and finite attenuating properties 
of the DU bar. For both the on-axis and off-axis cases, there is notable penetration 
through the bar, compared to zero flux for the perfectly blocked slit. The effect is 
somewhat mitigated in the off-axis cases as the bar is traversed at an angle from normal, 
increasing the effective thickness of the block. As discussed previously, this effect causes 
an overestimation of the measured penetration background that is most pronounced in the 
central region of the projection. The penumbral effect is only demonstrated for the off- 
axis points. This result is evidenced by the diminished intensity in the side lobes of the 
response profile. For the ideally blocked slit, there are highly peaked responses on either 
side of the block, associated with rays passing through region where the aperture knife 
edge meets the block. These responses are greatly diminished as rays incident in this area 
impinge on both the aperture edge and finite size DU bar. 
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Ideal, On-Axis * - - Ideal, Mid-Left -Ideal, Far Left 

DU. On-Axis -DU, Mid-Left -DU. Far-Left 


Figure 6.13. Comparative Point Source Profiles of Aperture Blocked with DU Bar 

and a Perfect Attenuator. 

The combined effects of block penetration and penumbra are better observed in 
figure 6.14 which shows the net difference between the DU blocked and perfectly 
blocked aperture. The net responses are presented for the on-axis source, followed by 
source positions to the far left of the field of view. Values greater than zero are attributed 
to penetration through the aperture block, causing an overestimation in the measured 
penetration data. Negative values correspond to penumbral effects between the aperture 
edge and aperture block and result in an underestimation of the measured penetration. 
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Profiles of DU Blocked - Perfect Attenuator 
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Figure 6.14. Simulated Assessment of Penumbral Effects and Aperture Block 

Penetration. 

The empirically derived Blackman window is explained by examining the ratio of 
integrals of the DU blocked slit data and the perfectly blocked slit data. By plotting these 
normalized ratios as a function of projection centroid, the points in figure 6.15 result. For 
the on-axis case, the ratio is approximately 1.1, in good agreement with the expected 10% 
transmission through the block. As source position is moved to the far left of the field of 
view, the ratio approaches 0.91, indicating the joint effects of penumbra and penetration 
result in off-axis positions being underestimated by 9%. The Monte Carlo data were fit to 
a three-parameter Blackman-Gaussian function which provided a 99.7% correlation 
with the data. This is expression is presented below and as the dashed line in the figure. 
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Where N are the number of projection elements, and ki=0.533, k 2 = 4.54, and 
k3=0.139. 



Blackman Window-Empirical Fit to MC o Monte Carlo 


Figure 6.15. Comparison of Monte Carlo Based Weight to Empirically Derived 

Blackman Window 

Comparing the Monte Carlo results with the plotted Blackman window 
demonstrates a maximum difference of 3.5%. The simulation based function was refined 
by convolving a Gaussian detector response (2 pixel a) to account for the detector 
blurring introduced into the simulation results. This had the effect of further improving 
the agreement between the Blackman to approximately 2.5%. A profile of the uniform 
object reconstructed with the Monte Carlo based response is shown in figure 6.8D. It is 
noted there is only slight difference between the Blackman and simulation based 
corrections. Notably, the simulation based correction is slightly higher peaked centrally 
than the Blackman reconstruction, in agreement with the form of the correction. 
Differences between the two weighting schemes are partially attributed to the limited 
accuracy in reproducing the SPRINT imaging geometry in the simulations, particularly 
magnification effects and blurring introduced by interactions within the scintillator. 
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An Iterative Technique for Improving the Proposed Aperture Penetration 
Correction Method 

The first order correction applied to the penetration data provides acceptable 
qualitative improvements in image quality, but does not actually correct for the non-ideal 
behavior of the block. If the blocked projection data was biased only from penetration 
through the block (i.e. no penumbral effects), further improvements in image quality 
could be achieved through a semi-iterative correction technique. The proposed system 
model for such a solution is; 

Yp ~ Poisson^P^[Ap^A + b)^ (a) 

Yg ~ Poisson^P^{aAp,yX + by^ (b) ^ ^ 

Where Yp is the measured projection data with the open aperture, Ajc are system weights 
derived for the geometric slit opening, and b is the background component composed of 
aperture penetration, scatter and environmental background. The blocked slit projection 
data, Yb, is modeled as the sum of aperture penetration plus penetration through the block 
taken as aAou^, where Adu are system weights generated from Monte Carlo simulations 
of the block penetration, and a is a scale factor. 

The penetration background projection b can be estimated from the following 
expression, which uses the spatially weighted blocked projection data to provide a first 
estimate of the object, 1 fbp- 

OCAdij'^ 

(FBP-SAGE) (6.9) 

This estimate of b can then be used in equation 6.8(a) to solve for the object distribution 
using a discrete estimation algorithm such as SAGE. Iteration between equations 6.8(a) 
and 6.9 provide improved estimates of the object distribution, 1. 

The potential for implementing this method was investigated through a series of 
Monte Carlo simulations of the SPRINT aperture blocked by a cylindrical DU bar. 
Projection data were obtained for the five point source locations used in the previous 
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simulations. Figure 6.16 shows the comparison between the cylinder blocked aperture 
and the perfectly blocked aperture. A comparison of these data with the profiles in figure 
6.13 shows that the cylinder is clearly a less effective attenuator, allowing more primary 
photons to be transmitted. 



Figure 6.16. Comparative Point Source Profiles of Aperture Blocked with 
Cylindrical DU Bar and a Perfect Attenuator. The cylinder is a less effective 
attenuator than the rectangular DU bar (ref fig 6.13). 

The advantage of using a cylindrical block becomes apparent when the net 
penetration through the cylinder is determined from the perfectly blocked aperture. This 
is shown in figure 6.17 as the heavy lined profiles for each of the five source locations. If 
these profiles are compared to those shown in figure 6.14, it is noted that the 
imderestimation of the penetration data is eliminated. The resulting overestimation of the 
penetration data has a form which can be readily incorporated into a sparse systems 
weights matrix, Adu- As an example, point sources at the same positions used in the 
simulations were forward projected using a SPRINT weights matrix generated with a 3 
mm slit. The resulting projections were normalized and overlaid onto the net cylindrical 
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penetration profiles, shown as the lighter profiles in figure 6.17. The good agreement 
between the two sets of data suggests that even the use of a geometric weights matrix, 
Atc, may provide good reconstructions using this iterative block penetration technique 
Optimizing weights matrix generation to match the Monte Carlo based results should 
provide very good reconstructions with improved resolution and sensitivity over the 
demonstrated technique. 



Figure 6.17. Comparison of the Net Penetration Through a DU Cylinder Compared 
to Forward Projected Point Sources Using SPRINT System Weights. Note that 
there is no penumbra with the net responses, and relatively good agreement with 

geometric weights. 

Summary 

Clear improvement in SPECT imaging of high energy gamma-rays can be 
obtained through the proposed penetration correction method. The method allows 
recovery of resolution to within one mm of the system geometric resolution, and 
significEintly improves contrast by removing the long penetration tails of the system point 
response function. Implementation of this method, combined with use of a high energy 
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slice collimator, allows SPRINT to achieve an axial and in-plane resolution of 
approximately 5 mm for 511 keV photons. This should allow good accuracy in 
quantifying organ radioactivity concentrations for large rat organs such as the heart, brain 
and liver. 

There are two important disadvantages imposed by this proposed correction 
method however. First, the protocol requires measuring two complete data sets to arrive 
at an estimate of the true signal. System sensitivity is thus effectively halved. Second, 
the variance in the estimated signal is twice the variance of the penetration background 
since this background is measured twice (reference Eqn. 6.4). If the penetration 
background were modeled and included in the reconstruction process, there would be a 
substantial increase in sensitivity and decrease in signal variance. Unfortunately, this is 
also a complicated approach. The full system matrix produced poses large memory 
requirements and dramatically slows the reconstruction process. Moreover, the limited 
accuracy of modeled penetration will introduce bias in the reconstructed image. 

Further improvements in the “blocked” aperture technique are predicted using 
iterative reconstruction methods that accurately model the imperfect attenuation of the 
depleted uranium block. Such methods should provide an improvement in both 
sensitivity and resolution. 

Although the effects are not as dramatic, pin-hole imaging of I labeled 
radiotracers is also significantly degraded by aperture penetration. Application of this 
method is expected be even more effective for iodine imaging due to better attenuation 
characteristics of the block for 364 keV photons, and the capacity to correct for high 
energy iodine emissions. The devised method should be germane to both planar and 
SPECT imaging and be well suited to clinical application on two-head and four-head 
commercial gamma cameras. 
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CHAPTER VII 


QUANTITATIVE HIGH ENERGY SPECT IMAGING USING SPRINT 

This chapter will address the calibration methods adopted for high energy imaging 
and factors examined to perform accurate image quantification using SPRINT. Methods 
and results are presented for qualitative and quantitative imaging of two-dimensional and 
three-dimensional phantoms and live rats. Finally, a comparison is made between 
SPECT images acquired using SPRINT, and those acquired under similar conditions 
using positron emission tomography. 

SPRINT Calibration and Factors Addressed for Image Quantification 

System Set-Up and Calibration 

SPRINT was originally designed for low to medium energy imaging up to 364 
keV. Calibration of SPRINT for high energy imaging required adjusting system gain so 
that the photopeak from 511 keV photons was at an energy (voltage) approximately equal 
to 140 keV photons. This was required because of the limited dynamic range of the 
position calculation circuits. To accomplish this, the photomultiplier tube high voltage 
was lowered to reduce the tube gain of all detector modules. Module specific energy gain 
factors and lower level discriminators were then adjusted to fine hme photopeak and 
spectrum agreement between energies. This relatively simple set of adjustments greatly 
simplified the overall calibration process and permitted use of existing system calibration 
software. For 511 keV imaging, a phototube voltage of 855 volts gave an output pulse 
height equal to the 140 keV photopeak acquired with a phototube voltage of 1000 volts. 
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Calibration of SPRINT after these adjustments followed procedures identical to low- 
energy imaging. This included adjusting each phototube gain within a module to a 
module specific target gain, generating energy correction maps for a module', and 
establishing energy windows that trigger event detection and position estimation. A 
detailed description of calibration procedures is available elsewhere^ Energy spectra 
observed for 511 keV photons were found to have photopeak centroids in good agreement 
with ^^"'Tc, and energy resolution of 13% FWHM compared to 12% for 140 keV photons. 

Event Position Estimation and Edge Effects at High Energy 

Gantry modules use maximum likelihood position estimation to establish the XY 
position coordinates of photons interacting within the scintillator. Weights used for 
position estimation are derived using a highly collimated source of 140 keV photons that 
is translated parallel to the module surface. Photomultiplier tube responses are recorded 
for each source position, forming a matrix of response versus position. Maximum 
liklihood position weights are computed as a function of this response. These weights are 
used for all photon energies imaged by SPRINT. Application of low-energy position 
weights to high energy imaging is expected to introduce errors in position estimation due 
to the differences in energy deposition within the scintillator. 

To determine the extent of position estimation error and assess the detector 
response for high energies, Monte Carlo simulations were used to compute the spatial 
energy distribution and centroid position as a function of incident photon energy, 
scintillator thickness, and scintillator block dimensions . Table 7.1 summarizes the 
salient results of these simulations for a highly collimated beam of photons incident on a 
5 cm block of Nal(Tl). Figure 7.1 shows the vertical energy deposition for events which 

' Energy correction maps are generated for each module to correct for small gain 
differences among phototubes, and differences in light collection at different locations 
within the detector. 
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result in a total energy deposition greater than 460 keV (10% energy window) for 12.7 
mm and 25.4 mm thick crystals. We observe from the figure that energy deposition for 
364 and 511 keV photons is roughly uniform throughout the crystal. The vertical 
centroid is in fact greater than half the scintillator thickness for both energies. Since 
photons can be incident on the detector modules at angles up to 22 degrees, a uniform 
energy deposition in the detector can introduce a maximum 3 mm blur parallel to the 
detector plane. 


Table 7.1. Simulation Results of Highly Collimated Photons Incident on a 5x5 cm 
Nal Slab Sampled at 0.5 mm Intervals with a Lower Energy Window of 460 keV. 


Incident 

Energy 

(keV) 

Total 

Slab Fraction 

Thickness of Energy 
(mm) Deposited 

Fraction 
On-Axis 
of Ray 

Fraction 
Off-Axis 
of Ray 

Vertical 

Centroid 

(mm) 

Horizontal 

Centroid 

(mm) 

Error* 

±(mm) 

141 

12.7 

0.93 

0.85 

0.15 

3.71 

0.31 

±0.0029 

364 

12.7 

0.28 

0.53 

0.47 

6.06 

1.75 

±0.0213 

511 

12.7 

0.17 

0.42 

0.58 

6.27 

2.29 

±0.0292 

141 

25.4 

0.98 

0.84 

0.16 

4.34 

0.35 

±0.0035 

364 

25.4 

0.51 

0.46 

0.54 

11.26 

2.50 

±0.0188 

511 

25.4 

0.35 

0.33 

0.67 

11.81 

3.08 

±0.0296 


* The standard deviation of the average radial centroid position. 


Additional blur results from scatter reactions that displace the centroid of incident 
photons radially from the line of incidence. Table 7.1 indicates that for 511 keV photons, 
approximately 60% of the energy is deposited in pixels off-axis of the incident beam, as 
compared to 15% at 140 keV. Simulation data indicates that the energy is deposited in 
an exponential profile around the central pixel (figure 7.2). The detector response 
observed in this figure has an approximate resolution of 3.5 mm off-axis for the 12.7 mm 
scintillator and 4.3 mm FWHM for the 25.4 mm scintillator. These tails cause additional 
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resolution loss as compared to 140 keV photons, which have much smaller centroids and 
resolution. 

Errors introduced in position estimation at high energy are indicated by figure 7.3 
which shows a comparison of SPRINT flood image profiles acquired using ^^‘"Tc and 
with a 20% energy window (±10%). The normalized fluorine data demonstrates a jagged 
pattern across a given module caused by an enhanced response at locations between tube 
locations. Position estimation is worst at the module edges where the light response 
function is distorted. Figure 7.3 shows a depressed response for 511 keV photons at the 
edge of most modules. This is better observed in figure 7.4 which shows the contrast 
enhanced difference image between a Tc flood and a normalized F flood. White 
areas in the image represent areas of diminished response for 511 keV photons. The left 
side of module ten demonstrates a region of gross signal deficit thought to be partly due 
to decoupling of the photomultiplier tubes from the scintillator. 
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Figure 7.1. Vertical Energy Deposition as a Function of Incident Energy. 
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Radial (Off-axis) Energy Distribution 



Figure 7.2. Horizontal Fractional Energy Deposition as a Function of Incident 
Energy. Approximate Resolution of Response is Presented, Relative to a Maximum 
Fraction of Off-Axis Energy Deposited of 0.06. 


Profile through ®®'"Tc and ’®F Flood Images 



0 100 200 300 400 500 


-Technetium-99m profile -Fliiorine-IS profile 


Figure 7.3. Comparison ofand ^*F Flood Profdes. Note the Regular Pattern 
in the 511 keV Profile and the More Severe Variation at the Module Edges Prior to 

Diminishing to Zero. 










126 



Figure 7.4. Difference Flood Image for and ‘*F. Bright Bands at Module 
Edges Indicate Low Count Areas For **F 


There are limited options available to overcome the blurring and distortion effects 
of high energy imaging. Ideally, detector position weights would be generated and 
applied for each imaged energy. However, position weight generation is time consuming, 
requires a highly collimated source and is performed for each module in a calibration set¬ 
up prior to installation in SPRINT. The difficulty in creating a highly collimated point 
source of annihilation radiation, limited flexibility in existing SPRINT software, and the 
complexity of the calibration process mitigates against a separate determination of 
position weights. Instead ^^"’Tc derived weights are applied, with some degradation in 
image quality. Edge affects on the other hand were found to have a severe effect on 
image quality since back-projection of erroneous low count data greatly increases 
structural noise. To minimize this effect, projection data was modified to exclude low 
count detector bars. While ^^""Tc imaging makes use of all 488 available bars, *^F 
imaging uses only 459 bars, reducing the geometric sensitivity by 6%. To partially 
compensate for the uniform energy deposition through the scintillator depth, the system 
weights matrix used in SAGE reconstruction incorporates a 6 mm mean interaction depth. 

Decay and Deadtime Correction 

Positron emitters have a relatively short physical half-life, on the order of 2 
minutes for to 110 minutes for **F. Further, the high photon energy of annihilation 
radiation results in poor photopeak sensitivity for conventional SPECT systems. These 
factors, together with the relatively low sensitivity of collimated systems, dictate that 
large amounts of activity be used with long imaging times to obtain images of acceptable 
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signal to noise. The use of long count times for SPRINT requires application of decay 
correction to the projection data, since significant decay of the radiotracer can occur 
during imaging. The use of high activities introduces the requirement for deadtime 
correction, since large activities result in large module count rates. 

The effect of deadtime in SPRINT is observed in figure 7.5, which presents the 
number of X-Y position board starts as a function of ^^“Tc source activity. The count 
rate reflects scintillator events above the lower level discriminator that triggers position 
estimation. The graph indicates module three has a significantly worse deadtime than 
other modules, and module six has a somewhat smaller deadtime than other modules. 

X-Y Position Board Starts Uncollimated Source Strength 

Module 



Figure 7.5. Observed Module Count Rate as a Function of Source Activity. 

Deadtime in the modules is principally nonparalyzable and assessed using the 
two-source method as discussed by Knoll^. Two sources of approximately equal strength 
were measured for equal times in the uncollimated SPRINT gantry. Sources were 
measured together and separately, and a background measurement was acquired. 
Acquired images were analyzed for total counts occurring in each module, and then 
divided by the acquisition time to provide a count rate. Deadtime was then calculated 
from the rate data. 










128 


Deadtimes are presented in Table 7.2, calculated from X-Y position starts and 
total image counts. The difference between image based deadtime and start pulse based 
deadtime is unexplained and being investigated. Deadtime correction is incorporated into 
the projection shuffling algorithm by scaling the projection data by: 

Dead Time Correction Factor = —— 

\ + nT 

where n is the total measured count rate in the image, including additional 
window images if acquired, and t is the module deadtime. For extended sources such as 
the micro-Jaszczak phantom where 40 mCi of activity were imaged, deadtime corrections 
were found to be as high as a factor of 1.9. Low activity line sources had deadtime 
corrections equivalent to one. The effect of deadtime correction on an image is 
demonstrated in figure 7.6 showing a uniform cylinder phantom reconstructed without 
and with deadtime correction. The associated profiles, though noisy, demonstrate a much 
more uniform image with correction. 

Table 7.2. Deadtimes for SPRINT Modules. 


Deadtime (s) 

Modul 

e 

X-Y Start 

Image (^^'"Tc) 

Image (**F) 

1 

2.42E-05 

4.07E-05 

6.02E-05 

2 

2.41E-05 

4.05E-05 

5.45E-05 

3 

4.34E-05 

7.67E-05 

0.000104 

4 

2.42E-05 

4.26E-05 

6.58E-05 

5 

2.45E-05 

4.43E-05 

7.95E-05 

6 

2.42E-05 

4.16E-05 

6.92E-05 

7 

2.34E-05 

4.41E-05 

7.81E-05 

8 

2.33E-05 

4.78E-05 

7.6E-05 

9 

2.38E-05 

4.25E-05 

6.57E-05 

10 

2.35E-05 

4.63E-05 

6.11E-05 

11 

2.37E-05 

4.47E-05 

5.93E-05 


The significance of decay correction is clearly dependent on the half-life of the 
nuclide and the count time. For 80 step acquisitions performed over a 55 minute period 
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(40 seconds per step), there is a 1.4 correction factor for the last projection compared to 
the first for and a 6.7 correction factor for *'C. Since one hour long acquisition times 
are typically employed with SPRINT, such large corrections make imaging of carbon-11 
radiotracers prohibitive. Decay correction is applied during the projection shuffling 
operation described in Chapter VI, directly after deadtime correction. A given projection 
is corrected to the start of the acquisition using the elapsed time between the start of the 
tomographic acquisition and the time of the projection measurement. The effect of decay 
during a given projection measurement is trivial for and longer lived radionuclides. 




-Uncorrected Profile-Profile w ith Deadtime correction 


c. 


Figure 7.6. Comparison Of Uniform Cylinder Imaged Without (A) and With (B) 
Deadtime Correction. Profile Through Center of Object (C) Shows Reduced Noise 

and More Uniform Appearance. 
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Flood Correction 

To correct for non-uniformity in the sliee collimator and the non-uniform response 
of the actual SPRINT detector modules, extrinsic flood correction maps are generated. 

The maps provide pixel-by-pixel relative sensitivity correction factors for each module, as 
well as normalization for the variation in sensitivity between modules. They are 
determined from long acquisitions of a uniform line souree with the slice collimator 
installed in the SPRINT gantry. Because of the slice collimator, this required use of a 
moderate to high activity line source of 500 pCi (18.5 MBq) for ^^""Tc and 700 pCi (30 
MBq) for imaged for 50 minutes. Resulting flood images for each of the eleven 
modules are shown in figure 7.7. 



A 99nij^ Flood Image acquired with Low-Energy Collimator 



18 

B. F Flood Image Acquired with High Energy Collimator 

Figure 7.7. Extrinsic Flood Images of’’""Tc and ^^F. Images are Generally Similar 

In Appearance. 


The flood images demonstrate good uniformity for both nuclides. Characteristic 
artifacts corrected by the maps include the visible tube pattern in several of the modules 
as well more subtle azimuthal bands. In the '*F flood, point defects can be observed at 
inter-tube locations associated with point errors in the module energy eorrection maps. 
Gross defects in modules 10 and 11 can also be observed. Module 10 is thought to have a 
defect in the scintillator or has experienced decoupling of the phototube from the light 
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pipe. Module 11 has an electronic defect in the X-position electronics. These artifacts 
are partially handled for ^^"’Tc, but are too severe for correction with *^F. 

The long flood images are converted to flood correction maps by determining the 
module deadtime correction factor and then averaging the central 400 pixels of each 
detector. Average values are divided by all deadtime corrected pixels to derive a 
correction factor up to a user specified maximum value. Generated maps are then applied 
to the raw projection data during reconstruction. 

Attenuation and Scatter Correction 

Attenuation correction and scatter correction must be considered for accurate 
quantification of organ concentrations'^. For SPRINT imaging 511 keV photons, the 
significance of attenuation and scatter within an extended object is shown in Table 7.3. 
For animal bodies and phantoms, the maximum diameter object imaged is limited by the 
field-of-view to be ~8 cm. For photons emitted at the center of a tissue equivalent 
cylinder of this diameter, there is an interaction probability of 25% to 32% based on a 
total attenuation cross-section. Since attenuation and scatter directly affect the measured 
line integral through the object, a correction for each is desirable in the reconstruction 
algorithm to provide images that accurately reflect activity concentration. 

Attenuation correction is conventionally performed in filtered backprojection 
using a first order or second order Chang correction^ An attenuation coefficient for this 
method is derived from a theoretical value for the attenuating medium and the energy of 
the gamma ray. If scatter correction is first applied, the total attenuation coefficient of 
0.095 cm'* for 511 keV photons in water is used^. Effectively, a first order Chang 
correction method derives a pixel by pixel correction factor for attenuation at point (xo,yo) 
equal to: 
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4^0 ,yo) = - 

where M is the total number of projections taken in a 360° scan measurement and le is the 
distance along the ray between point (xo,yo) and the boundary point of the object at 
projection angle 9. Alternatively, iterative reconstruction allows direct integration of 
object attenuation into the system weights matrix. In this case, element-by-element 
multiplication is performed using a two-dimensional attenuation map of the object (which 
can be other than uniform) and the system weights matrix^. 

Table 7.3. Effective Transmission and Interaction Probabilities for Annihilation 
Radiation in Various Tissue Thickness.’* 


Thickness 

(cm) 

Non-interacting 

Transmission 

Interacting 

Probability 

Scatter/Primary 

Ratio (SF)' 

1 

91% 

9% 

0.6% 

2 

82% 

17% 

1.4% 

3 

75% 

25% 

2.1% 

4 

68% 

32% 

3.2% 

5 

62% 

38% 

3.9% 

6 

56% 

44% 

4.4% 

7 

51% 

49% 

4.9% 

8 

46% 

54% 

5.6% 


'The ratio of 511 photons emerging from the phantom which have been scattered to those 
which escape without interaction. 

Scatter is by comparison a more complicated process to correct. Originally, we 

had planned to use the dual energy window method, originally proposed by Jaszczak ^ 

In essence, the scatter corrected projection data is given by: 

P = p -kP 

^ con ^ pp scat 

where Pcorr is the corrected projection, Ppp is the photopeak measured projection and Pscat 
is a projection measured from a lower energy window in the Compton scatter region of 
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the spectrum. The factor k is empirically derived using the method described by Koral et 
al}\ where k is assumed to be independent of position in the object. Unfortimately this 
method could not be adapted to high energy imaging. Spectra acquired of an line 
source in air and within a 7.3 cm phantom filled with water are shown in figure 7.8. 

There is a observed a very high Compton background both with and without the presence 
of a scattering medium. By comparison, the baseline scatter level is 83% of that when the 
phantom is used. Much of this high Compton background is attributed to Compton 
interactions within the scintillator followed by an escape of the scattered photon. 
Approximately 83% of the interactions in the scintillator with 511 keV photons are from 
incoherent and coherent scatter". Scatter is also contributed from interactions with the 
slice collimator and slit aperture. 


Comparison of Energy Spectra for Line Sources in Air and in 
Center of Water Filled Phantom 



Figure 7.8. Comparison of Energy Spectra from Normalized Line Sources in Air 
and in Water Filled Phantom. Note Comparable Spectra and High Compton 

Regions for Both Cases. 


" Based on NNDC cross-section values for Nal(Tl) and a 12.5 mm scintillator thickness. 
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Consequently, use of a secondary energy window does not accurately reflect 
Compton interactions within the object. As an alternative, a first-order attenuation 
correction was investigated to correct for both attenuation and scatter using an effective 
attenuation coefficient that is less than the theoretical value A value of 0.06 cm'* 
was empirically derived from projection data of a line source in the center of a water 
filled phantom. When this value was applied to FBP reconstructions of the line source in 
the phantom, total counts from the line image were found to agree well with that for an 
uncorrected line source in air. Comparative profiles of the two reconstructed images are 
shown in figure 7.9. Regions of interest constructed around the lines indicate that the line 
in air was detected with an average sensitivity of 110 ± 2.5 cpm/pCi, whereas the 
attenuation corrected line source within the phantom was detected with a sensitivity of 


108 ± 2.39 cpm/pCi. 

Horizontal and Vertical Profiles through Line Source 
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Figure 7.9. Line Source Profiles, in Air and Attenuation Corrected Within 
Phantom. Profiles Have Been Normalized by Only Decay Correction. 
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The use of an effective attenuation coefficient does not truly correct for scatter. 

To determine the errors in quantification which could be expected fi-om scattered 
radiation, Monte Carlo simulations were performed. The geometry assessed was a 
cylindrical tissue equivalent phantom filled with uniformly distributed activity. Results 
using the SIMIND code*"^ are presented in Table 7.3. Note that approximately 5% of 
photons emanating from a 7 to 8 cm cylinder have undergone a scattering reaction. Since 
laboratory rats imaged with SPRINT have a diameter usually smaller than 8 cm, it is 
reasonable to expect that scatter will introduce small errors in image quantification. 

Correction for Background Radiation 

Background radiation adds counts to acquired projection data and can affect 
contrast in the reconstructed image if background count rates within the photopeak 
window are significant. Background in SPRINT stems from environmental background 
and that from the depleted uranium bars mounted in the slit aperture. Flood images 
acquired without any external sources in the gantry indicated that background introduces 
a uniform exposure of the modules with no significant spatial pattern observed. Figure 
7.10 shows the measured energy spectra from environmental background with and 
without the aperture present. During SPECT imaging, typical activities used generate 
photopeak coimt rates ranging from several hundred to well over ten thousand counts per 
second. It is apparent fi-om the background count rates shown in the figure that 
background radiation is typically less than 1% of the measured projection data, and is 
consequently ignored. 



136 


Energy Spectra of Background w/ and w/o Aperture 



-Open Slit -Open Slit (Photopeak) 

-Covered Slit -Covered Slit (Photopeak) 

.Aperture Removed - - - -Aperture Removed (Photopeak) 


Figure 7.10. Energy Spectra of Background Radiation. 


Summary 

Raw projection data are manipulated through the application of numerous 
corrections to enhance both image quality and quantitative accuracy. These factors are 
summarized in Table 7.4. The most substantial correction was that devised for aperture 
penetration from high energy gamma rays. This was found to have a dramatic impact on 
both image contrast and resolution. Other significant corrections included the exclusion 
of low-count module bars, deadtime correction and decay correction. These corrections 
improved the noise characteristics of the reconstructed image as well as accuracy. 
Similarly, flood and energy correction maps serve to provide modest improvements in 
image quality and uniformity. Attempts to correct for the biasing effects of scatter, 
attenuation and background were also examined. Because of the difficulties encountered 










137 


in applying a dual window scatter correction technique, a simple Chang attenuation 
correction was applied using an effective attenuation coefficient for FBP reconstructed 
images. This effective attenuation coefficient was also integrated into the system weights 
matrix (assuming a uniform object) for SAGE iterative reconstructions. Since scatter 
effects were predicted to be small, this correction was thought adequate for the purposes 
of this study. 


Table 7.4: Corrections Applied to SPRINT Projection Data and Incorporated into 

Image Reconstruction. 


Correction Factor 

Primary Issue and General Effect and Limitations 

Aperture Penetration 

Corrects for high energy photons penetrating the aperture. 

Results in improved resolution and contrast. 

Module Edge Exclusion 

Removes low count detector data to improve image contrast 
and remove structured noise. 

Deadtime 

Corrects for loss of counts due to finite pulse processing time. 
Particularly important at high count rates. Improves image 
uniformity and quantitative accuracy. 

Decay 

Corrects for loss of counts as radionuclide decays during 
imaging. Improves quantification accuracy. 

Energy Mapping 

Corrects for small gain differences among phototubes, and 
differences in light collection at different locations within the 
detector. Improves system uniformity and energy resolution 

Flood Mapping 

Corrects for pixel to pixel and module to module variation in 
detector and electronic sensitivity. Improves system 
uniformity. 

Attenuation/S carter 

Corrects for attenuation and scatter of primary photons emitted 
within the object. Uses an effective attenuation coefficient for 
uniform object. The method is simple and not does not really 
“correcf ’ for scatter. Improves quantification accuracy. 
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Quantitative Imaging of 2-D and 3-D Phantoms 

SPRINT quantitative performance was assessed initially using a two-dimensional 
and a three-dimensional phantom. This provided known object activities that could be 
directly compared to results obtained by region-of-interest analysis of the reconstructed 
SPECT images. 

The two-dimensional phantom and three-dimensional phantom used are shown in 
figures 7.11 and 7.12, respectively. The phantom body consists of a 6.8 cm inner 
diameter, 5.7 long cylindrical vessel (left) with a screw-on o-ring sealed lid (right). The 
empty phantom has a volume of 236 ml. On the surface of the phantom holder are two 
sets of ten cylindrical holes each with a 10 pi volume. One set of holes fo rms a straight 
line along the phantom body at 6 mm intervals. The second series follows a helical 
pattern around the body, with hole positions in the same axial plane as the straight line 
set. Holes were filled with a Gd-DTPA solution during MR imaging, and an '^F solution 
during SPECT imaging to serve as fiducial markers for image comparison. The base and 
lid of the phantom holder are fitted with injection ports to allow for controlled 
introduction of radioactive material and complete filling of the phantom. 

The two-dimensional phantom insert is seen in the center of figure 7.11. It 
consists of a base fitted with two holes to accept the hub of two 10 ml plastic syringes 
trimmed to a filled volume of 6 ml. It is attached to the phantom body with nylon screws. 
Ten milliliter syringe bodies were chosen because of their ease of filling with a precise 
volume, and because their inner diameter (14 mm) was greater than two times the system 
resolution (~5 mm). This was required for accurate quantification. The phantom design 
allowed for one rod to be filled with a “hot” concentration of radioactivity approximately 
five times the concentration of the “warm” background solution in the body of the 
phantom. The other rod does not contain activity and is used to assess scatter. 


Figure 7.11. Two-Dimensional Syringe Phantom, with Holder (Note Fidueiary 
Markers) and Lid (Note Injection Port). 

The three-dimensional phantom insert is observed in the center of figure 7.12. A 
base was constructed to accept two plungers fi-om of a 5 ml syringe. These served as 
support posts for two hollow plastic “eggs” each with a 3 ml volume, and linear 
dimensions of 15 mm by 24 mm. These objects were again chosen because of their 
relatively precise volume, ease of filling, and acceptable linear dimensions. Eggs were 
sealed after filling and affixed to the support stems with small dabs of hot melt glue. As 
with the 2-D phantom, one egg is filled with a concentration five times the background 
activity concentration, and one egg is “cold”. 








Figure 7.12. Three-Dimensional Phantom Consisting of Two 3 ml Egg Shaped 
Objects. One was Filled with Five Times the Surrounding Activity Concentration, 

One was Cold. 

SPECT Imaging Procedures 

Fluorine- 18 in water with a minimum activity of 40 mCi and a minimum 
concentration of 6 mCi/ml was acquired from the University of Michigan PET Facility. 
Phantom objects were filled and sealed such that the “hot” object had an activity 
concentration five times that of the surrounding background activity. Object activities 
were assessed using a Capintec CRC-4R dose calibrator. The other object in each 
phantom was “cold” and was filled with either epoxy resin (rod) or distilled water (egg). 
Fiducial markers were filled with residual '^F solution to provide an activity of 50 to 100 
pCi per marker. Fiducial marker holes were covered with electrical tape to prevent 
leaking. After assembling the phantom, the body was filled with distilled water and 
activity was added using the injection ports. The sealed phantom was then shaken, not 
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stirred. Total initial activity within the egg phantom was 30 mCi (1110 MBq) as 
compared to 25 mCi (925 MBq) for the rod phantom. 

Phantoms were imaged in SPRINT directly after the phantom was assembled. 
SPRINT acquisition parameters for the 2-D phantom (one acquisition) and the 3-D 
phantom (8 acquisitions) are presented in Table 7.5. The phantoms were centered 
radially and axially within the gantry, and planar imaging was performed to verify 
position. Two energy windows were set prior to acquiring projection data: a photopeak 
window set to ±10% the peak channel; and a scatter vvdndow set three channels below the 
photopeak window and equal to the photopeak window width. The phantoms were then 
imaged using SPRINT tomographic acquisition software. For the rod phantom a single 
55 minute (27 minute effective count time) acquisition was performed. The egg phantom 
was imaged sequentially over a 10 hour period using two acquisition times, 55 minutes 
and 10 minutes (4 minute effective count time). These were the maximum and minimmn 
acquisition times permitted by existing software when using aperture penetration 
correction. The decay of the phantom permitted observation of image degradation as total 
image counts decreased by a factor of 18 over this time. The short acquisition times gave 
additional data with a factor of 10 reduction in image coimts. 

It was found that SPRINT phototubes drift substantially over time. To partly 
correct for this effect, energy correction maps were generated within 12 hours of imaging, 
and photopeak windows were set just prior to imaging. Flood correction maps were then 
generated from flood images acquired immediately following (or just prior to) phantom 
imaging using a uniform line source with the slice collimator in-place. 

Image reconstruction was performed using filtered backprojection and the SAGE 
iterative algorithm. Raw acquisition data was first reorganized to create open and 
weighted blocked slit projection data sets. Decay correction and deadtime correction 
were applied at this time. This data was processed to form flood corrected slice 
projection sets (sprintograms) which had low count bins rejected. Finally reconstruction 
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was performed using the two algorithms. In filtered backprojection, penetration 
correction was performed, followed by filtering, backprojection and attenuation 
correction. In SAGE, attenuation was addressed within the system weights. A 
regularization constant (P) was selected to provide a desired level of smoothing, and the 
system model directly incorporated the weighted penetration data as the random 
background. 

Accurate quantification of the resulting images required calibration of SPRINT 
with a known activity source. Two different calibration factors were used: one based on 
the known concentration in the 3-D phantom’s uniform background; the other based on 
total activity in a capillary tube filled with a known lineal activity. Image reconstruction 
was performed again using FBP and SAGE, with attenuation for the phantom, and 
without attenuation correction for the line. Acquisition and reconstruction details are 
provided in Table 7.5. Reconstructed line source images were analyzed by scribing large 
regions of interest (ROIs) around the line to measure the total counts in the line image. 
The counts were divided by the total acquisition live time to arrive at a count rate. The 
calibration factor was then calculated using the decay corrected lineal activity of the 
source and the slice thickness of the reconstructed image: 


Calibration Factor (cpm/ uCi) = 


ROI Countrate (cpm) 


Lineal Activity (pCi/ mm)S\\CQ Thickness (mm) 
Images of the 3-D phantom were analyzed by scribing large regions of interest 
within the uniform background region of several 6 mm thick slices. This provided the 
average number of counts per voxel, and from the total acquire time, the average count 
rate per voxel. Activity concentration of the phantom was determined from the known 
injected activity and the phantom volume. The calibration factor was then determined as: 


Calibration Factor(cpm/pixel/ uCi / ml) = 


ROI Countrate^'c;?^/ pixel) 
Activity Concentration^" uCi/ml) 
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Table 7.5. Acquisition Parameters for Rod (2-D) and Egg (3-D) SPRINT Imaging. 



ROD Imaging 

Egg Imaging 

Line Source 

No. Steps 

80 

80 

80 

Time/Step 

40 s 

40 or 6 s 

30 s 

Total Time 

55 m 

55 or 10 m 

48 m 

Live Time 

26.67 m 

26.67 or 4 m 

20 m 

Reconstruction: 

SAGE/FBP 

SAGE/FBP 

SAGE/FBP 

FOV 

86x86 mm 

86x86 mm 

86x86 mm 

Image Size 

100x100 

100x100 

100x100 

Slice Thick: 

6 mm 

3 mm 

6 mm 

No. Slices 

1 

10 

12 

Filters (FBP) 

Ramp/Butterworth 

Ramp/Butterworth 

Ramp/Butterworth 


0.9 cm'' cutoff 

0.9 cm'' cutoff 

0.9 cm'* cutoff 

Beta/Convergence 

Lvl (to halt iterations) 

2'^’/0.01 

2'^VO.Ol 

2'^VO.OI 


Image analysis of the phantoms was conducted by scribing small regions of 


interest well within the boundaries of the phantom objects. This was necessary to avoid 
the affects of finite spatial resolution in the reconstructed images. Typical ROIs for the 
rods and eggs contained less than 20 pixels. For the hot object in each phantom, the 
maximum single pixel counts was also recorded. An average pixel count rate (cpm/pixel) 
was calculated from the average counts per pixel within a region and the total acquisition 
time. The activity concentration of a given region was then calculated using both of the 
following equations with the appropriate calibration factor: 

n j. ^ ROI Specific Countratefcpw/mxe/J 

RadioactivityConcentration( uCi / ml) = - - - — - - -^- 

Cal Factor^ cpm/ uCi) ■ Voxel Size^' ml/pixel) 

„ j. ^ , 1 , ROI Specific Countratefcpw/nixc/) 

RadioactmtyConcentration( uCi / ml) = - - - — - - -^ 

Cal Factorf^ cpm/pixel / uCi /ml)- 
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The error as a percent difference from each single measurement was determined 
from the above calculated concentration and the known concentration of the phantom 
regions. 

To provide a relative indication of noise in the uniform regions of an image, the 
standard deviation in the average pixel counts was calculated. It must be emphasized that 
this standard deviation, and the resulting deviation in that region’s concentration, should 
not to be confused with the actual standard deviation for this region resulting from the 
measurement process. The mathematical operations of backprojection, ramp filtering and 
other filtering (smoothing) operations introduce correlation between the pixels of a 
reconstructed image. The nature of statistical iterative reconstruction introduces similar 
covariance between pixels, dependent on the degree of regularization. Consequently, it is 
not possible to differentiate the bias nor variance in the measurement process from a 
single measure. The measured standard deviation in a given single region is used as only 
an estimate of “noise” in the region, and the percent difference calculated for a region is 
taken as an approximate “error” in the assessed concentration. 

Actual bias and standard deviation were estimated from the serial acquisitions of 
the 3-D phantom. Since there is considerable decay occurring between each acquisition, 
it was not possible to treat the estimated concentrations from each acquisition as identical 
independently distributed random processes. However, without knowing the distribution 
of estimates from each acquisition, it is possible to derive weighted least square 
estimators for the mean and the variance of these observations*: 

t^Fri 

- 

DF, 

(=1 

Where DFi is the decay correction factor for acquisition i, ^ o is the estimated 
activity concentration at t=0 for measurement i, and N is the number of observations. The 
decay corrected activity concentration for acquisition i is given by: 
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ROI Countrate; f 

=- - where DF, = exp\ 

Cal Factor • DF, V 


■ln(2)ATi 


M/2 / 


Substituting the expression for qo into the equations for the estimated population 


mean and variance results in: 

^ ROI Countratej 
^ CalFactor 


N 

t.DF, 

/=1 


and 


r 


qo 


N 


1 ^ DF^ 

-IttDF, ' 


ROI Countrate 


^ ROI Countratej 
^ CalFactor 


i ;=I 


CalFactor • DF, 


I,DF, 


i=\ 


The subsequent bias and STD are estimated as: 

A 

bias^a^„-q^„^„ 

STD^^, 

It is recognized that due to deadtime corrections, the process is not truly Poisson, 
and hence the above expressions are only approximate. 


MRI Imaging Procedure 


To assist in image analysis, phantoms were imaged using nuclear magnetic 
resonance to provide a high resolution physical template of the phantom cross-section. 
The instrument used was an Omega CSI2T Animal Imaging System manufactured by 
General Electric. Phantom imaging in this system required construction of a small bore 
receiver probe to improve the object signal-to-noise ratio. This probe was designed with 
an 8.25 cm irmer diameter, which is 5 mm larger than the outer diameter of the phantom. 
This provides an improvement in signal by a factor of 3.5 as compared to using the 
system’s large probe. 

Image acquisition using the CSI 2T for the egg phantom and rod phantom was 
performed a minimum of 24 hours after SPECT imaging to allow complete decay of the 
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original ’^F. Fiducial markers were emptied and replaced with a 0.1 mM solution of Gd- 
DTPA which provides a bright signal during MR imaging"'. 

The prepared phantom was placed in the small bore receiver probe, which was 
then installed in the CSI2T core. The phantom was oriented so that the linear markers 
were superior. The probe was tuned to establish the proper excitation frequency, and a 
single coronal image was acquired through the fiducial markers to determine slice 
locations. System parameters were then adjusted to acquire eight 3 mm thick slices of the 
object, with each slice centered on a pair of fiducial markers (Table 7.6). To minimize 
acquisition time, a T1 weighted spin-echo imaging sequence was employed, which 
featured a relatively short repetition time (TR). Acquired images were 256x128, 
corresponding to 256 frequency encoding steps and 128 phase encodings. MR images 
were interpolated to create a 256x256 image and then scaled to provide a pixel size 
equivalent to SPECT images. 

Table 7.6. Parameters Used for MR Phantom Imaging on CSI 2T System. 


MRJ Parameters 

Rod Phantom 

Egg Phantom 

FOV 

150 mm 

100 mm 

Slice Thickness 

3 mm 

3 mm 

Slice Separation 

6 mm 

3 mm 

Image Size 

256x128 

256x128 

Pixel Size 

0.585 mm 

0.391 mm 

Image Type 

T1-spin-echo 

T1-spin-echo 

TR 

600 ms 

1000 ms 

TE 

20 ms 

18 ms 

Elapsed Time 

2:33 min 

6:24 min 


Gadolinium (Gd) is a paramagnetic MR contrast agent which induces a large decrease 
in the spin-lattice relaxation time of surrounding tissue (or water). In T1-weighted spin 
echo imaging, the result is a very bright signal in areas where Gd is present. 
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Results 

Reconstructed SPECT and MR images of the rod phantom are presented in figure 
7.13. The images are shown without contrast enhancement. In figure 7.13A we observe 
a single MR image of the phantom with the fiducial markers clearly visible. The epoxy in 
the “cold” rod emits no MR signal but still provides a scattering medium for the SPECT 
measurements. The first SPECT image represents the rod reconstructed with filtered 
backprojection. It is observed that considerable noise exists in the image despite being 
reconstructed with over 50x10^ counts in the slice. Backprojection artifacts occur from 
the hot and cold rods, appearing as a spreading of activity in horizontal bands. Images 
7.13C and 7.13D show the same projection data reconstructed using the SAGE iterative 
algorithm. The difference between these images is the weighting applied to the 
penetration projection data. The data for the left image data are weighted by the 
Blackman window, and the right image data are weighted using the Monte Carlo based 
window. Both images are qualitatively superior to the FBP reconstruction but the 
Blackman windowed reconstruction has several dark line artifacts running through the 
object which are of indeterminate origin. Because these artifacts are effectively 
eliminated when using the Monte Carlo based weights, this window was used 
preferentially in the reconstruction of the three dimensional phantom as well as for animal 
imaging. 

The profile plots presented at the bottom of figure 7.13 represent two horizontal 
slices through figure 7.13D, centered on the hot rod and cold rod. The profiles 
demonstrate the relatively low noise provided by SAGE reconstruction. The finite 
resolution of SPRINT is also observed by the sloped edges of the hot and cold objects. 

The counts for the cold rod approaches zero within the object, supporting the simulation- 
based observation that scatter is a small effect. The “hot” rod has an absolute magnitude 
four to five times that of the surrounding activity, as expected. 
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Figure 7.13. MR and SPECT Images of Rod Source (A) MR Image, (B) FBP 
SPECT Image (C) SAGE w/ Blackman Weighting (D) SAGE w/ MC Weighting, 
and (E) Horizontal Profiles Through D. 
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Quantitative results for these images are presented in Table 7.7, with the applied 
calibration factor presented at the bottom. For the hot object, two assessments were 
made, one for the pixel of maximum counts and the other for a small region of interest 
(<20 pixels). The region-based assessment results in a small increase in error compared 
to the single pixel. The activity of the hot region is underestimated from the ROl analysis 
by 16% to 26%, and the background concentration is underestimated by 12%. The 
positive error of the cold rod is about five percent of the background concentration. The 
results are generally consistent when using either FBP or SAGE for image reconstruction, 
although SAGE provides slightly improved precision and accuracy. 

Table 7.7. Quantitative Analysis for 2-D Phantom, FBP and SAGE Reconstruction. 


FBP Reconstruction 

cpm/pixel ^ 

pCi/cc ±pCi/cc 

True 

Cone. 

pCi/cc 

Error (%)* 

Hot (pixel)^ 

327 

588 

703 

-16.2% 

(ROI) 

304 ±15.63 

548 ±90.2 


-22.1% 

Background 

72.6 ±21.8 

129 ±36.5 

136 

-4.9% 

Cold Rod 

5.48 ±7.28 

9.9 ±13.2 

0 

7.6% 

SAGE Reconstruction 

cpm/pixel 

pCi/cc ±pCi/cc 

True 

Cone. 

pCi/cc 

Error (%)’ 

Hot (pixel)^ 

43770 

543 

703 

-22.73% 

(ROI) 

41590 ±1422 

516±83.1 


-26.58% 

Background 

9631 ±1397 

119 ±25.5 

136 

-12.03% 

Cold Rod 

449 ±466 

5.57 ±5.85 

0 

4.66% 


Calibration Factors 

FBP: 

0.555 ±0.087 


(cpm/pixel/ixCi/ml) 

SAGE: 

80.6 ±12.6 


* Value presented for cold rod is fraction of background activity measured. 

^ No standard deviation reported because a single (maximum) pixel is used. 

^ Standard deviation equal to the variance in average pixel counts within single ROL 
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The single measurement standard deviation in the assessed background 
concentration is high. Although this effect can be reduced by increased filtering of the 
image, this will result in increased error. The cold rod has a standard deviation in excess 
of 100% for both reconstruction methods attributable to low count statistics in this region. 
This is reasonable, considering that if there were no biasing due to scatter, the variance in 
the cold region would be infinite. 

The imaging conditions for the two-dimensional phantom are considered ideal in 
that high activity was used in objects large compared to the resolution of the system, and 
the object did not vary axially. This removes the effects of a finite axial resolution. A 
more realistic assessment of SPRINT image performance is obtained by imaging the 3-D 
phantom over a range of object activities. Image quality for the highest activity 
acquisition is shown in figure 7.14 with matched MR images for reference. Slices are 
contiguous with a 3 mm thickness. MR and SPECT axial alignment is approximate, as 
evidenced by only partial agreement of the fiducial markers between the matched images. 
The fifth image of the sequence (slice E) is taken as the center of the object for purposes 
of quantification. 

Each of these SPECT images was reconstructed from approximately 8 million 
open slit projection counts and 4 million blocked slit counts using the SAGE algorithm 
and regularization of P=2' , chosen as the best compromise between image resolution 
and noise for this acquisition. Despite high image counts, considerable noise is observed 
in the images. The figure also shows the effect of the limited axial resolution provided by 
the slice collimator. Activity from the “hof ’ egg and “cold” egg is seen to blur slightly 
into slices above and below the actual object. 
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Figure 7.14. MR and SPECT Images of 3-D Egg Phantom. Images are Contiguous 
3 mm Thick Slices. Slice E is Approximately in the Center of the Phantom. 

The effect on image quality of reducing total image counts is shown in figures 
7.15 and 7.16. The first figure shows FBP reconstructions of the center slice of the 
phantom, and the second figure shows SAGE reconstructions of the same slice. All 
images were reconstructed using the parameters shown in Table 7.5. The top row in each 
figure shows the effect of a 2 hour decay period between 55 minute acquisitions. This 









152 


results in neighboring images having total image counts reduced by a factor of 2.6. The 
lower sets of images were acquired approximately one hour after the start of the overlying 
image, and using a 10 minute acquisition time. This reduced the image counts between 
top and bottom images by an approximate factor of 10. We observe that at the end of the 
imaging sequence, the hot and cold objects are just discernible against the high 
background noise of the image. 

A requisite for accurate quantification of phantom objects and animal organs is 
that the system be linear over a wide range of object activities. To assess the linearity of 
SPRINT, total projection counts and total image counts for the center phantom slice were 
plotted as a function of total phantom activity times the image acquisition time. This 
latter function is directly proportional to the total flux emitted from the phantom during a 
given image acquisition. The results are shown in figure 7.17. The figure shows that for 
the first three acquisitions (as labeled), the number of image counts and projection counts 
are consistently less than expected based on a linear fit of the last four acquisitions when 
the phantom activity is lowest. The bias is better observed in table 7.8 which shows 
percent difference in measured counts compared to a linear fit of the data. The effect is 
believed to originate from paralyzable deadtime effects in module signal processing 
which are not addressed by the simple non-paralyzable model described previously. 

Results of the region of interest analysis are shown in Table 7.9. The table shows 
the average count rate per pixel for the various phantom regions reconstructed using the 
SAGE algorithm. Results using EBP were found to be of poorer accuracy and precision 
than those obtained using SAGE, and are not presented. The hot egg concentration was 
determined from both a small region of interest (<40 pixels) as well as the pixel of 
maximum counts within the object. Background concentrations were assessed using 
large regions (several hundred pixels) averaged over the center slice and slices above and 
below the hot and cold objects. The cold region was assessed from a small region drawn 
within the cold egg boundaries. 



Figure 7.15. Serial Acquisitions of Center Slice of Egg Phantom Reconstructed with 
FBP and Butterworth Filtering (0.9 cm'^). Adjacent Images are Decreased in 
Counts by 2.6, and Bottom Images are Factor of 10 Less in Counts Than Overlying 
Images. Hot egg Concentrations vary from 853 |aCi/ml to 28.4 |aCi/ml. 



Figure 7.16. Serial Acquisitions of Center Slice of Egg Phantom Reconstructed with 
SAGE and p=2'^^. Phantom Activity and Concentrations are in Tables 7.8 and 7.9. 
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Figure 7.17. Total Image Counts as a Function of Total Flux Emitted by the 
Phantom. The Product of Total Phantom Activity and Time the Phantom was 
Imaged (fxCi/hr) is Proportional to Total Flux. Note the Bias for the First Three 

Acquisitions. 


Table 7.8. Bias Introduced by Uncompensated System Deadtime. 



Total 

Phantom 

(mCi) 

Acquire 
Live Time 
(min) 

Net Projection 
Counts (Op-Bl) 
(Counts X 10^) 

Bias from 
Linear 

SAGE 

Image 

(Counts X 10^) 

Bias from 
Linear 

1 

30.26 

27 

5430 

-21.6% 

1060 

-20.7% 

2 

21.26 

4 

572 

-16.7% 

111 

-15.7% 

3 

11.39 

27 

2040 

-5.16% 

398 

-5.49% 

4 

7.10 

4 

191 

-2.06% 

37.2 

-1.74% 

5 

4.29 

27 

769 

-0.14% 

150 

0.06% 

6 

2.81 

4 

75.7 

0.87% 

14.7 

2.82% 

7 

1.62 

27 

290 

0.60% 

56.4 

-2.61% 

8 

1.01 

4 

27.1 

-4.35% 

5.27 

0.08% 
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The activity concentration was assessed in these regions using the background 
derived calibration factor. This factor was based on images acquired with phantom 
activities less than 7 mCi to avoid the paralyzable deadtime effects discussed previously. 
Results were compared to the known object activity concentration at each acquisition 
time, providing the error in the single measurement as a percent difference. The error in 
the cold egg concentration was determined as a fraction of the background concentration. 

The most notable observation made from these data is the change in bias of the 
hot egg and background regions as the activity within the phantom decays. When the 
phantom is at its highest activity, there is a -28% error in the assessed hot egg 
concentration. When total phantom activity diminishes to less than 7 mCi (260 MBq), 
the average error approaches 3% for the hot egg. This same trend is noted for the ROI 
analysis of this region, but activity concentration is consistently underestimated, even for 
low phantom activities. Similarly, the background error is initially -20%, which then 
increases to a residual positive error of -3%. There is generally good agreement between 
the initial acquisition results for the 3-D phantom, and the single acquisition result for 
2-D phantom. The measured activity concentration within the cold egg shows a fairly 
uniform activity concentration that is 3% to 20% of the background. 

As with the two-dimensional phantom, it must be recognized that the standard 
deviation presented for each measurement is not the actual variance of the measurement. 
This is again because the reconstruction process and effects of filtering or regularization 
correlate the image pixels. No standard deviation for the count rate of the hot egg is 
reported since only a single pixel value was used. The measure is observed to increase as 
source concentration diminishes and acquisition time is reduced (reference Table 7.8 for 
true acquire times), as expected. The standard deviation in the cold region is again 
extreme, approaching several hundred percent in some instances due to low counts. 
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Table 7.9. Quantitative Analysis for 3-D Egg Phantom, SAGE Reconstruction. 


Hot Egg Region Region count rate 

Measured Concent. 

True Cone 

Error (%)* 

Acquisition 

cpm/pixel ±cpm/pixel^ pCi/ml ±pCi/ml 

pCi/ml 


1 (max pixel) 

24600 ^ 

611 

853 

-28.3% 

(ROI) 

24000 ±530 

595 ±94.2 

853 

-30.2% 

2 (max pixel) 

19400 

482 

599 

-19.6% 

(ROI) 

17600 ±799 

437 ±71.2 

599 

-27.2% 

3 (max pixel) 

10984 

273 

321 

-15.14% 

(ROI) 

10500 ±218 

261 ±40.2 

321 

-18.64% 

4 (max pixel) 

7660 

190 

200 

-5.01% 

(ROI) 

7070 ±477 

176 ±30.0 

200 

-12.32% 

5 (max pixel) 

4520 

112 

121 

-7.21% 

(ROI) 

4330±110 

107 ±17 

121 

-11.20% 

6 (max pixel) 

3580 

88.7 

79.3 

11.93% 

(ROI) 

3000 ±357 

74.6 ±14.6 

79.3 

-5.88% 

7 (max pixel) 

1690 

41.8 

45.5 

-8.11% 

(ROI) 

1540 ±96 

38.3 ±6.44 

45.5 

-15.97% 

8 (max pixel) 

1450 

36.0 

28.4 

26.74% 

(ROI) 

1000 ±338 

24.9 ±9.26 

28.4 

-12.43% 

Background Region 

1 

5383 ±676 

134 ±27.0 

167 

-20.0% 

2 

4114±1286 

102 ±35.4 

117 

-12.9% 

3 

2297 ±327 

57 ±12.1 

63 

-9.30% 

4 

1525 ±702 

38 ±18.4 

39 

-3.43% 

5 

950 ±193 

24 ±6.05 

24 

-0.41% 

6 

621 ±432 

15.4±11.0 

15.5 

-0.70% 

7 

361 ±121 

8.9±3.31 

8.9 

0.41% 

8 

204 ±220 

5.1 ±5.52 

5.6 

-8.86% 

Cold Egg Region 

1 

680 ±425 

16.9 ±10.9 

0 

11.7% 

2 

371 ±562 

9.21 ±14.0 

0 

8.9% 

3 

492 ±190 

12.2 ±5.12 

0 

19.6% 

4 

49.0 ±130 

1.22 ±3.23 

0 

2.96% 

5 

173 ±108 

4.29 ±2.75 

0 

17.4% 

6 

73.3 ±136 

1.82 ±3.38 

0 

10.3% 

7 

59.8 ±59 

1.48±1.48 

0 

18.3% 

8 

14.1 ±40 

0.35 ±0.99 

0 

5.83% 

Calib. Factor: 

40.3 ±6.30 

cpm/pixel/pCi/ml 




* Value presented for cold rod is fraction of background activity measured. 
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^ No standard deviation reported because a single (maximum) pixel is used. 

^ Standard deviation equal to the variance in average pixel counts within single ROI. 

Estimations of the bias and standard deviation (STD) for the regions of the three 
dimensional phantom are presented in Table 7.10. To avoid the effects of unrecovered 
deadtime, only the last five acquisitions were used in the determination. Both weighted 
least squares, as discussed earlier, and unweighted least square estimates are presented. 


Table 7.10. Bias and Standard Deviation Assessed for Hot Egg and Uniform 
Background of Three Dimensional Phantom. Weighted and Unweighted Least 
Square Estimates are Presented (n=5). 


1^,0 

Qknown 

% Bias 

STD 

%STD 

Hot Egg 

(pCi/ml) 

(p,Ci/ml) 


(qCi/ml) 


Weighted 

843 

853 

-1.14% 

132 

±15.7% 

Unweighted 

884 


3.67% 

130 

±14.7% 

Weighted^ 

756 

853 

-11.3% 

143 

±18.8% 

Unweighted 

754 


-11.6% 

148 

±19.7% 


Background 

Region 


Weighted 

158 

167 

-5.04% 

9.72 

±6.13% 

Unweighted 

161 


-3.72% 

6.65 

±4.14% 

Cold Egg 

Weighted 

16.5 

0 

12.3%^ 

31.3 

±190% 

Unweighted 

18.7 


18.3% 

23.5 

±126% 


Based on maximum pixel 
^ Based on region of interest 
^ Reported as percent of uniform background 
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Discussion 

Results obtained from phantom imaging show the limitations of high energy pin¬ 
hole or slit aperture SPECT. Filtered backprojection reconstructions show very high 
noise and artifact characteristics. This is due to the low intrinsic detector sensitivity 
(11%) and the use of high resolution collimators, the large penetration component of the 
projection data and the projection subtraction technique employed to correct it. 

Additional noise is introduced in the FBP reconstructions from the arithmetic operations 
of back-projection*^ and the use of almost a pure ramp filter that enhances high frequency 
noise. Because projection data are not weighted by their statistical significance (as is 
done in SAGE) noisy projection introduce streak artifacts in FBP images. Image 
reconstruction performed using the regularized SAGE algorithm provides significant 
improvements in image quality over FBP, but limitations posed by system sensitivity still 
require use of extremely high activities to provide high quality images. Although 
phantom images are noisy, they also demonstrate very good resolution for this energy. 
This is achieved through high resolution collimation, and a magnifying imaging geometry 
combined with the application of aperture penetration correction. 

Table 7.9 indicates that for uniform background regions, a limiting activity 
concentration which could be imaged for 26.7 minutes and provide images with less than 
a 10% standard deviation in pixel counts exceeds 160 pCi/ml (5.6 MBq/ml). By 
comparison, all 4 minute acquisitions provided images with background standard 
deviations greater than 30%. Clearly, obtaining these activity concentrations routinely in 
small animals for all target organs of interest would require prohibitive levels of 
radioactivity to be used. 

Examining the error in each of these measurements indicates that for high 
phantom activities, the activity concentration is underestimated on the order of 20% to 
30%. Therefore, the best qualitative images are acquired at the cost of using very high 
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activities and result in a negative bias in quantified hot objects originating from deadtime 
effects not corrected by a non-paralyzable model. As total phantom activity decays to 
less than 7 mCi, the residual error becomes small but noise in the image becomes large. 

The bias and standard deviation observed in Table 7.10 indicates that it is possible 
to measure the activity concentration in large organs (relative to system resolution) with 
fair precision and good accuracy. It is recognized that this particular assessment does not 
have high statistical power, but experimental conditions limited the number of samples 
that could be acquired. 

The measured activity concentration was found to be dependent on numerous 
factors. Use of a line source based calibration factor was found to result in an 
overestimate of measured activity concentrations by up to 14%. The effect is attributed 
to errors in the assumed pixel volume, which is required when assessing concentrations 
using the line source based calibration factor. Further, it was found that if background 
measurements are made in the central slice of the phantom, the background is 
overestimated by about 10%. This is thought due blurring and scatter of the hot region 
activity into the neighboring background region. 

In assessing the hot egg concentration, the use of other than the maximum count 
pixel in the region causes in an underestimation of the concentration by 10%. The egg 
object had linear dimension of 14 mm by 24 mm in-plane and 14 mm axially (as 
positioned in the phantom). In-plane resolution of SPRINT, assessed using SAGE 
reconstruction, was 4.4 mm FWTM and 8.0 mm FWTM. Axial resolution was 
approximately 5 mm FWHM and 8 mm FWTM, but with considerable variance. We note 
that for these dimensions, even small (4 to 6 mm diameter) regions of interest will be 
affected by the limited axial and in-plane resolution. The negative bias observed for this 
object thus originates from the blurring of activity out of the finite sized ROI. It is also 
likely that scatter out of the region contributes to the observed negative bias, since the 
method adopted for attenuation and scatter correction was not ideal. The hottest central 
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pixel of the region serves as the best estimator of region concentration (minimum bias) 
because it is located at approximately the systems FWTM from the object boundary and 
thus experiences the least impact from the systems limited resolution. 

Results for the cold object in each phantom show the effects of scatter as well as 
the limited axial resolution provided by the high energy slice collimator. For the rod 
phantom where the object is uniform axially, contribution from scatter causes a measured 
concentration approximately 5% of the surrounding background. This is in good 
agreement with Monte Carlo simulation results presented earlier. For the 3-D phantom, 
limited axial resolution contributes additional activity to this region. The resulting error 
of 10% is thus reasonable. In theory, a better estimate of cold region activity would be 
achieved by only measuring the single central pixel of the object. However, the 
significant noise demonstrated in this region made such a measure unreliable. For all 
cold region measurements, the high noise is again attributable to the low count 
characteristics of this area. 
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Quantitative Rat Organ Imaging 

To establish the accuracy and precision of SPRINT as applied to high energy 
animal imaging, two studies were conducted using Sprague-Dawley rats injected with 
[‘^FJfluorodeoxyglucose (^*FDG). Images were acquired of the head and thoracic regions 
of each animal, and quantification was attempted for the heart, brain and lungs. Results 
were compared to invitro biodistribution measurements obtained after imaging 

Methods 

Two female Sprague-Dawley rats with weights of 250 to 275 gm were sedated 
with 40 mg/kg pentobarbitol under ether anesthesia. After a 10 minute waiting period, 
the animal was loaded into a specially designed animal holder compatible with SPRINT 
and MR imaging. The holder was fitted with Teflon nose and tail rings to hold the animal 
securely during imaging. Fiducial markers along the holder walls (10 pi volume, 8 mm 
axial separation) were filled with 2 mM Gd-DTPA. A picture of the sedated animal 
prepared for MR imaging is presented in figure 7.18. 

The animal was then imaged using the CSI2T MR scanner described previously. 
Image acquisition parameters are presented in Table 7.11. As with the phantom, a special 
MR probe was constructed to precisely match the outer diameter of the animal holder in 
an effort to maximize the detected RF signal. The animal was positioned in the receiver 
probe with the linear array of fiducial markers facing upward, with the first fiducial 
marker located just within the sensitive field of the probe. The animal and probe were 
then inserted into the MR coil. A coronal image was acquired to assess slice location, 
and then a Tl-weighted multi-slice spin-echo imaging sequence was performed. Figure 
7.19 depicts the animal within the CSI 2T, and supported by the receiver probe. 
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After MR imaging, [^*F]fluorodeoxyglucose was obtained from the University of 
Michigan PET Facility. The available specific concentration dictated the maximum 
radioactivity that could be injected in the animal. One ml of solution was prepared for 
each animal, allowing injection 13 mCi for the first animal, and 3.84 mCi for the second. 
Injected doses and residual syringe activity were assessed using the Capintec dose 
calibrator. After animal injection, fiducial markers were emptied of the Gd-DTPA 
solution and filled with residual ^^FDG, resulting in 66 pCi per marker for the first 
animal, and 100 pCi per marker for the second. A thirty minute waiting period was 
permitted between injection and the start of SPRINT acquisitions to allow metabolism of 
the ^*FDG. The animal was then placed within SPRINT using a holder that precisely 
centered the animal within the gantry. The mounting system can be observed in figure 
3.2. Imaging was then performed in a manner analogous to the phantom studies using the 
parameters shown in Table 7.11. Two images were acquired of each animal, using a 55 
minute and 10 minute acquisition time. An animal (tail) during an acquisition is observed 


in figure 7.20. 

At two hours post injection, the animal was sacrificed and biodistributions 
measured for the brain, heart, lung, liver and blood. Total organ weights were recorded 
for the heart and brain. Organ activities were measured using a Packard Tricarb 578 
autogamma system with a 3”x3” Nal(Tl) detector. Results were reported as %dose per 
gram for all organs, and %dose per organ for the heart and brain. From the known 
injected activity, the time of imaging, and assuming the biodistribution was constant 
between imaging and sacrifice, the total organ activity and activity concentration at the 


time of imaging were estimated: 

rr ^ 1 A %Dose Inj.Actmty(jua) 

Total Actmtyi jsCi) -- 

Organ- 100 


- ln( 2 )*DecayTime(min) 


ActivityConcentration( juCi /ml) 


%Dose Inj.Aclmly(pCi) 
100 


gm- 
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Table 7.11. Acquisition Parameters for Rat MR and SPRINT Imaging. 


MRI Parameters 

SPRINT Parameters 

Image Type 

T1-spin-echo 

No. Steps 

80 

TR 

600 ms 

Time/Step 

40 s (1,3) 6 s (2,4) 

TE 

20 ms 

Total Time 

55 m(l,3) 10m(2,4) 

NEX 

2 

Live Time 

26.67 m and 4 m 

Phase Encoding 

128 

Reconstruction: 

SAGE/FBP 

Frequency Encoding 

256 

Image Size 

100x100 

FOV (mm) 

80x80 

FOV (mm) 

86x86 

Pixel Size 

0.3125 mm 

Slice Thick: 

3 mm 

Acquire Time 

2.33 m 

Sum Interval 

2 

Slice Thickness 

3 mm 

No. Slices 

14 

Slice Separation 

8 mm 

Beta/ Convergence 

2'^’ and 2-^Vo.Ol 



Figure 7.20. Animal (Tail) in SPRINT During SPECT Acquisition. 
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Image reconstruction was performed using filtered backprojection and the SAGE 
iterative algorithm. Image analysis was conducted by drawing regions of interest around 
and within the brain and heart, and within the approximate regions for the lungs. Three 
measures were obtained for the brain and heart: the total counts in a large region 
surrounding the organ, counts per pixel for small regions within the organ, and the total 
counts of the maximum pixel. This allowed assessment of total organ activity and organ 
activity concentration using the calibration factors derived from the 3-D phantom and line 

source, and the following expressions: 

, ^ X ^ / ,> ROI Specific Countratef^cp/nZ/j/ve/^ 

KadioactivityConcentration( uCi/ ml) =-;-—-;—;— 

Cal Factor(^ cpm/pixel/ uCi/ml) 


Total Activity( uCi) 


ROI Countrate(^cpwj 
Cal Factor cpm / uCi) • 


For the heart and brain, measurements were made in the two or three contiguous 


6 mm slices where activity from these organs could be observed. For the brain and lung, 
the standard deviation in average pixel counts was measured from the region of interest. 
This again provided a relative indication of noise in these regions. From the activity 
concentration and total activity of these organs determined by invitro measurements, the 
error (as %difference) in the SPECT measurement was calculated. 


Results 

SAGE reconstructed images are presented in figure 7.21 for the animal injected 
with 13 mCi and imaged for 55 minutes, with matched MR images to the left. The rodent 
had an approximate total body activity of 10 mCi at the time of imaging. Images are 
presented without contrast enhancement. Because of slice separation, there are three MR 
images presented for every four SPECT images. MR images were only acquired for the 
head and upper thoracic region. SPECT images were acquired using a larger field of 
view to include the cardiac region, as observed in figure 7.22. 
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Figure 7.22. SPECT Images of Rat Cardiac Region, 6 mm Thick Slices Continued 

from Figure 7.21. 


Each image shows well visualized fiducial markers, one superior to the animal, 
and the other rotating in a counter-clockwise fashion around the animal holder. The 
central regions of the rat brain are observed in SPECT slices D, E and F of figure 7.21, 
assessed principally on the anatomical information provided by the matched MR images. 
Determination of brain location exclusively from SPECT images leads to the erroneous 
conclusion that the brain is principally in slice C. This high activity region will be 
discussed in greater detail in the following section. Since MR image information was 
not available for the thoracic cavity, cardiac regions of interest were based purely on 
SPECT image contrast. 

Contrast and resolution of the MR brain images allows differentiation of the 
cortical structure, corpus callosum, and left and right caudate putamen (slice E). In slice 
C, the left and right olfactory bulbs are seen and in slice F the cerebellum is found with its 
pyramids and parafloccular bulbs. Matched SPECT images show relatively uniform 
uptake in the these regions. In contrast, cardiac uptake was quite high in both animals 
studied, while the assessed lung uptake was very low. This results in the appearance of a 
well defined cardiac region in slices K and L of figure 7.22. These images are continued 
spatially from the last SPECT image of figure 7.21. Limited axial resolution and partial 
volume effects permit cardiac activity to be visible in slice J of figure 7.21 and M of 
figure 7.22. Cardiac motion prevents any anatomical structure (e.g., left ventricle) from 
being observed in these images. 
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The central SPECT slice for the brain and heart for all four animal acquisitions is 
presented in figure 7.23. Each SAGE image has been contrast enhanced to overcome the 
suppression in detail caused by the hot fiducial markers. The first two images of each 
row are for the animal injected with 13 mCi ’®FDG, and imaged for 55 minutes and 10 
minutes, respectively. The second two images are for the second animal injected with 
3.84 mCi and imaged for the same periods. Table 7.12 shows, for each of these images, 
the invitro assessed organ activity and activity concentration decay corrected to the 
acquisition time, and live acquire time used for each acquisition. Comparing the images 
to the tabulated data, the effects on image quality of decreased acquisition time and organ 
content are observed. 



A. Brain Images, Acquisitions 1-4 



B. Cardiac Images, Acquisitions 1-4 


Figure 7.23. Contrast Enhanced SPECT Images of Rat Brain (Top Row) and 
Cardiac Regions (Bottom Row). First Two Images: Rat Injected with 13 mCi 
^*FDG, Second Two Images: Rat Injected with 3.8 mCi ^*FDG. 
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Table 7.12. Measured Invitro Biodistribution and Estimated Organ Activity 
Concentrations for the Rat Heart and Brain. Each Column Shows the Decay 
Corrected Organ Concentration and Total Activity of FDG at Start of an Image 
Acquisition, Based on the Invitro Measured Value. 




Animal 1, imaged 30 and 90 minutes 
post injection with 13 mCi '^FDG 

Animal 2, imaged 30 and 90 minutes 
post injection with 3.84 mCi **FDG 



Injected; 

Acquired: 

Acqu.Time 

13:20 

11:04 12:06 

55 min 10 min 

14:18 

14:46 15:42 

55 min 10 min 


Invitro 

Biodistribution 

%dose/g 

Decay Corrected^ 
Organ Concentration 
pCi/g 1 

Invitro 

Biodistrib. 

%dose/g 

Decay Corrected^ 
Organ Concentration 
pCi/g 

Brain 


0.840 

88.1 59.6 

1.20 

38.5 27.1 

Heart 


6.14 

644 436 

3.65 

117.4 82.5 


Invitro 

Biodistribution 

%dose/org 

Decay Corrected^ 
Organ Activity 
pCi pCi 

Invitro 

Biodistrib. 

%dose/org 

Decay Corrected^ 
Organ Activity 
pCi pCi 

Brain 


1.68 

176 119 

2.31 

74.2 52.2 

Heart 


5.60 

588 398 

2.92 

94.0 66.1 


No precision is reported for biodistribution measurements since only one animal was 
used in each assessment. 

^ Invitro biodistribution data is decay corrected to the start of image acquisition. 


Results of the quantitative assessment of these images is presented in Table 7.13. 
For each organ, activity concentration (pCi/ml) and total activity (pCi) is presented. 
Results are compared to the invitro measured organ activity and concentration to provide 
a percent difference measure of error. The total activity assessed in the brain and heart 
had modest errors. Similarly, the activity concentration for the brain varied from -3% to 
+7% The concentration in the heart was assessed from the pixel of maximum counts. 
This value still resulted in the estimated organ concentration to be approximately 50% 
less than the invitro measured organ concentration. For the assessed lung concentration, 
results were inconsistent and had errors ranging from ±20% between animals. A standard 
deviation in assessed organ concentration is presented only for the lung and brain, where 
the variance in measured pixel counts for relatively large regions is reported. 
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Table 7.13. Quantitative Image Analysis of Rat Heart, Brain and Lnng. 


Activity 

Small region count rate 

Calculated 

Estimated^ 

Error 

Concentration 

(brain & lung) and 

Concentration 

Concentration 

(%) 

Assessment 

single hot pixel (heart) 

(ROI cpm/pixel/CF) 



Brain 

cpm/pixel ±cpm/pixel 

fiCi/cc ±|j,Ci/cc 

pCi/cc 


Acquis. 1 

7480 ±1030 

93 ±15.4 

88.1 

5.40% 

Acquis. 2 

5380 ±685 

66.8 ±10.4 

59.6 

12.1% 

Acquis. 3 

2880 ±244 

35.7 ±4.43 

38.5 

-7.41% 

Acquis. 4 

2110±177 

26.2 ±3.23 

27.1 

-3.36% 

Heart 

cpm/pixel 

pCi/cc 

pCi/cc 


Acquis. 1 

29300 

364 

644 

-43.5% 

Acquis. 2 

23200 

287 

436 

-34.0% 

Acquis. 3 

4170 

55.2 

117 

-53.0% 

Acquis. 4 

3880 

53.6 

82.5 

-35.0% 

Lung 

cpm/pixel ±cpm/pixel 

pCi/cc ±pCi/cc 

pCi/cc 


Acquis. 1 

2780 ±556 

34.5 ±7.58 

28.9 

19.6% 

Acquis. 2 

1760 ±439 

21.9 ±6.79 

19.5 

11.9% 

Acquis. 3 

1230 ±215 

15.2 ±3.00 

18.8 

-19.1% 

Acquis. 4 

837 ±260 

10.4 ±3.36 

13.2 

-21.4% 

Calibration Factor (CF): 80.6±12.6 cpm/pixel/pCi/ml 


Total Activity Large region count rate^ 

Calculated Activity 

Estimated^ 

Error 

Assessment 


(ROI cpm/CF) 

Organ Activity 

(%) 

Brain 

ROI cpm 

pCi 

pCi 


Acquis. 1 

3040000 

189 

176 

7.28% 

Acquis. 2 

1930000 

120 

119 

0.73% 

Acquis. 3 

1160000 

72.2 

74.2 

-2.73% 

Acquis. 4 

857000 

53.6 

52.2 

2.65% 

Heart 

ROI cpm 

pCi 

pCi 


Acquis. 1 

9660000 

600 

588 

1.99% 

Acquis. 2 

7390000 

400 

398 

0.57% 

Acquis. 3 

1430000 

89.5 

94 

-4.81% 

Acquis. 4 

1140000 

71.2 

66 

7.73% 

Calibration Factor (CF): 16100 ±600 cpm/pCi 


* ROI analysis based on small to moderate size regions within the organ boundaries to 
provide an estimate of organ concentration. 


^ ROI analysis based on large regions which completely encompass the observed organ. 
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^ Estimated organ activity concentration based on invitro measured organ concentration at 
time of sacrifice, decay corrected to start time of imaging. 

^ Estimated organ activity based on invitro measured organ activity (%dose/organ) decay 
corrected to the start time of image acquisition 

Discussion 

For the limited number of acquisitions made, there was generally good agreement 
between the invitro biodistribution measurements of the brain and heart, and those 
measured from SPECT images. An obvious exception is the assessed concentration of 
activity in the rat heart, where the invivo result is half that expected. This is attributed to 
heart motion and limited SPRINT resolution which combine to decrease the estimate of 
myocardial concentration. Specifically, the myocardial wall thickness of 1 to 3 mm is 
blurred by motion of the heart wall. This blurred concentration is further distributed in 
the reconstructed image by the four to five mm resolution of SPRINT. The consequence 
is a significant underestimate of the myocardial concentration compared to invitro 
measurements. The inconsistent results for the lungs are attributed to the lack of 
anatomical information available when drawing regions of interest. Since lung regions 
were estimated purely from cardiac SPECT position, it is likely that region counts were 
affected by activity in the mediastinum or even by including areas outside of the animal. 

One of the purposes of this work was to determine total organ activity invivo to 
determine organ time-activity profiles more accurately. The error observed in organ 
concentration measurements are compounded if an inaccurate determination of organ 
mass or volume is made. Since it appears possible to accurately assess the organ activity 
concentration for organs having a substantial uptake, invivo measurements can be used to 
determine the time-activity biodistribution pattern of a given compound. The total organ 
volume can then be assessed ex-vivo following the imaging sequence. By measuring the 
biodistribution of the organ or tissue at the end of imaging, an organ specific calibration 
factor be also be derived using image data from the last acquisition. 
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1 R 

For both the hrain and heart, significant local uptake of FDG was experienced 
compared to surrounding tissue. Brain uptake may have in fact been mitigated by the use 
of CNS depressant to anesthetize the animal. In these cases, large regions of interest that 
completely encompass the organ provide an estimate of total activity that appears to he 
close agreement to invitro measurements. Success was dependent on summing activity in 
all slices where the organ was present, and included three contiguous slices for the brain 
and two slices for the heart. Unfortunately, with the possible exception of the rodent 
testes, the brain and the heart are the only organs which could be assessed in this manner 
due to their location and high glucose metabolism. Determining the biodistribution of 
organs which have significant activity surrounding them (e.g. organs in the abdominal 
region), must be assessed from small regions within the organ to estimate activity 
concentration. 

As was the case for phantom studies, significant noise and limited resolution 
compromise the ability to quantify organ concentration. For the organs measured, ROI 
based standard deviations were reasonable, considering the relatively high activities 
present. However, it is evident that assessing all significant animal organs for uptake of a 
given compound is not possible. Limited resolution and poor sensitivity effectively limit 
accurate and precise assessment to animal organs that are large compared to system 
resolution and demonstrate substantial uptake of the compound. 

It is again emphasized that the reported errors and standard deviations in Table 
7.13 are only approximate indicators of the accuracy and precision of this measurement 
process due to correlation between image pixels. Although improved estimates of the 
bias and variance could have been obtained for the organs examined by imaging a 
sequence of animals, this was considered beyond the scope of the dissertation. 
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Acute ^®FDG Uptake by the Harderian Gland 

A casual observation of figure 7.21C, without anatomical reference provided by 
the MR images, would lead to the conclusion that fluorodeoxyglucose uptake in the brain 
was significant and well differentiated between hemispheres. The two bright foci could 
easily be interpreted as local uptake in the parietal cortex or caudate putamen. It is only 
after the fiducial marked MR images are overlaid with the SPECT images that these foci 
appear to be lateral and superior to the animal brain. It is noted that SPECT slice C lies 
between adjacent MR slices B and C. Figure 7.24 presents the effect of splitting slice C 
into two adjacent 3 mm slices, and overlaying the peak activity regions onto the adjacent 
MR slices (figure 7.24 E and F). We observe in MR slice A, the right eye and a white 
region that is medial and inferior to both eyes. Because the MR data was acquired with 
T1 weighting, this region is recognized as probably having a high fat content. The 
MR/SPECT overlaid image indicates that activity is concentrated in this area, and that it 
extends posterior to this region. Gross dissection of a rodent head suggests that this high 
regional uptake may be attributed to the Harderian gland. This gland is horseshoe shaped 
and encircles the optic nerve, extending medially from the eye to the dorso-lateral region 
of the orbit^^. The gland is found in most mammals, and is particularly well developed in 
swimming mammals, responsible for an oily secretion to lubricate and protect the eye*’. 
This is in good agreement with both MR position and composition information. 

To verify that the Harderian gland has preferential uptake of **FDG or a fluorine 
labeled metabolite, a biodistribution study was performed on three Sprague-Dawley rats. 
Animals were prepared in a fashion identical to the imaging studies and injected with 
approximately 80 p,Ci **FDG each. At two hours post injection the rat brain, pituitary 
gland, olfactory bulbs, eye, temporalis muscle and Harderian gland were dissected and 
measured for radioactivity content as described previously. 
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Figure 7.24. Assessment of Focal ‘*FDG Uptake Medial, Inferior and Posterior to 
Rodent Eye. Figure A and B: Adjacent 3 mm Thick MR Images with 8 mm 
Separation; Figure C and D: Contiguous 3 mm SPECT Images Between MR 
Images. Figure E and F: Overlay of SPECT Peak Information onto MR Images. 
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Results of the study indicated that the Harderian gland had a biodistribution of 
5.39±0.19 %dose/gm (1.29±0.09 %dose/organ), compared to 1.13±0.14 %dose/gm 
(2.14±0.22 %dose/organ) for the brain. The small organ mass of this gland (0.12 
gm/gland) results in the bright foci observed in the SPECT images. ROI assessed gland 
activity for the two imaged animals results in fair agreement with these biodistribution 
measurements, with an error of ±14%. Specific activity determined from ROIs was a 
factor of two lower than the invitro results. The differences can be attributed to many 
factors including: differences in biodistribution for the imaged animals and those 
assessed invitro; and finite resolution of SPRINT compared to the size of the gland. 

The acute uptake expressed by the Harderian gland emphasizes the significant 
physiological and anatomical differences between research animals and humans. Side 
anatomical information is required during the invivo assessment of organ activity, not 
only to assist in region of interest analysis, but to differentiate and reveal physiological 
data that may confound measurements and affect their application to humans. 
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Comparison of Image Performance with Positron Emission Tomography 

It is useful to provide a baseline by which the qualitative and quantitative results 
obtained with SPRINT can be compared. This was accomplished by imaging a simple 
uniform phantom and a single rat injected with FDG using a conventional two- 
dimensional positron emission tomography (PET) system. The PET scanner used was a 
model PCT4600 produced by the Cyclotron Corporation, which is available for animal 
imaging. Images were acquired of a line source and a cylinder of uniformly distributed 
activity to determine system sensitivity and image noise. Images were then acquired of a 
Sprague-Dawley rat injected with [ Fjfluorodeoxyglucose. 

Methods 

The PET scanner featured three rings of bismuth germinate (BGO) detectors, 
forming three slices axially. Individual rings and adjacent rings are in coincidence, 
allowing five slices to be reconstructed with a 10 mm slice thickness and an 11.5 mm 
slice separation. Sensitivity and resolution were determined by imaging an F line 
source with a lineal activity of 0.84 pCi/mm for 15 minutes. Sensitivity was verified and 
image uniformity was evaluated using a cylinder phantom filled with a imiform activity 
concentration (1.85 pCi/ml) and again imaged for 15 minutes. Projection data was 
reconstructed using filtered backprojection with a pure ramp filter with a 1.8 cm'* cutoff, 
equal to the Nyquist frequency of the system. Corrections for non-uniformity, deadtime, 
decay, attenuation, and wobble"'' were applied during reconstruction. 

Phantom imaging was followed by imaging the brain and thoracic region of a 275 
gm female Sprague-Dawley rat. The animal was anesthetized with 40 mg/kg 
pentobarbital and then placed in the animal holder described previously. Fiducial 

The detector ring is wobbled during acquisition to increase the spatial sampling of the 
system. 
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markers were filled with 4 mM Gd-DTPA. T1-weighted spin-echo MR images were 
acquired of the animal using the CTI 2T system. Acquisition parameters are shown in 
Table 7.14. The rat was then injected with 1170 pCi of ‘^FDG, and distribution was 
permitted for 30 minutes during which time the fiducial markers were replaced with 
*^FDG (~1 pCi/marker). The animal was imaged for three 15 minutes scans on the 
PCT4600, as observed in figure 7.25, using the parameters in Table 7.14. The rat was 
sacrificed at two hours post injection, with brain, heart, lungs, liver and blood being 
harvested and activity measured using the Packard-Tricarb well counting system. 


Table 7.14. Acquisition Parameters for Rat MR and PET Imaging. 


MRI Parameters 

PET Parameters 

FOV 

80 mm 

FOV 

160 mm 

Slice Thickness 

3 mm 

Slice Thickness 

9-11.5 mm 

Slice Separation 

8 mm 

Slice Separation 

11.5 mm 

Size 

256x128 

Size 

128x128 

Pixel Size 

0.312 mm 

Pixel Size 

1.25 mm 

Image Type 

T1 weighted 

Reconstruction 

FBP 

TR 

600 ms 

Filters 

Ramp Only 

TE 

20 ms 

Cutoff 

1.8 cm"' 

Elapsed Time 

2:33 min 

Elapsed Time 

15 min 

FI 

85.495 MHz 




Calibration of the PCT 4600 was performed using appropriate ring sources and 
blanks for normalization and attenuation correction, and using the cylindrical phantom 
and line source to obtain system calibration factors. A concentration based calibration 
factor was determined from a region of interest analysis of the uniformly distributed 
cylindrical phantom. Activity concentration was assessed using a 0.5 ml aliquot of the 
phantom measured in a Nal(Tl) well counting system. Activity concentration was 
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verified against the known injected activity and phantom volume. Calibration factors 
were calculated in a manner identical to those for SPRINT. 



Figure 7.25. Rat in Animal Holder During PET Imaging. 


Animal image analysis was performed by drawing regions of interest around and 
within organ boundaries of the reconstructed images. A pixel count rate (cpm/pixel) was 
calculated from the average pixel counts found within a region and the total acquire time. 
The activity concentration was then assessed by dividing by the concentration based 
calibration factor. Alternatively, total organ activity was assessed using the line source 
calibration factor and total region count rates for ROIs that encompassed the organ, as 
was done for SPECT image analysis. 
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Results and Discussion 

The central slice of the line source is shown in figure 126A, with associated 
horizontal and vertical profiles in figure 7.26B. The reconstructed image of the 
cylindrical phantom is shown in figure 7.26C, with its profiles in 7.26D. Analysis of the 
line source indicated a system resolution of 7.6 mm FWHM and 14.79 mm FWTM. This 
is two to three mm wider than that achieved by SPRINT using a 3.175 mm slit aperture 
and aperture penetration correction. The lineal activity concentration in the line source 
was used to obtain an absolute system sensitivity from the total counts obtained in a given 
slice during a 15 minute acquisition. Sensitivity was measured at 10600±270 cpm/pCi 
(4750±120 cps/MBq), approximately 160 times that achieved with SPRINT equipped 
with the high energy parallel collimator, a 3.125 mm slit, and use of aperture penetration 
correction. It should be noted that PET sensitivity is not a function of resolution, as it is 
in SPECT based systems. 

A direct comparison of the noise characteristics in reconstructed images between 
PET and SPRINT is difficult because of the strong dependence on the reconstruction 
algorithm and filters used. However, the dramatic difference in sensitivity between these 
systems indicates that for equal activities and imaging times, PET images will have 
significantly better noise characteristics than SPRINT images. Conversely, the sensitivity 
difference indicates that PET imaging requires 160 times less activity than SPRINT to 
obtain an equivalent quality image when using the same reconstruction method. As an 
example, the standard deviation in the PET image of a uniform cylinder was 
approximately 3.8%, reconstructed with an average of 2600 count per pixel (5.8x10^ total 
counts), FBP and a pure ramp filter. The image was composed of 1.56 mm pixels and 
the activity concentration in the phantom was 1.85 pCi/ml. A SPRINT image of the 
uniform 2-D phantom region had a standard deviation of 24%, and was reconstructed 
with an average of 1800 counts per pixel (8x10® total image counts), using SAGE with 
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P=2'^^. The image was composed of 0.77 mm^ pixels and the background activity 
concentration was 135 pCi/ml. We observe that despite using a more noise tolerant 
algorithm (SAGE versus EBP) there is considerably more noise in the SPRINT image. 
This is partially do to the poorer count statistics (fewer counts in smaller pixels), and 
variance added by the penetration correction method employed for SPRINT. 



Profile through Line Source Image 
7.63 +0.16 mm FWHM 
14.69 + 0.28 mm FWIM 



Horizontal-Vertical 

. ^ B. 



Profile through 6.6 cm Cylinder of Uniform 
Activity 



-Horizontal-Vertical 


D. 


Figure 7.26. Image and Profiles for Uniform Cylinder Source and Line Sources 
Imaged by Positron Emission Tomography. 


The reconstructed animal images are shown in figure 7.27, covering the region 
from the top of the animals head to the upper thoracic region. Matched MR images are 
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shown above the PET image for anatomical reference. The first PET image shows 
Harderian gland uptake and regional brain uptake. The second slice shows only brain 
uptake. The fourth image covers at least a portion of the cardiac region. Comparing these 
images to those in figure 7.21, the PET images have considerably worse resolution, both 
axially and transaxially. However, the EBP images are of generally low noise, with the 
area of the Harderian gland, brain and heart well defined and of high contrast compared to 
surrounding tissue. 

Table 7.15 shows the relatively poor results obtained when attempts were made to 
quantify image data. Regions of interest were scribed within and around the total brain 
(image B) and cardiac region (image D). Biodistributions at two hours post injection 
were 2.16 %dose/organ for the brain and 4.87 %dose/organ for the heart, in reasonable 
agreement with the animals used for SPECT imaging. These levels correspond to brain 
and heart activities at the time of imaging of 17 and 48 pCi, respectively. For both these 
organs, there is a considerable underestimate of both tissues concentration and total organ 
activity. For the brain, concentrations were a factor of two less than that measured 
invitro, and total activity was a factor of four less. For the heart, there is a consistent 
factor of four difference for both concentration and total activity between the invitro and 
invivo measurements. Since the calibration process was accomplished using two 
different methods, it is unlikely that errors in calibration account for these differences. 
Examination of the corresponding MR images indicates that for both brain and heart, 
there may be serious inaccuracies in quantification because of poor axial resolution, and 
limited field of view of this system. For the brain, the large slice thickness may cause 
scribed regions to include areas of low uptake outside the brain. This effect, compounded 
by the larger pixel size and poorer in-plane resolution as compared to SPRINT cause 
ROIs to measure an activity concentration considerably lower than what is actually 
present. Similarly, the determination of total activity requires summing the total count 
rate for all region in which the brain is observed. An underestimate of brain activity is 
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likely due to excluding the first slice, which had considerably activity contributed by the 
Harderian glands. If this slice is included, the brain activity is overestimated by several 
hundred percent. For the heart, a significant portion of the heart may have been outside 
the sensitive region of the system, resulting in the negative biases observed for 
concentration and total activity. This argument is supported by examining the MR and 
PET images of figure 7.27. The last available PET image appears to be above the 
primary cardiac region as noted by the MR images. 


Table 7.15. Quantitative Results Obtained from PET Image Analysis. 


Activity 

Concentration 

Small region count rate 
(brain and heart) 

Calculated 

Concentration 

Estimated^ 

Concentration 

Error (%) 

Assessment 


(ROI cpm/pixel/CF) 



Brain 

cpm/pixel ±cpm/pixel 

pCi/cc ±pCi/cc 

pCi/cc 


Acquis. 1 

65560 ±5540 

6.86 ±0.60 

8.74 

-21.5% 

Acquis. 2 

57693 ±4718 

6.04 ±0.51 

7.90 

-23.6% 

Heart 

cpm/pixel ±cpm/pixel 

pCi/cc ±pCi/cc 

pCi/cc 


Acquis. 1 

97933 ±10467 

10.25 ±0.17 

46.5 

-77.9% 

Acquis. 2 

94667 ±12867 

9.91 ±0.19 

42.0 

-76.4% 

Calibration Factor: 1160000 ±41200 cpm/pCi 


-- ^ 

Total Activity Large region count rate' 

Assessment 

Calculated Activity 
(ROI cpm.CF) 

Estimated^ 

Activity 

Error (%) 

Brain 

cpm 

pCi 

pCi 


Acquis. 1 

5383333 

4.65 

18.0 

-74.1% 

Acquis. 2 

4392000 

3.80 

16.3 

-76.6% 

Heart 

cpm 

pCi 

pCi 


Acquis. 1 

10200000 

8.82 

40.5 

-78.3% 

Acquis. 2 

10486667 

9.07 

36.7 

-75.3% 

Calibration Factor (CF): 9550 ±212 cpm/pixel/pCi/ml 


* ROI analysis based on small to moderate size regions within the organ boundaries to 
provide an estimate of organ concentration. 

^ ROI analysis based on large regions which completely encompass the observed organ. 

^ Estimated organ activity concentration or total activity based on invitro measured organ 
concentration at time of sacrifice, decay corrected to start time of imaging. 
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Although there is poor accuracy in these PET images, the PET images show much 
lower noise when using EBP with a pure ramp filter, than similar SPRINT images 
reconstructed with SAGE (P=2'^^), which devalues noisy data. For PET, the animal was 
injected with three to ten times less activity then used in SPRINT imaging, and image 
acquisition times were almost half (15 minutes for PET compared to 27 minutes for 
SPRINT). The lower noise is again due principally to the much higher sensitivity of PET. 

Summary 

Qualitative and quantitative analysis of SPRINT images indicated that good 
accuracy and fair precision could be obtained of large objects in simple phantoms and of 
large animal organs that have significant radiotracer uptakes. Accuracy was dependent on 
numerous corrections to the raw projection data, as well as the methods employed for 
image reconstruction, image analysis and system calibration. 

The very poor sensitivity of SPRINT for 511 keV photons, and variance 
introduced by aperture penetration correction caused all SPECT images to be noisy. 
Minimum standard deviations of 20% or greater were observed in the uniform 
background regions of very high count phantom studies, increasing up to several hundred 
percent for low activity cold regions. Reconstruction using SAGE provided considerable 
improvement in reducing noise when compared to filtered backprojection reconstruction 
where minimum filtering was used. A result of the iterative algorithm weighting the 
projection data by its statistical significance. Low noise SPRINT images were only 
acquired by using very high activities (30 to 40 mCi of ’^F) and imaging for long periods, 
particularly in comparison to PET. Unfortunately, high activities also introduced bias in 
quantified results due to uncorrected deadtime effects. 

An important disadvantage posed by using long acquisition times (in addition to 
the requirement for decay correction) is loss in temporal resolution of animal 
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biodistributions. Since imaging times are long compared to the short uptake and 
biodistribution of most PET agents, is impossible to observe specific biodistribution 
patterns. The poor sensitivity of SPRINT thus precludes the assessment of labeled 
compounds, since images could not be collected fast enough (<5 min) and still provide 
images of high signal to noise. 

Despite SPRINT’S poor sensitivity, the ability to achieve high resolution images 
that can provide good quantitative results has been demonstrated. This resolution is 
achieved through a magnifying slit aperture geometry, high resolution collimation and 
correcting acquired projection data for aperture penetration by 511 keV photons. The 
resulting in-plane resolution was as low as 4.4 mm for SAGE reconstructed images, 2 to 
3 mm less than that obtained with a commercial 2-D SPECT system. SPRINT also 
demonstrated good axial resolution provided by the specially designed high energy slice 
collimator. These developments indicate that quantitative high energy animal imaging 
can be performed, but only if the object to be assessed has sufficient volume and uptake, 
and poor sensitivity can be tolerated. 
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CHAPTER VIII 


SUMMARY AND CONCLUSIONS 

This study has examined two methods of hiodistribution assessment in rats for 
ligands labeled with positron emitting radionuclides. The first objective of this work was 
to explore ways to simplify conventional invitro methods of biodistribution assessment. 
For the case of ‘’C labeled tracers, evidence has been presented that indicates examining 
the hiodistribution within the first 15 minutes post injection is sufficient to characterize 
the overall hiodistribution for solid organs. A reduced sacrifice protocol based on two 
measurements during this period provides estimates of cumulated activity that are either 
equivalent or conservative compared to those obtained using multiple sacrifice intervals. 
For the gallbladder and urinary bladder, a correction is required to the assessed cumulated 
activity to compensate for delayed organ uptake. For the ten radiotracers examined, the 
maximum activity which could be injected and not exceed a five rad absorbed dose limit 
to a given organ was found to range from 410 MBq (11 mCi) for [*'C]FNZPAM to 3441 
MBq (93mCi)for[^^C]MTBZ. 

The use of 3D PET imaging systems provides greater relevance to this proposed 
method. These new systems have sensitivities that are six to ten times that achievable 
with 2D systems, allowing preliminary studies of new radiotracers to be conducted with 
significantly lower dosages. Thus, even the more conservative measures of absorbed 
dose assessed with this method result in allowable injected doses which can readily be 
used for clinical evaluation of new compounds. 

The second objective of this work was to examine the capabilities and limitations 
of using SPRINT to assess biodistributions invivo of positron emitting radiotracers. It 
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has been demonstrated that there are severe limitations in the use of SPRINT for imaging 
511 keV photons, especially when using high resolution collimation. Specifically, system 
sensitivity at high energies is simply inadequate to acquire projection data with the 
necessary time resolution to trace organ specific time-activity profiles of short half-life 
emitters. 

Although sensitivity is low for positron imaging, it has been established that 
SPRINT can provide high quality images if activity concentrations on the order of 100 
pCi per ml can be attained and long acquisition times are employed. The system can 
produce images of 4 to 5 mm resolution axially and in-plane, with modest noise in 
uniform regions (<20% standard deviation). Rodent imaging with injected dosages of 3 
to 10 mCi [**F]fluorodeoxyglucose provided images with good contrast and resolution of 
the brain and heart. The capacity to image 511 keV photons at high resolution is 
attributed to the magnifying slit-aperture geometry of SPRINT, the development of a 
simple method to overcome aperture penetration by high energy photons, and a precision 
high energy slice collimator. The compromise made in achieving high system resolution 
was a loss in system sensitivity. For the slice collimator, an effective open area of 66% 
resulted in a point source sensitivity of only 127 cpm/pCi, when a 3.125 mm slit opening 
was used. This already low sensitivity is reduced by more than half when implementing 
the blocked slit method for correcting aperture penetration. 

Very low system sensitivity in turn affects the ability to successfully estimate 
organ concentration from image data. Low noise images were only observed in very high 
count images which were negatively biased due to uncompensated deadtime. 
Alternatively, low count images were very noisy (>50% SD in uniform regions), 
providing quantitative estimates with low precision. Despite these issues, “hot” phantom 
objects, and the rodent brain and heart were assessed with errors ranging from +12% to 
- 8 %. 
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Comparison of the sensitivity of SPRINT (configured with a 3.175 mm slit and 
implementation of penetration correction) with that of the PCT 4600 PET system 
demonstrated a 160 fold reduction in sensitivity, chiefly due to the use of electronic 
collimation and the high intrinsic efficiency of the BGO detector elements. This 
difference approaches a factor of 1600 with septaless 3D PET systems. The resolution of 
clinical PET systems is generally worse than that achieved with SPRINT, on the order of 
6 to 8 mm FWHM. (Note that if SPRINT were designed for 8x8 mm resolution, 
sensitivity would be improved by a factor of 2.5). However, new research in the area of 
PET imaging has reported the development of high resolution, small field-of-view PET 
systems designed specifically for small animals * ^ 3 4 5 6 jjjggg systems report in-plane 
resolutions ranging from 1.8 mm to 4.6 mm FWHM (3.1 mm to 5.4 mm axially) and 
sensitivities ranging from 3300 cps/|j,Ci/ml (198,000 cpm/pCi/ml) to 20700 kcps/pCi/ml 
(1.24x10^ cpm/pCi/ml), which are comparable to 2D clinical systems. Further, these 
instruments have already been applied to the measurement of glucose metabolism and “C 
biodistribution assessment in the rat ’ * The greater sensitivity and higher resolution 
afforded by these small animal PET systems make these instruments superior for invivo 
biodistribution assessments of positron emitters compared to SPRINT and other SPECT 
based imaging systems. 
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Recommendations for Future Work 


It is evident that increasing system sensitivity is highly desirable if further 
qualitative (and quantitative) improvements are to be achieved with SPRINT imaging at 
high energies. Thicker scintillators with higher stopping power than Nal(Tl) are an 
effective solution if resolution losses can be minimized. Current research suggests that 
lutetium oxyorthosilicate (LSO) may provide considerable improvements in detection 
efficiency, while still providing a very fast decay time and good light output’*^. 
Characteristics of this scintillator in comparison to Nal are provided in Table 8.1. The 
use of a 25 mm thick LSO scintillator can provide a 50% photoelectric detection 
efficiency at 511 keV, as compared to 7% for 13 mm ofNal(Tl). Disadvantages if this 
material is that it is presently very expensive and it is mildly radioactive due to a naturally 
occurring contaminant. Resolution losses in thick scintillators can potentially be limited 
by use of digital light pipes", energy specific position weights, and implementation of 
surface altering technologies which would shape the light distribution within the 
scintillator . 

Table 8.1. Characteristics of LSO as a High Efficiency Scintillator. 




Photoelect. 

Photoelect. 

Relative 

Decay 




Density 

Attenuat. 

Interaction 

Scintillation Time 

Radioactive? Hygro.? 


(g/cm^) 

Coefficient Probability 

Efficiency 





(cm’) 

(25 mm) 

% 

(ns) 



Nal 

3.67 

0.0585 

14% 

100 

230 

No 

Yes 

LSO 

7.4 

0.27 

50% 

75 

40 

Yes 

No 


If SPRINT detectors were made 100% efficient for 511 keV photons, overall 
system sensitivity would increase by an approximate factor of 10*. This is still 
considerably less than PET systems because of the different sensitivities provided by 


* This is based on the current 12.7 mm Nal SPRINT modules having an 11% detection probability assessed 
from Monte Carlo simulations (7% photoelectric effect detection probability). 
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electronic versus mechanical collimation. To improve SPECT sensitivity further, either 
the means of mechanical collimation must be optimized, or a method for electronic 
SPECT collimation must be employed. The latter serves as the basic operating premise 
for the Compton aperture camera which uses electronic coincidence between a primary 
scattering detector and a secondary scintillation detector (e.g. SPRINT) for position 
estimation'^ Theoretical performance of these systems suggests high spatial resolution 
and good sensitivity can be achieved for high energy photons. Currently, these systems 
are being aggressively investigated for their application in high energy medical and 
industrial imaging. 

Reengineering the mechanical collimation of SPRINT through use of a six pin¬ 
hole aperture (instead of slits) and removing the slice collimator could also improve 
sensitivity. Installing a pin-hole aperture would provide three-to-one axial magnification 
increasing sensitivity for a reduced axial field of view. Other researchers have 
investigated the use of pin-hole SPECT systems and cone-beam EBP algorithms to 
conduct high resolution low-energy animal imaging of ^^"’Tc labeled compounds 
19 20 performance for two of these systems as compared to SPRINT are listed in 

Table 8.2. In general, resolutions and sensitivities between systems are similar. 
sprint’s higher sensitivity originates from its full ring design and six apertures, offset 
by the reduced sensitivity of a longer source to aperture distance** and its present fan- 
beam geometry. 

Aperture penetration in such geometries still requires correction if penetration 
blurring of the PSRF is to be minimized. This could be accomplished by blocking 
alternate pinholes and correcting the projection data for block penetration and penumbra 
using a two-dimensional weighting function. Investigating an optimal block design to 

** Pinhole collimator efficiency decreases inversely with the square of the aperture to 
source distance. 
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minimize penumbral effects and developing an iterative solution for the penetration 
corrected image may provide additional qualitative improvements when using this 
“block” technique. 

The optimal solution to aperture penetration is to accurately model the penetration 
effect and incorporate it into the system response matrix. A statistical iterative 
reconstruction technique, such as SAGE, can then be used to provide images that should 
have the best achievable sensitivity and resolution. The usefulness of this method is 
dependent on the accuracy by which penetration can be modeled for the entire system, 
and the availability of computer platforms capable of handling the resultant non-sparse 
weights matrix. 

Successful development of an aperture penetration correction algorithm suitable 
for pin-hole geometries would allow high resolution small field-of-view pin-hole SPECT 
to be performed using clinical imaging systems. Since low energy pin-hole SPECT has 
already been investigated using a clinical four head gamma camera^' application of the 
developed correction method to clinical multi-head systems should be relatively straight¬ 
forward. 

Table 8.2. Comparison of Pinhole Collimator Resolution and Sensitivities for ’’™Tc 

for Several Imaging Systems. 

Weber etal.^‘’>* Jaszcak^'*>* SPRINT^ 

(Source to Aper. 4.0 cm) (Source to Aper. 4.0 cm) (Source to Aper. 6.5 cm) 


(FOV 

-3.0 cm) 


(FOV=6.0 cm) 


(FOV=9 cm) 



Resolution 

Sensitivity 

d 

Resolution 

Sensitivity 

d 

Resolution 

Sensitivity 

(mm) 

(mm) 

(cps/pCi) 

(mm) 

(mm) 

(cps/pCi) 

(mm) 

(mm) 

(cps/pCi) 




0.6 

1.5 

0.9 




1.0 


1.1 

1.2 

1.9 

2.6 




2.0 

2.8 

4.0 

2.0 

2.8 

5.7 

2.0 

2.9 

11 

3.3 


10.9 




3.1 

4.4 

17 


Single head gamma cameras 
^ Full ring tomograph 
^ Aperture opening 
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Improvements in SPRINT sensitivity would allow it’s application to other high 
energy imaging applications, including qualitative and quantitative tumor imaging. 

Tumor metabolism could be assessed using [ FJfluorodeoxyglucose, and the diagnostic 
and therapeutic efficacies of iodine labeled compounds such as [*^^I]MIBG and ’^*1 
labeled monoclonal antibodies could be examined. The importance of invivo assessment 
of tumor radioactivity originates from the unique nature of xenografted tumors in 
laboratory animals. In contrast to the typically small variance in biodistributions seen 
between organs in normal animals (reference Chapter II), the biodistribution in tumors are 
extremely variant due to the unique physiology and anatomy of a given tumor. 
Consequently, there is a clear advantage in determining the time activity profile of a 
radiopharmaceutical or radiotracer within the same animal. Further, invivo assessment 
provides long term tracking of tumor progression and can allow comparison of the 
pharmacological affect of different compounds in the same subject. 

Since SPRINT acquisition software now allows multi-window image acquisitions, 
the potential does exist to monitor within the same animal the biodistributions of FDG 
and an iodine labeled radiotracer simultaneously. Figure 8.1 demonstrates the difference 
in spectra resulting from I and F imaged separately and jointly with the high energy 
parallel collimator in place. Both peaks have a resolution of approximately 13 percent 
when imaged separately, with the iodine peak showing some degradation in resolution 
when imaged with '^F. For equal source activities, the annihilation peak area is 
approximately 1.3 times the iodine peak (this is explained in Chapter V). Most notably 
when both sources are imaged, approximately 37% of the iodine peak is attributable to 
scattered annihilation radiation, whereas 7% of the annihilation peak is attributable to 
scattered high energy iodine energies. Dual energy imaging would therefore require 
integrating methods to correct for photopeak contamination from opposing nuclides. 
Methods include spectrum stripping and the use of injected activities of iodine and 
fluorine which balance the level of contamination in the adjacent photopeaks. 
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Comparative Spectra from Iodine-131 and Fluorine-18, 
42 microCuries each 

25 - 



-18FOnly.1311 Only-Joint Spectra 


Figure 8.1. Module Spectra of Simultaneously Acquired and Sources. 
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Conclusions 

The major focus of this work has investigated the application of SPRINT to 
quantitative animal imaging for rats injected with positron emitting radioligands. It has 
been shown that measuring the biodistribution of very short half-life compounds is not 
possible due to the limited sensitivity of SPRINT and the constraints made on injected 
activity in the animal. Further, it has been shown that the general application of SPECT 
to 511 keV imaging is extremely unfavorable compared to the use of coincidence 
imaging. However, good qualitative and quantitative results have been observed when 
using high activities and long acquisition times for SPECT imaging of F. Implementing 
technologies to improve system sensitivity will allow more rapid imaging and provide 
improved image quality. A more thorough evaluation of such factors as deadtime and 
scatter should allow higher quantitative accuracy in image assessment and allow SPRINT 
be a more effective tool in assessing animal biodistributions of high-energy emitting 


radiotracers. 
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